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This article reviews recent theoretical and experimental advances in the fundamental understand- 
ing and active control of quantum fluids of light in nonlinear optical systems. In presence of effec- 
tive photon-photon interactions induced by the optical nonlinearity of the medium, a many-photon 
system can behave collectively as a quantum fluid with a number of novel features stemming from 
its intrinsically non-equilibrium nature. We present a rich variety of photon hydrodynamical ef- 
fects that have been recently observed, from the superfluid flow around a defect at low speeds, 
to the appearance of a Mach-Cherenkov cone in a supersonic flow, to the hydrodynamic forma- 
tion of topological excitations such as quantized vortices and dark solitons at the surface of large 
impenetrable obstacles. While our review is mostly focused on a class of semiconductor systems 
that have been extensively studied in recent years (namely planar semiconductor microcavities in 
the strong light-matter coupling regime having cavity polaritons as elementary excitations), the 
very concept of quantum fluids of light applies to a broad spectrum of systems, ranging from bulk 
nonlinear crystals, to atomic clouds embedded in optical fibers and cavities, to photonic crystal 
cavities, to superconducting quantum circuits based on Josephson junctions. The conclusive part 
of our article is devoted to a review of the exciting perspectives to achieve strongly correlated 
photon gases. In particular, we present different mechanisms to obtain efficient photon blockade, 
we discuss the novel quantum phases that are expected to appear in arrays of strongly nonlinear 
cavities, and we point out the rich phenomenology offered by the implementation of artificial gauge 
fields for photons. 
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I. HISTORICAL INTRODUCTION 

In the last decades, the study of the physics of quantum 
fluids has attracted a tremendous interest in a variety of 
different many particle systems, ranging from liquid He- 
lium ( ,Leggetti 2004 Pines and Nozieres,|1998}, electrons 



in solid-state materials (Mahan 1990 Schrieffer [1999 
Tinkh am','200 4t, trapped gases of ultracold atoms fBlochJ 



et al.^ 2008; Dalfovo et al.\ |1999| [Giorgini et ai, ,2008^, 



quark- gluon plasma in colliders 
et al. 



2005), nuclei (Peter Ring [2004). When the ther- 



Satz et al. 2010 Yagi 



mal de Broglie wavelength becomes comparable or larger 
than the average interparticle spacing, the Bose vs. Fermi 
statistics of the constituent particles starts playing a cru- 
cial role in determining the properties of the fluid: in non- 
interacting Fermi gases, the Pauli principle is responsible 
for the rigidity of the Fermi sphere and the appearance of 
a Fermi pressure down to zero temperature, while in Bose 
gases a macroscopic fraction of the particles accumulate 
into the lowest energy single-particle state, the so-called 
Bose- Einstein condensate (BEG). The situation is even 
richer when quantum degeneracy combines with signif- 
icant inter-particle interactions to produce a variety of 
spectac ular effects such as superconductivity an d super- 
fluidity (|Leggett[ |1999[ [Tilley and Tilley[ |1994[) and the 



Historically, most of the theoretical and experimen- 
tal activities in this field of many-body physics have ad- 
dressed systems of material particles such as atoms, elec- 
trons, nucleons, or quarks. However, in the last decades, 
a growing community of researchers has started wonder- 
ing whether in suitable circumstances light can be consid- 
ered as a fluid composed of a large number of corpuscular 
photons with sizable photon-photon interactions. Even if 
this point of view is perfectly legitimate within the wave- 
particle duality in quantum mechanics, it is somehow at 
odd with our intuitive picture of light: the historical de- 
velopment of our understandings of matter and light have 
in fact followed very different paths. 

The idea of matter being formed by a huge number 
of elementary corpuscles that combine in different ways 
to form the variety of existing materials dates back to 
the ancient age with Demokritos' atomistic hypothesis, 
while the wavy nature of particles was put forward only 
in 1924 by de Broglie and experimentally demonstrated 
by Davisson and Germer in 1927. On the other hand, 
the long-standing debate between Newton's corpuscular 
and Huygens' undulatory theories of light appeared to be 
solved in the early nineteenth century with the observa- 
tion of fringes in Young's double slit experiment and of 
the remarkable Arago's white spot in the shadow of a cir- 
cular object. With the microscopic support of Maxwell's 
theory of electromagnetism, the undulatory thory was 
able to explain most experimental observations until the 
beginning of the twentieth century when the corpuscular 
concept of a photon as a discrete quantum of light was 
revived by Einstein's theory of the photo-electric effect. 
Within the wave-particle duality, our standard interpre- 
tation of light then consists of a dual wave/particle beam 
that is emitted by the source and then freely propagates 
through optical devices until it is absorbed. 

While this intuitive picture of light is perfectly suffi- 
cient to describe most cases of interest, still it is missing 
a crucial element, namely the possibility of frequent colli- 
sions between photons that allow for collective fluid-like 
behaviors in the many photon system. While photon- 
photon interactions have been predicted to occur even 
in vacuum via virtual excitation of electron-positron 
pairs (Heisenberg and Euler 1936), the cross section for 



such a process is so small that it can hardly be expected 
to play any role in realistic optical systems. On the other 
hand, the nonlinear polarization of nonlinear optical me- 
dia is able to mediate significant interactions between 
photons ( Boyd, ,2008; jButcher and Gotter , |2008[): upon 



elimination of the matter degrees of freedom, third-order 
X^^'^ nonlinearities correspond in the language of Feyn- 
man diagrams to four-legged vertices describing, among 
other, binary collisions between a pair of photons. 

Among the many different configurations that have 
been studied in the last few decades for nonlinear op- 
tical applications, systems in the so-called strong light- 
matter coupling regime have turned out to be particularly 
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promising in order to obtain the relatively strong nonlin- 
ear interactions that are necessary for collective behavior. 
In this regime, the photon is strongly mixed with mat- 
ter degrees of freedom, which gives rise to a new mixed 
quasi-particle, the po^ariton (Hopfield 19581. Pictorially, 



Staliunas 1993), where the time-evolution of the coher- 



the polariton can be seen as a photon dressed by a mat- 
ter excitation: a reinforced optical nonlinearity then ap- 
pears thanks to the relatively strong interactions between 
matter excitations. This strong coupling regime can be 
achieved in a number of material systems, from atomic 



gases ( Berman et al. 1993 Fleischhauer et al. 2005 Rai 



mond et al. 2001 ) to semiconducting solid state media 
both in bulk (IKlingsh irnH2007| |Yu and Cardona[ |20 05D 
and in cavity ( DcvcaudJ' 2007 Weisbuch et al.. 1992 



ge- 



ent electromagnetic field in a laser cavity with large Fres- 
nel number was reformulated in terms of hydrodynamic 
equations for the many photon system analogous to the 
Gross-Pitaevskii equation for the superfluid order param- 
eter. The local light intensity corresponds indeed to the 
photon density and the spatial gradient of its phase to 
the local current; the collective behavior originates from 
the effective photon-photon interactions stemming from 
the nonlinear refractive index of the medium as well from 
saturation of gain. In the following years, the transverse 
dynamics of the electromagnetic field in cavity devices 
has attracted a lot of attention, in particular the phe- 
nomena related to the spontaneous formation of trans- 



ometries, to circuit-QED systems based on superconduct- verse patterns (Denz et al. 2010 Staliunas and Morcillo 



ing Josephson junctions (Schoelkopf and Girvin 



2008 



2003) and to the generation and control of dissipative 



You and Nori 2011). In the following of the review, we cavity solitons (Ackemann et al. 2009) 



will consider both photon and polariton excitations, de- 
pending on the non-resonant or resonant character of the 
electronic excitation dressing the photon within the ma- 
terial medium. 

To create a stable luminous fluid, it is also crucial to 
give a finite effective mass to the photon. A simplest 
strategy to this purpose involves a spatial confinement of 
the photon by metallic and/or dielectric planar mirrors. 
In a planar geometry with a dielectric medium of refrac- 
tive index hq and thickness £z enclosed within a pair of 
metallic mirrors, the photon motion along the perpen- 
dicular z direction is quantized as Qz — ttA/ /£z , M being 
a positive integer. For each longitudinal mode, the fre- 
quency dispersion as a function of the in-plane wavevec- 
tor k has the form 



OJc 



no 2.mca^ 



(1) 



where the effective mass Wcav of the photon and the cut- 
off frequency oj°^^ are related by the relativistic-like ex- 
pression 



hn^qz 



2 • 



(2) 



Using suitable values of the effective mass mcav and the 
cut-off frequency 0J°g_.^, extracted from microscopic calcu- 
lations (Savona 19991, the generic form ([T]) of the disper- 
sion can be extended to the case of dielectric mirrors. In 
the presence of some electronic excitation resonant with 
the cavity mode, the elementary excitations of the cavity 
have a polaritonic character with a peculiar dispersion 
law that reflects their hybrid light-matter nature. An 
example of such dispersion is shown in the central panel 
of Figjl] in spite of the complex light-matter interaction 
dynamics, the bottom of the lower polariton branch is 
still well approximated by a parabolic dispersion of the 
form ([T]) with an effective mass niLP and the cut-off fre- 
quency w£p. 

Historically, the first mention of the concept of photon 



In this pioneering literature on hydrodynamics of the 
photon fluid, a special attention was paid to the phase 
singularities of the photon field, that were immediately 
interpreted as quantized vortices (Coullet et al. 19891. 



After their first experimental observation in bulk nonlin- 
ear crystals ( Swartzlander and Law 1992), most of the 



following literature addressed the physics of optical vor- 
tices in the context of the transverse dynamics of laser 
or photorefractive oscillators: the drift of vortices un- 
der the hydrodynamical effect of the Magn us force due 
to buoyancy was experimentally studied in (Vaupel and 



Weiss 1995). The hydrodynamical shedding of vortices 



in a moving photon fluid hitting a large defect was first 
predicted in the optical context in (IStahunasl 119931) and 



a pioneering attempt of experimental investigation of this 
crucial effect of superfluid hydrodynamics was reported 



soon after in (Vaupel et al. 1996) 



From a different perspective, the close analogy between 
a laser threshold and a second order phase transition was 



recognized as early as in (DeGiorgio and Scully 1970 



Graham and Haken 19701: above the laser threshold. 



the electromagnetic field acquires a well defined phase 
by spontaneously breaking a U{1) symmetry as it hap- 
pens to the matter Bose field in a Bose-Einstein con- 



densate of material particles ( Gunton and Buckingham 



1968 Huang 1987). This interpretation of lasing as the 



result of a kind of Bose-Einstein condensation of photons 
is clearest in spatially extended devices such as verti- 
cal cavity surface emitting lasers (VCSELs), where the 
onset of a coherent laser emission is associated to the ap- 
pearance of long-range spatial coherence along the cavity 
plane according to the Penrose-Onsager criterion for off- 
diagonal long-range order (Huang 1987 Pitaevskii and 
Stringari[ [2004| 



lim 



Et(r)E(r'))^0. 



(3) 



Of course, the non-equilibrium nature of the laser device 
introduces crucial differences with respect to standard 
equilibrium BEG as explained in statistical mechanics 



fluid dates back to the work of (Brambilla et al. 1991 textbooks (Huang 19871: the steady state of the laser 
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device is in fact not determined by a thermal equilibrium 
condition, but rather follows from a dynamical balance 
between the pumping and losses (Haken 1975[ ). As we 
shall see in the following of the review, this feature is 
responsible for a number of new effects. 

From the experimental point of view, the close link be- 
tween BEC and spontaneous coherence effects in optical 
systems has been fully recognized only in the last decade, 
starting with the literature on the so-called BEC of 
exciton-polaritons in semiconductor microcavities (Baas 



eFal] [2006, Baumberg et aLj |2000| [Kasprzak et ai:||2006 



Stevenson et al. 2000[ ). In~this context, questions related 



to the analogies and differences between laser operation 
and photon/polariton BEC have attracted a strong inter- 
est from the community, with a special attention to ther- 
malization issues. An experimental observation of micro- 
cavity polaritons coherently accumulating in the lowest 
energy states of a harmonic trap potential according to a 
Bose distribution was reported in (Balili et al. 20071: 



thermalization of the polariton gas was attributed to 
polariton-polariton collisions within the gas. A similar, 
apparently thermalized photon distribution was however 



observed in (Bajoni et al. 2007) also in a weak-coupling 



regime where photon-photon interactions are very weak, 
which raises fundamental questions about the nature of 
fluctuations on top of a photon/polariton condensate. 
These observations are to be contrasted with the strongly 
non-equilibrium regimes of laser operation observed in 
a vertical cavity surface emittin g laser (VCSEL) device 
in ( Scheuer and Orenstein 1999 1 : thermalization into the 



lowest state is completely ineffective and the condensate 
mode displays a complex structure with array of vortices. 

In the last years, the quest for condensation effects in 
photon gases has successfully explored a few other inter- 
esting avenues. An early mention of the possibility of 
a gas of bare photons thermalizing to a Bose condensed 
state via co llisions media ted by the optical nonlinearity 
is found in ( Navez , 2003) . BEC in a thermalized gas of 



photons was experimentally observed in (Klaers et al. 



2010a) using a macroscopic optical cavity containing a 
dye solution: because of the rapid decoherence time of 
the dye molecules, photons are only weakly coupled to 
the electronic excitations and their thermalization is be- 
lieved to occur via repeated absorption-emission cycles, 
which determine the temperature and the chemical po- 
tential of the gas in a grand-canonical picture (Klaers 



et al. 2012 2010b I. A kinetic condensation of purely 



et al. 2008) and in actively- mode- locked lasers (Weill 
erZ[[20T0a|bl ). 



Simultaneously to these studies on Bose-Einstein con- 
densation and spontaneous coherence effects, a revived 
interest has been devoted also to the hydrodynamic prop- 
erties of the photon gas. The concept of the Bogoliubov 
dispersion of elementary excitations on top of a photon 
condensate was first investigated for a planar cavity ge- 



ometry in the pioneering works (Chiao and Boyce 1999 



[Tanzini and Sorella 1999): thanks to the spatial con 



finement, photons have a massive dispersion of the form 
([T]) and the photon-photon interactions responsible for 
the collective behavior of the fluid are provided by the 
X*^"^) nonlinearity of the cavity medium. None of these 
works however addressed the crucial role of dissipation 
in the physics of photon fluids: fir st steps in this direc- 
tion appeared in (Bolda et al. 2001 ) where the nucleation 



of vortices in a moving photon fluid past an impenetra- 
ble, cylindrical defect in a planar cavity was theoretically 
investigated. Differently from previous works on oscilla- 
tors (Staliunas 1993), a coherent pumping was proposed 



as a way to create the photon superfluid in the cavity: 
this resulted in a different form of the Gross-Pitaevskii 
equation as originally studied in ( [Lugiato and Lefever[ 



1987) and, more importantly, in the need to switch off 



the coherent pump to unlock the condensate phase be- 
fore vortices can appear. 

A completely different approach to the study of super- 
fluidity properties of light involved the paraxial propa- 
gation of a light beam through a bulk nonlinear crystal, 
which can be recast into a superfluid hydrodynamic form 
under the replacement of the time coordinate with the 
longitudinal coordinate along the propagation direction. 
After the pioneering experiment in ( Swartzlander and 



Law 1992 ), many authors have theoretically investigated 



a number of hydrodynamic features in lig ht propagation, 
from stable liquid- like solitonic st ructures ( Josserand and 
RicaHl 997'; 'Mi chinel et aLH2 006), to vorti ces ([Firth and 



Skryabiug ,1997; [Paz-Alonso and Michinelt [2005[) , to th e 



scattering on a defect potential (Khamis et al. 20081, 



to superfluid motion (Leboeuf and Moulieras 2010), to 
dispersive shock waves ( El et aL] "2007l, to the genera- 
tion of optical analogs of acoustic black holes (Fouxon 



et al. 2010). In spite of all this theoretical activity, not 



many nonlinear optics experiments have specifically ad- 
dressed hydrodynamic features of light. Among the few 
exceptions, the experimental study of vortex and soliton 



classical light waves was observed in the remarkable ex- stability and interactions ( Krolikowski et al. . 1998 Ma- 
periment of (Sun et al. 2012): as theoretically discussed maev et al. 1996a|b I, the generation, propagation anc 
in, e.g. , ( Connaughton et al. 2005 1 , turbulent wave mix- interaction of dispersive shock waves in fluids of light ( Jia 



ing by the optical nonlinearity leads to a redistribution of 
energy among the different modes and, eventually, to its 
accumulation into the lowest, condensate mode. A most 
remarkable feature of this experiment is the complete ab- 
sence of quantum features, which emphasizes the funda- 
mentally classical origin of the Bosc-Einstein condensa- 
tion phenomenon. More complex condensation phenom- 



et al. 2007 Wan et al. 2007 ) , the propagation of a one- 



dimen sional fluid through a small barrier potential ( Wan 
|2010) , a study of the generation of vortex pairs 



et al. 



in the wake of an obstacle in a two-dimensional geom- 
etry (Wan et al. 20081, a stu dy of the Raylei gh- Taylor 

and the 



instability in stratified fluids (Jia et al. 



2012 



realization of a trans-sonic flow of light in a 



J aval noz- 

ena have been reported also in disordered lasers ([Conti[ zle configuration ([Elazar et al.\ [2012|. An alternative 
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interesting perspective on many-body physics of photons 
was developed in the pioneering Uterature on quantum 
soUtons in nonhnear optical fiber using a quantum non- 
linear Schrodinger equation as well as Bethe ansatz tech- 



niques ( 


Drummond et al. 


1993j|Kartner and Haus 1993 


Lai and Haus 




1989a|b 


I- 



The research on exciton-polaritons in semiconductor 
microcavities approached the physics of luminous quan- 
tum fluids following a rather different pathway. For many 
years, an intense activity has been devoted to the quest 
for Bose-Einstein condensation phenomena in gases of 
excitons in solid-state materials ( |Grifltin et~ar. 1996): ex- 



citons are neutral electron-hole pairs bound by Coulomb 
interaction, which behave as (composite) bosons. In spite 
of the interesting advances in the direction of exciton 
Bose-Einstein condensation in bulk cuprous oxide and 
chloride 



cuprous 



tov 



bila yer electron systems ( Eisenstein 
and MacDonald 2004), and coupled quantum wells (Bu- 
2007, High et al. 2012), so far none of these re- 



search axes has led to extensive studies of the quantum 
fluid properties of the alleged exciton condensate. The 
situation appears to be similar for what concerns con- 
densates of magnons, i.e. magnetic excitations in solid- 
state materials: Bose-Einstein condensation has been ob- 



served (Demokritov ei aL 2006 Giamarchi ei a/. 2008), 



but no quantum hydrodynamic study has been reported 
yet. 

The situation is very different for exciton-polaritons 
in semiconductor microcavities, that is bosonic quasi- 
particles resulting from the hybridization of the exci- 



ton w ith a planar cavity photon mode ( Weisbuch et al. 
1992). Following the pioneering proposal by Imamoglu 
et al.\ |1996[ researchers have successfully explored the 
physics of Bose-Einstein condensation in these gases of 
exciton-polaritons. Thanks to the much smaller mass of 
polaritons, several orders of magnitude smaller than the 
exciton mass, this system can display Bose degeneracy at 
much higher temperatures and/or lower densities. 

Historically, the first configuration where spontaneous 
coherence was observed in a polariton system was based 
on a coherent pumping of the cavity at a finite angle, 
close to the inflection point of the lower polariton dis- 
persion. As experimentally demonstrated in (Baumberg 



et al. 2000 Stevenson et al. 2000), above a threshold 



valu e of the pump intensity a sort of paramet ric oscilla- 
tion (■CmtreTari 120001 120011 IWhittakerl 120011) occurs in 



the planar microcavity and the parametric luminescence 
on the signal and idler modes acquires a long-range co- 
herence in both time and space (Baas et al. 2006). As 



theoretically discussed in (Carusotto and Ciutil 20051, 



the onset of parametric oscillation in these spatially ex- 
tended planar cavity devices can be interpreted as an 
example of non-equilibrium Bose-Einstein condensation: 
the coherence of the signal and idler is not directly in- 
herited from the pump, but appears via the spontaneous 
breaking of a C/(l) phase symmetry. 

The quest for Bose-Einstein condensation in a thermal- 
ized polariton gas under incoherent pumping required a 
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FIG. 1 Figure from 'Kasprza k et al] |2006| Upper panel: 
Sketch of a planar semiconductor microcavity delimited by 
two Bragg mirrors and embedding a quantum well (QW). The 
wavevector in the z direction perpendicular to the cavity plane 
is quantized, while the in-plane motion is free. The cavity 
photon mode is strongly coupled to the excitonic transitions in 
the QWs. A laser beam with incidence angle 6 and frequency 
uj can excite a microcavity mode with in-plane wavevector 
fc|l = '^sinO, while the near-field (far-field) secondary emis- 
sion from the cavity provides information on the real-space 
(k-space) density of excitations. Central panel: The energy 
dispersion of the polariton modes versus in-plane wavevector 
(angle). The exciton dispersion is negligible, due to the heavy 
mass of the exciton compared to that of the cavity photon. 
In the experiments, the system is incoherently excited by a 
laser beam tuned at a very high energy. Relaxation of the 
excess energy (via phonon emission, exciton-exciton scatter- 
ing, etc.) leads to a population of the cavity polariton states 
and, possibly, Bose-Einstein condensation into the lowest po- 
lariton state. Lower panel: Experimental observation of po- 
lariton Bose-Einstein condensation obtained by increasing the 
intensity of the incoherent off-resonant optical pump. 
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few more years to be achieved: after some preliminary 
claims ( Dang et al. [1998 Deng et al. 2002 Richard 



et al. 2005a|b ), a conclusive demonstration was reported 



by |Kasprzak et arj|200 6': the onset of Bose-Einstein con- 
densation in a gas of exciton-polaritons was assessed both 
in k-space from the macroscopic accumulation of parti- 
cles into the low-energy states and in real space from 
the appearance of long-range coherence. The principle of 
these experiments is illustrated in Figjlj the laser pump 
injects hot electron-hole pairs, whose excess energy can 
be dissipated via phonon emission and then Coulomb 
scattering processes. For sufficiently high pump pow- 
ers, the density of the incoherently injected polaritons 
exceeds the critical density for Bose-Einstein condensa- 
tion and a coherent condensate appears. Even if the sys- 
tem is still a driven-dissipative one, a quasi-thermal equi- 
librium state seems to be achieved once the loss rate is 
sufficiently slow as compared to the thermalization time 
of the polariton gas. For a more detailed discussion of 
the issues related to the relaxation mechanisms and the 
thermalization of the polariton gases, we refer the reader 
to the review papers that are already available on th e 
subject, e.g. ( |Deng et al\ |2010| [Keeling et~oI\ |2007[ ). 



While most of this research is being carried out in devices 
fabricated with inorganic GaAs- or CdTe-based alloys, a 
first observation of spontaneous polariton coherence in 
a non-resonantly pumped organic single-crystal micro- 
cavity was very recently reported in (Kena-Cohen and 



Forrest 20101 



The possibility of using exciton-polariton gases for 
studies of many-body physics and, more p recisely, super- 
fluid hydrodynamics was first proposed in ( Carusotto and 



Ciuti 


20041: 


et al. 2001 



s adopted. In contrast to the case of 
off-resonant and incoherent excitation, in the resonant 
driving configuration it is possible to perform an ah 
initio theoretical description in terms of a generalized 
Gross-Pitaevskii equation without having to cope with a 
phenomenological description of complex relaxation pro- 
cesses. In spite of the simplicity of the system, the coher- 
ently injected condensate shows peculiar superfluid prop- 
erties when it hits a defect in the cavity. The shape of 
the resulting density perturbation could be interpreted 
in terms of the celebrated Landau criterion of superflu- 
idity using the generalized Bogoliubov dispersion of ex- 
citations in the non-equilibrium condensate: at low flow 
speeds, superfluidity manifests itself in the suppression 
of the real-space modulation around the defect and, cor- 
respondingly, in the disappearance of the Rayleigh scat- 
tering ring in fc-space. For larger flow speeds peculiar 
patterns appear in both the density and momentum dis- 
tribution of the condensate polaritons. Experimental ver- 



ification of these predictions was reported in ( Amo et al. 
2009a) and is summarized in Fig[2] where the transition 



from a dissipative flow (panels I and IV) to a superfluid 
one (panels III and VI) is apparent in both the real- and 
the momentum-space images. Remarkably, in the same 
work it has been directly shown that photon-photon in- 



1 

4 


II 

n 


III 






VI 


-0.5 0.5 


O -0 5 □ 3.5 


-0,5 □ 0.5 1,0 


1 



JC LLrn 


t 




IV 

W , 


V 


VI 

• 


-0.5 D.S 


.0 -0.5 □ 0.5 


0-0 5 □ 0.5 1.0 
A, (tun"') 



FIG. 2 Figures from |Amo et al.\ |2009a| Polaritons are co- 
herently injected into the microcavity by a nearly resonant 
laser field: in contrast to the non-resonant and incoherent 
pumping scheme of Fig. [T] this pumping scheme allows to 
precisely control the density and the in-plane flow speed of 
the polariton fluid by changing the parameters (intensity, fre- 
quency, incidence angle) of the driving laser. Upper group 
of panels: experimental images of the real- (panels I-III) and 
momentum- (panels IV- VI) space polariton density extracted 
from the near-fleld and far-fleld emitted light from the cav- 
ity. The different columns to increasing values of the polari- 
ton density from left to right refer: for the highest density 
value, polariton superfluidity is apparent as a suppression of 
the real-space density modulation (panel III) and the corre- 
sponding disappearance of the Rayleigh scattering ring (panel 
VI). Lower group of panels: corresponding theoretical results 
obtained by numerically solving the non-equilibrium Gross- 
Pitaevskii equation of (Carusotto and Ciuti 20041. 



teractions are responsible for the appearance of a sound 
mode in the polariton fluid, as attested by the appearance 
of a Cherenkov-Mach cone at "supersonic" flow speeds. 
Following experimental works have extended the study 
to strong defects, shown the hydrodynamic nucleation 



of vortex-antivortex pairs ( Nardin et al. 


2011a i 


janvitto 


et al. 


2011 


and of dark solitons( Amo et al. 


2011 


Grosso 


et al. 


2011 


] in the flowing superfluid. A remarkable the- 



oretical development of the last few years is to use light 
superfluids in planar geometries to study quantum hydro- 
dynamics effects, in particular the analog Hawking radia- 



tion from acoustic black hole configurations ( 


Gerace and 


Carusotto 


2012 


Marino 2008 


Solnyshkov et al. 


20111. 



In all these works, interactions between single photons 
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FIG. 3 Top panel: figure from |Hartmann et aL\ |2006| In 
an array of coupled optical cavities, photon hopping occurs 
thanks to the spatial overlap of the photon modes of adja- 
cent cavities. The strong optical nonlinearity is induced by a 
coherently driven atomic gas present in the cavities. Middle 
panel: figure from |Greentree et a/. | [2006} showing a schematic 
diagram of a two-dimensional array of photonic crystal cavi- 
ties, each cavity containing a single two-level atom. Bottom 
panel: a Jaynes-Cummings-Hubbard system obtained with 
superconducting quantum circuits. Each cavity consists of 
a superconducting transmission line resonator embedding a 
superconducting qubit, that is an artificial two-level atom. 
Figure from |Koch ei al.\ [20T0l 



are weak and the liydrodynamic behavior of the photon 
gas originates from the collective interactions of a large 
number of coherent photons sharing the same orbital 
wavefunction. Almost in the same period, a series of pio- 



neering theoretical works ( 


Angelakis et al. 


2007 IGreen- 


tree et al. 


2006 


Hartmann et al. 2006 


I made the first 



steps in the theoretical investigation of a completely dif- 
ferent regime of strongly interacting photon gases, where 
interactions between single photons are large enough to 
induce sizable quantum correlations in the photon gas. 
The starting point of these proposals is the so-called pho- 



ton blockade phenomenon predicted in (Imamoglu et al. 



1997) and then experimentally observed in (Birnbaum 



et al. 2005): when the optical nonlinearity of a single- 
mode cavity is large enough for a single photon to shift 
the resonance frequency by an amount larger than the 
cavity linewidth, a resonant laser is able to inject pho- 
tons in the cavity only one at a time and the stream of 
transmitted photons are strongly antibunched in time. 



When this strong optical nonlinearity is inserted in a lat- 
tice geometry with many coupled cavities, one can expect 
that the photon gas will show the rich physics of the 



Bose-Hubbard model (Fisher et al. 1989), including the 



superfluid to Mott-insulator transition recently observed 
in atomic gases ( Greiner et al. 2002 ) . The most promis- 



ing systems to experimentally address this physics are 
illustrated in the different panels of FigjSj macroscopic 
cavities filled by an optically dressed atomic gas; an array 
of photonic crystal cavities embedding two-level emitters; 
superconducting circuits embedding artificial two-level 
atoms based on Josephson junctions. Soon after, these 
proposals were generalized to other many-body states, in 
particular Tonks-Girardeau gases of impenetrab le pho- 
tons in a one-dimensional hollow fiber geometry ( Chang 



et al. 2008 ) 



et al. 2008) and fractional quantum Hall states (Cho 



In contrast to the developments in superfluid hydro- 
dynamics of dilute photon fluids, none of these pioneer- 
ing works on strongly correlated gases speciflcally ad- 
dressed the consequences of the losses that are unavoid- 
ably present when dealing with light. Dissipation was in 
fact considered only as a hindrance limiting the available 
time for manipulation and observation of the quantum 
state, with a most detrimental effect on the Mott insu- 
lator state. A new perspective on strongly correlated 



photon gases was introduced in ( Gerace et al. 2009 ) for 
the case of a two site photonic Josephson system and, 
soon later, in ( Carusotto et al. 2009 ) for the case of a 



Tonks-Girardeau gas: a full inclusion of the interplay 
of driving, dissipation and strong interactions into the 
model offers the opportunity to observe novel dynamical 
features in such a driven-dissipative photon gas, and also 
suggests new tools for its experimental manipulation. In 
addition to this, the driven-dissipative regime allows for 
novel mechanisms of photon blockade to be exploited in 



coupled cavity systems ( Liew and Savona 2010), (Bamba 



and Ciuti 2011 Bamba et al. 2011). Very recently, the 



analog y with on-going devel opments in ultracold atomic 
gases (Dalibard et al. 2011) has opened a new research 



line in the direction of creating synthethic gauge fields for 
photons: with a suitable tailoring of the photonic envi- 
ronment, the motion of a photon can be made to expe- 



rience an effective vector potential ( 


Hafezi et al. 


2011b 


Koch et al. 


2010 Umucalilar and Carusotto 2011 


Wang 


et al. 


2009 


) and novel quantum states of the photon fluid 



analogous to the fractional quantum Hall effect can be re- 



alized ( Cho ei aL 2008 Umucalilar and Carusotto 2012). 



This review article is organized as follows. In Sec. |TTj 
we briefly summarize the main tools of elementary quan- 
tum field theory that are used to theoretically describe 
the non-equilibrium and quantum physics of the photon 
fluid in nonlinear optical devices. Even though a partic- 
ular attention is paid to planar microcavity geometries 
where most of the recent experimental observations were 
obtained, most of the concepts are fully general and can 
can be transferred to other systems. The generalized 
Gross-Pitaevskii equation describing the dynamics of the 
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photon fluid at mean-field level is introduced in Sec III 
particular emphasis will be devoted to the new features 
stemming from the non-equilibrium nature of the sys- 
tem. The most significant static and dynamic proper- 
ties of the coherent photon fluid under the different ex- 
citation schemes are reviewed in the following three sec- 
tions: Sec |IV| deals with the case of a coherent and quasi- 
resonant driving, Sec |V] deals with the optical parametric 
oscillation case when the cavity is pumped in the vicin- 
ity of the inflection point of the lower polariton branch, 
and Sec |VI| deals with the case of a incoherent pumping. 
We will review the consequences of the non-equilibrium 
condition on the condensate shape both in real and mo- 
mentum space as well as its impact on the elementary 
excitation spectrum. The superfluidity properties of the 



photon (polariton) fluid are reviewed in Sec VII among 
the several signatures of superfluidity, most emphasis will 
be devoted to the density modulation induced in the pho- 
ton fluid by a weak impurity, for which impressive experi- 
mental observations have been recently obtained. Subtle 
issues related to the generalization of the Landau cri- 
terion and the role of the interaction-induced speed of 
sound will be extensively discussed. More complex hy- 
drodynamic effects involving the nucleation of solitons 
and vortices in the wake of a large and strong defect are 
reviewed in Sec |VIIl| The review of the emerging fleld 
of strongly correlated photons is the subject of Sec. |IX| 
we will discuss in detail the different microscopic mech- 
anisms that can be used to achieve an efficient photon 
blockade regime and the rich new physics that has been 
predicted for lattices of strongly nonlinear cavities, in 
particular in the presence of an synthetic gauge field. Fi- 
nally, conclusions and future perspectives are drawn in 
SecH 



II. QUANTUM FIELD DESCRIPTION OF NONLINEAR 
PLANAR CAVITIES 



In this Section, we will review a second quantization 
formalism approach to describe the physics of 2D nonlin- 
ear cavities. Our goal is to present the essential quantum 
field theoretical tools that will be useful to understand 
the discussion of quantum fiuid effects in the following 
sections. While we will give particular emphasis to semi- 
conductor microcavity systems where most of the exper- 
imental observations have been carried out, much of the 
theoretical concepts reviewed here can be applied to ar- 
bitrary planar optical resonators embedding a nonlinear 
slab and are easily generalized to other geometries. We 
will be careful in pointing out when some properties are 
specific to a given system or not. Comprehensive intro- 
ductions to semiconductor microcavity systems can be 



found in the specific literature, for instance ( 


Deveaud 


2007 


Deveaud et al. 


2003 


2009 


Kavokin et al. 


20081. 



A. Free cavity fields and input-output formalism 

1. The two-dimensional photon field 

Since a planar cavity is by definition invariant under 
in-plane translations, the in-plane wavevector k is a good 
quantum number for the free photon dynamics, which 
can be described by an Hamiltonian of the form 



^cav 



/iwcav(k)aj. ,,(k)ac,^(k) (4) 



where the ac^o-(k) and aj^^(k) operators respectively 
destroy and create a cavity photon with an in-plane 
wavevector k and a polarization state a. The creation 
and destruction operators satisfy standard Bose commu- 
tation rules 

[aa.(k),a[,,,,(k')] = (2^)^ 5(k - k') <5,,.' (5) 

[aa.(k),ac,a'(k')] = 0. (6) 

Two-dimensional real-space cavity photon field operators 
^c,cr(i") and ^{y) are defined as the Fourier transform 
of ac,cr(k). 



dk 

(2^ 



(7) 



and again satisfy Bose commutation rules 

[*c,.(r),*^,.,(r')] = <5(r-r')<5.,.', (8) 

[*c,.(r),*c,.'(r')] = 0. (9) 

While the nco-Cr) = ai^) ^c,cr(i") operator can be in- 
terpreted as a sort of two-dimensional photon density on 
the cavity plane, the physically observable electric field 
at the three-dimensional position (r, z) is expressed in 
terms of the field operators by 



E(r,^) 



d^k 



(2^) 



2 



e*'^'-E,(k, z)ac,, (k)+h.c. (10) 



where the z-dependence of the photon mode wavefunc- 
tion ECT(k, z) keeps track of the metallic or dielectric na- 
ture of mirrors. In the simplest case of a planar cavity 
of refractive index no enclosed between metallic mirrors 
spaced by a distance ^z^ the z-dependent profile ECT(k, z) 
of the lowest mode has a sinusoidal shape extending for 
< z < 4, 



E^(k,z) = i. 



■ e„ sm 



TT z 



(11) 



the boundary conditions at the mirrors fix the position 
of the nodes, and the global amplitude of the field is de- 
termined by the energy of a single photon state hujcav (k) . 
The unit vectors are a basis of light polarizations. A 
detailed discussion of the more complex case of dielectric 
cavities enclosed between a pair of Distributed Bragg Re- 



flector (DBR) dielectric mirrors can be found in ( Savona 
T999|. 
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2. Polarization effects and effective spin-orbit interaction 
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FIG. 4 Top panel: snapshots of the propagation of cavity 
excitations created by a tightly focused and linearly polar- 
ized laser pulse in a semiconductor planar microcavity. Left 
and right panels show the optical intensity in the same (left) 
and in the orthogonal (right) polarization. The polarization- 
dependent patterns a re d ue to the TE/TM splitt ing of the 
planar cavity modes (1121). Top panel taken from Langbein 



et al. 



2005 Middle panels labeled as (a-c): the effective opti- 
12|) is equivalent to a k-dependent 
that rotates the pseudo spin 



cal spin-orbit interaction 
magnetic field (orange arrows) 
associated to the photon polarization state on the Poincare 



2005a I 



sphere. Lower panel: as anticipated in ( Kavokin e£ a/, 
the combination of the effective spin-orbit interaction with 
disorder-induced resonant Rayleigh scattering of a linearly po- 
larized pump leads to the optical analogue of the spin Hall ef- 
fect that is visible in the k-space pattern of the z-polarization 
shown in the bottom panel. Middle and lower panels taken 
from |Leyder et fflZ.|[2057l 



In writing the Hamiltonian we have implicitly as- 
sumed that the polarization states are degenerate. While 
this assumption is exact at fc = 0, the reflection ampli- 
tudes off a dielectric mirror for the TE (Transverse Elec- 



tric) and TM (Transverse Magnetic) linear polarization 
states are generally different at k 7^ 0, introducing a fre- 
quency splitting of t he TE/TM modes prop ortional to fc^ 
in the small k limit (Panzarini et al. |1999 ). 

Physically, the resulting TE/TM splitting can be inter- 
preted as a kind of spin-orbit interaction term coupling 
the orbital (the wavevector k) and pseudospin (the polar- 
ization) degrees of freedom of the cavity photon and can 



be described by a Hamiltonian term of the form ( Kavokin 



et al. 2004 Shelykh et al. 2010 1 



J2 ^so{k) \{kl~kl)ai^J,+2k^kya 



iv) 

aa' 



d^k 



X a[, ,^(k,i)ac,a/(k,t) = 
^ fc2 j^^^(^) rcos(2(^) a^^l, + sin(2^) a^^j, 



X a^^^(k,Oac,.'(k,i) (12) 

where a'^^''^^ are the Pauli matrices, the polarization in- 
dices cr, a' run here over the circular polarizations basis 
a± and f2so(fc) quantifies the fc-dependent magnitude of 
the TE/TM splitting. In the alternative form shown in 
the last line, is the angle between the wavevector k and 
the X axis on the plane ^ 

If one represents the two a± polarization states of the 
cavity photon as a 1/2-pseudospin, then the spin-orbit 
interaction (12) is equivalent to an effective momentum- 



dependent magnetic field, which induces a precession of 
the pseudospins. As shown in Fig. [4j this kind of op- 
tical spin-orbit interactions can give rise to spectacu- 
lar effects such as a strongly anisotropic polarization- 
dependent propagation (iLangbein et al\ 120051). When 



combined with disorder-induced scattering, this effect is 
responsible for an optical analogue of the spin Hall ef- 
fect pCavoHirTFan |2005a| iLe^^ More 



recently, this same p 



photonic cavity in ( |Maragkou et al. 



lysics was investigated in a purely 
2011) . 



It is worth noting that the term spin Hall effect of light 
is sometimes used in the photonic literature to denote a 
different family of effects stemming from the spin-orbit 
coupling experienced by photons propagating in bulk op- 
tical media with weak spatial inhomogeneities of the re- 
fractive index (|BUokh and Bliokh[ |2004a|b[ |2006| [Liber 



man an d Zcl'd ovich[|1992[|Onoda et a/.H2004f ^ecently, 



this coupling was shown to have remarkable consequences 
such as a sizable lateral shift of the trajectory of a light 



beam beyond the geometrical optics prediction (Bliokh 
et al. 



2008 Hosten and Kwiat 



2008| ). 



^ At a naive glance it may appear that the spin-dependent term 
\12\ breaks rotational symmetry: the point is that the Pauli 
matrices act here on the space of <t± polarized states of a spin-1 
(and not 1/2) object such as the photon, which explains their 
angular momentum being equal to 2 and t hen the factor 2 in the 
trigonometric functions on the last line of (|12|l. 



10 



3. Pumping and losses: input-output theory and master 
equation 

A quantum description of the driving of the cavity pho- 
ton mode by an incident coherent laser beam can be ob- 
tained using the input-output theory of optical cavities. 
A complete discussion of this theory for single mode cavi- 
ties can be found in quantum optics textbooks (Gardiner 



andZoller 2004 Walls and Milburn 2006). Its extension 



to the case of planar microcavities with a continuum of 



in-pla ne k modes can be found in ( Ciuti and Carusotto 
[2006) : the in-plane translational symmetry of the device 
guarantees that the in-plane k vector is conserved, so 
that an external radiation of frequency uj and incident 
angle ^inc only couples to the cavity field component of 
in-plane wavevector k = lo/c sin^inc- Correspondingly, 
the cavity field component of in-plane k can only decay 
into external radiation emitted at an angle ^out satisfying 
the analogous condition k = ui/c sin 0out : this latter con- 
dition is schematically illustrated in the right-most part 
of the upper panel of Fig[T] 

The Hamiltonian term describing the external driving 
of the cavity by a coherent incident field of amplitude 
E™'^{Y,t) can be written in k-space as 



-,7r(k)^r*(k,i)aa.(k)V (13) 



where £^J^"'^(k, t) is the Fourier transform of E^J^^{r,t) 
and the r/^J (k) coefficient is proportional to the transmis- 
sion amplitude of the front mirror for light with in-plane 
wavevector k. 

The finite transmittivity of the front and the back mir- 
rors of the cavity is responsible for the re-emission of light 
from the cavity with an amplitude proportional to the in- 
cavity field operator. Combining this secondary emission 
with the direct reflection of the coherent laser light off 
the front mirror, input-output theory leads to the form 



&^\k,t) 



.,back 



ik)acA^,t) (14) 
{k,t) + K'^(k)acA^^t) (15) 



for the quantum operators describing the transmitted 
and reflected fields, respectively. Analogously to the 
input t^ctC^) coefRcients in (13 1, the output coefRcients 



,fr,back 



(k) are proportional to the transmission ampli- 
tude of the mirrors and can depend on the in-plane 
wavevector k and the polarization a. An implicit as- 
sumption of this formalism is that transmission through 
the mirrors is almost instantaneous, so that the relations 



between the external and in-cavity fields ( 13 ) are local 
in time. In frequency space, this means that the 77 and 
K coefficients do not depend on the frequency uj. The 
validity of this Markovian assumption is intuitively un- 
derstood for thin metallic mirrors; a discussion for di- 



electric DBR mirrors was recently reported (Sarchi and 



Carusotto 2010) 



As usual, the emission of light by the cavity is ac- 
companied by a radiative damping of the cavity field at 
a wavevector- and polarization-dependent rate 7^^'^(k), 
proportional to the sum of the mirror transmittivities 
|K^'(k)p -I- |/tj^^'^'*(k)p. Tracing out the radiative modes 
of the field outside the cavity, dissipation results in addi- 
tional terms to be included in the master equation for the 



evolution of the density matrix of the cavity field (jBreuer 



and Milburn 2006), 



and Petruccione 2002 Gardiner and Zoller 2004 Walls 



dp 
'di 



ih 



m,p] + c[p] 



(16) 



Under the assumptions that the temperature of the ra- 
diative modes outside the cavity is much lower than the 
frequency of the cavity mode fc^T <C ^cav and that the 



radiative coupling is Markovian ( Breuer and Petruccione 



2002), the super-operator C[p] accounting for the dissi- 
pative effects has the zero temperature Lindblad form 



rad r 



2ac,aik)pal,^(k)- 



(2^ 2 

acAk)dc,a{k)p- pal.Ak)dc,aik) . (17) 



with a frequency-independent decay rate 7^'*'^(k). Addi- 
tional, non-radiative decay channels due, for instance, to 
absorption in the cavity material can be included in the 



model by including into the master eq uati on ( 16 ) addi- 
'^'^[p] of the same form ( [l7| and propor- 
^^). In the fol- 

rad(k)+7^''^'l(k) 



tional terms [p] of the same form 
tional to the non-radiative loss rate 7™ 
lowing we shall indicate with 7^ (k) = 7, 
the total decay rate of a cavity photon. 

A complete calculation of the value of the 77, k and 
7 parameters for specific configurations goes beyond the 
scope of this review and generally requires a microscopic 
solution of Maxwell equations for the field propagation 
across the device; in a planar geometry, a tool commonly 
used for this kind of calculations is the so-called transfer 
matrix method, reviewed e.g. in (Burstein and Weis- 



buch 1995 Savona 19991. Useful general relations can 



be mentioned for the case when the front and back mir- 
rors of the cavity are identical: under this assumption, 
r]^^''^ (k) , (k) = t4'"=^{k) and the radiative de- 
cay rate Ya^'ik) = 2 77,(k) K<,(k) 



e(k) 

Furthermore, explicit 
expressions for the ij^ik) and Ka{k) coefficients can be 
given for the case of a planar cavity enclosed by loss- 
less metallic mirrors of transmittivity t <^ 1 separated 
by a distance i^, namely 77 = [ctjll^^^ l-^jixhujc and 

K — ~t y^Trhujc/^z- 



B. Optical nonlinearities and effective pfioton-pfioton 
interactions 

So far, we have discussed the dynamics of the cavity 
field at the level of linear optics, where the quantum dy- 
namics of the non-interacting cavity field reduces to the 
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classical wave equation stemming from Maxwell's elec- 
trodynamics in material media. The situation changes 
when the cavity layer (or the mirrors) embeds a material 
with a sizable optical nonlinearity: the non-linear depen- 
dence of the matter polarization on the applied electric 
field is responsible for a number of wave-mixing processes 
coupling different cavity modes and generating new fre- 
quency components (Boyd) |2008 Butcher and Cotter 
2008). Moreover, it can lead to strong modifications of 



the quantum fluctuation properties of the cavity field, 
such as a reduced noise on some field quadrat ure or even 
the generation of entangled states for the field (Walls and 



Milburn 2006[ ). The present review being devoted to the 
quantum fluid aspects of the photon gas, we will con- 
centrate our attention on third-order nonlinearities pro- 
portional to the x^^^ nonlinear polarizability that can be 
described in terms of binary interactions between pairs 
of photons. 

Under the standard rotating-wave approximation (that 
is valid here provided the photon mass is larger than all 
other energy scales, e.g. losses, kinetic energy, interac- 
tions), the total number of photons is conserved and the 
nonlinear process can be described by a four-operator 
Hamiltonian term of the form, 



U 



cav— cav 



1 r d'^k 

2 J (2^ J (27r)2 J (27r)2 



cav— cav 



(k,k',q)x 



xaj._^(k + q)a^^^,(k'-q)ac,.-(k')aa.(k) : (18) 

where the matter degrees of freedom have been traced out 
and summarized into the effective photon-photon interac- 
tion potential V^''Jr''''^(k, k', q). The polarization index 
runs over the circular polarization states a — a± and 
spin angular momentum is implicitly assumed to be con- 
served in the collision process. Total momentum is also 
conserved in the process of two photons of initial in-plane 
wavevector k and k' scattering into the new wavevector 
states of in-plane momentum k -I- q and k' — q. 

As the typical length scale of the electron dynam- 
ics in typical bulk material media and in semiconduc- 
tor hctcrostructures (Bastard 1991 ) used for quantum 



fluid effects is much shorter than the optical wavelength 
along the plane, the interaction potential v^'ij,^'^^^ can 
be approximated with its zero-momentum value V°^, for 
k = k' = q = 0: in real-space, this corresponds to as- 
suming that photon-photon interaction occur via a local 
potential, 



U 



cav— cav 



x^c..a'ir)'^cA^). (19) 

In two and three dimensions, the use of a strictly lo- 
cal interaction potential beyond the Born approxima- 
tion often leads to UV divergences: standard techniques 
to make the theory regular involve renormalization of 



the interaction potential on a discrete lattice ( Mora and 
Castin [2003 ) or the use of suitably defined pseudo- 
potentials (IHuangl 119871). 



In the simplest case when the photon frequencies that 
are involved in the photon fluid dynamics are very far 
away from electronic resonances in the nonlinear opti- 
cal material, the optical transitions are virtual and the 
population of the excited electronic states remains neg- 
ligible. In this regime, the photon-photon potential can 
be expressed as V°^, — — ""^(^"cav) j,^ terms of the 
nonlinear refractive index n„/ normalized in a way such 
that the effective refractive index is ricg — nQ + n^i |f p. 
A more sophisticated theory of photon-photon interac- 
tions in a planar cavity device embedding a collection of 
anharmonic atoms appeared in (Chiao et al. 2004). In 
the case of dielectric cavities, the cavity layer thickness 
has to be replaced by the effective thickness of the cavity 
mode £ cav ( Savona 1999). 



While sitting far from resonance generally allows to 
minimize undesired absorption losses in the material, 
non-resonant optical nonlinearities have the serious draw- 
back of being generally very small and requiring a large 
number of photons to be observable. While this is not too 
much a concern for quantum fluid experiments, strong 
nonlinearities are required to observe photon blockade 
effects and generate strongly correlated photon fluids. In 
the next Subsection we shall extend our discussion to 
systems where the cavity mode is strongly coupled to 
narrow transitions in the optical medium: in this case, 
the photon inherits the strong nonlinearity of the matter 
excitations. 



C. Strong light-matter interaction and cavity polaritons 

In the last two decades, the physics of strong light- 
matter interaction has flourished in many interest- 
ing domains, including the fields of atomic cavity 



QED (Bcrman et al., 


1993 


Fleischhauer et al. 


2005 


Rai- 


mond et al. 2001 


), semiconductor microcavities and su- 



perconducting circuit QED ( Schoelkopf and Girvin 2008 



You and Nori 201 1 1 . In this section of our review we shall 



concentrate our attention on the case of planar semicon- 
ductor microcavities embedding one or many quantum 
wells (Deveaud[ 2007 Deveaud et al. 



2003 2009 Ka 



vokin et al. 2008 ). A brief account of other systems will 



be given later in Sec IX 



1. Quantum well excitons coupled to the cavity mode 



As detailed in more specialized reviews such as |Bas-| 
|tard[ [T991| and [Deveaud et al.\ |2003[ a quantum well 
(QW) consists of a thin semiconductor layer (a few nm 
thick) embedded in a different semiconductor compound 
acting as 'barrier' material. The chemical composition of 
the well is chosen to have the bottom of the conduction 
(the top of the valence) band at a lower (higher) energy 
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than the surrounding material, thus producing quantum 
confinement of both electrons and holes. The lowest en- 
ergy optical transition corresponds to the excitation a 
two-dimensional (hydrogen- like) electron-hole pair con- 
fined in the QW layer, the so-called exciton. In typical 
microcavity samples, one or more QWs are embedded in 
the cavity layer, with their plane parallel to the cavity 
plane. To have a strong and quasi-resonant coupling of 
their electronic degrees of freedom with the cavity mode, 
the QWs are placed at the antinodes of the cavity field 
and the cavity mode frequency u!°^^ is tuned in the vicin- 
ity of the lowest QW exciton with frequency Wexc- 



to the cavity field can be written as 



At the level of linear optics ( Bastard 19911 Klingshirn 



2007), the contribution of the quantum well exciton to 



the optical properties of the cavity can be described in 
terms of a resonant contribution to the dielectric polar- 
izability of the structure of the form 



Siz- ZQw) /; 



2D 



,(k) - - i7x/2' 



(20) 



where the quantum well is approximated as a very thin 
layer located at zqw and the excitonic transition has 
an oscillator strength surface density f2D- The exciton 
resonant frequency is weakly dependent on the in-plane 
wavevector k, W(,xc(k) — w"^^, -I- hk^ /2mcxc with an exci- 
ton mass TOoxc of the order of the electron mass, i.e. or- 
ders of magnitude larger than the effective photon mass 
Wcav The linewidth "fx accounts for all non-radiative 
decay channels of the exciton. 

In a quantum picture, the quantum well exciton can 
be described in terms of destruction (creation) opera- 
tors ax,aO^) (axcr(k)) that destroy (create) an exciton 
with total momentum k and their real-space counterparts 
^x,a and i>l( ^ defined as their Fourier transform as in 



([t]). In standard QWs ( Bastardj [1991 ), electrons can 
have spin projection tJe = ±1/2 along the growth axis 
z, while holes can have ah = ±3/2. Among the four 
exciton states with spin projection a = ±1,±2 that ex- 
ist in a QW, here we shall restrict our attention to the 
a = ±1 that are coupled to the cavity mode; the ones 
with cr ± 2 are optically inactive. Provided the inter- 
particle distance remains much larger than their Bohr 
radius, excitons behave as bosonic particles, whose cre- 
ation and destruction operators satisfy Bose commuta- 
tion rules of the same form as equations (5][6) for cavity 
photons. The corrections to bosonic behavior due to the 
composite nature of the exciton have been theoretically 
addressed in ( Combescot et al. |2008| ) . Discussions of ex- 
citon physics with a special eye to planar microcavities 
can be found in the recent review papers by Deveaud] 
eTarj [20031 [Deveaud et [2009| and |Deng et aL[|2010| 



In a second quantized formalism, the Hamiltonian de- 
scribing in k space the exciton dynamics and its coupling 



|^'^cxc(k)a^,^(k)ax,a(k) + 
'a^xJk)acA^} + al jk)axA^)] } (21) 



The last term describes the coherent conversion of an 
exciton into a cavity photon at a Rabi frequency flu pro- 
portional to the product of the electric dipole of the quan- 
tum well exciton transition times the local amplitude of 
the cavity photon electric field. In terms of the classical 



dielectric properties ( 20 ) of the quantum well, the Rabi 



frequency flu in a planar metallic cavity can be related 
to the exciton oscillator strength f2D by 



47ra;" 



2D 



1/2 



Ejz) 
E 



(22) 



where the final fraction accounts for the displacement of 
the quantum well position from an antinode of the cavity 
mode. In the general case of DBR cavities, the distance 
iz between the mirrors has to be repl aced by the ef fec- 
tive thickness £cav of the cavity mode (Savona 1999). If 
Nqw quantum wells are present in the cavity, Nqw — 1 
linear combinations of the exciton states are dark, while 
the single bright one is coupled to the cavity mode with 

an enhanced couphng Q.fi = \/X)ilT^ I^i?^]^- When the 
QWs are located in equivalent positions, the collective 

1/2 

coupling enhancement factor is Nq^. Note that, in con- 
trast to the case of a quantum well in free space (An- 



dreani et al. 1991 Tassone et al. 1992), the presence 



of the microcavity eliminates the direct coupling of the 
quantum well exciton to the external radiative modes and 
radiative decay of the exciton can only take place via the 
lossy cavity mode. 

In addition to these effects stemming from its cou- 
pling to radiation, the exciton is also subject to non- 
radiative recombination processes. In a master equation 
formalism, these processes can be described by a Lind- 
blad term of the form ( 17) with a (generally weak) decay 
rate ^ x.a{k), as well as deph asing processes at a rate 
7xct^ (Liew and Savona 2011) due to, e.g., interactions 
with carriers and spatial inhornogeneities in the quantum 
well thickness. The corresponding term in the master 
equation reads: 



dcph /I \ 
deph _ ^X,a \'^) 



c 



2a^ (k)a^,^(k)pa^ (k)ax,^(k)- 



- aT^^^(k)ax,a(k)a5f ,^(k)ax,^(k)p + 

- pa^^,(k)ax,.(k)a^^^(k)ax,.(k)] (23) 

In view of the on-going theoretical and experimental de- 
velopments, it is important to remind that the form \21\ 
of the light-matter coupling term is based on the so-called 
rotating-wave approximation where anti-resonant terms 
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proportional to ax,aO^)a,c.a{)<) are neglected. This ap- 
proximation is very accurate as long as flu <C WcavjWexc, 
which is generally the case for microcavity systems under 
examination here where is of the order of 10 meV for 
typical III-V based samples and of the order of 20 meV 
for II- VI based samples. Significant deviations appear 
in the opposite ultra-strong coupling regime where flji is 
comparable or larger than ojcavjWcxc- Such a regime has 
been recently theoretically predicted and experimentally 
observed in systems where the cavity mode is strongly 
coupled to intersubband electronic transitions in doped 



quantum wells ( 


Ciuti et al. 


2005 


), ( 


Anappara et al. 


2009 


Guenter et al. 


2009 




Todorov et al. 


20101, cy- 



clotron transitions of two-dimensional electron gases ( Ha- 



genmiiller and De Liberato 


2010 Scalar! etal. 


2012 1 and 


in circuit-QED systems ( Devoret et al. |2007 1 , 


( Bourassa 


etal. 2009 Nataf and Ciuti 


20111, (Fedorov et a?. 2010 


Niemczyk et al. 20101. In addition to going beyond the 



rotating-wave approximation in the light-matter coupling 



term in (21), a consistent theoretical description of sys- 
tems in the ultra-strong coupling regime also requires a 
more sophisticated modeling of the frequency dependence 



of the dissipative baths ( 


Carusotto et al. 


2012 


Ciuti and 


Carusotto 


2006 


Ridolfo et al. 


20121. 



2. Exciton-polaritons 

The coupled dynamics of photons and excitons in a 
microcavity is described by the Hamiltonian terms Q 
and (21) that involve the product of two field operators: 



as a result, the quantum dynamics exactly recovers the 
classical Maxwell wave equations inserting a thin layer 
of material with the susceptibility (20). Furthermore, 
simple linear combinations of the ac,a-0^) and ax,aO^) 



operators (Hopfield 



1958 ) can be used to transform the 



Hamiltonian into a diagonal form 



n cav cxc 
-72 



fuvLP,aik) a{p^^{k) aLP,„(k)-^ 
+ fiw;7P,dr(k)a|^p^(k)aj/p,o-(k), (24) 



crossing form 



^(UP,LP),a{'^) 



± 



2 ( ■^cav,cr(k) Wcxc,CT(k) 



1/2 



(25) 



that is shown in the middle panel of Fig{l] the bare 
cavity photon and exciton branches (dashed lines) are 
mixed by the light-matter coupling term into the polari- 
ton branches; the minimum splitting 2D,ji of the two po- 
lariton branches is obtained for a;cav(k) = Wcxc(k). The 
so-called strong light-matter coupling regime is defined 
when rifj > 7, i.e., light-matter coupling exceeding losses. 

The photonic and excitonic content of the polariton 
modes is quantified by the real-valued Hopfield coeffi- 
cients, which relate the cavity mode and exciton field 
amplitude to the ones in the lower and upper polariton 
modes, 



aa<T(k) 



.LP I 



.LP I 



(k) olpA^) + < (k) clupA^) (26) 
(k) aLP,.(k) + 7/^^(k) ac;p,,(k).(27) 



Close to the crossing point, both polariton modes 
have approximately equal photon and exciton content 
~ 1/2, while farther away the two polariton ac- 
quire a purely excitonic or photonic character. Differ- 
ently from the general case of (Hopfield 1958), the ab- 
sence of the a|^p jjp-^ ^ creation operators in the Bo- 
goliubov transormation ( 26|p7 ) is a consequence of the 
rotating-wave approximation. 

The effective masses mi^p up of the lower and upper 
polariton branches are calculated from the curvature of 
the dispersion (25). In the most remarkable case where 
uj°^^ « 10°^^, the condition mcav ^ "t-oxc implies that 
^LP,UP - ^?av,oxc T ^B. and niLp ~ m^/p ~ 2mcav; 
furthermore, the two polaritons have equal photon and 
exciton content as long as hk"^ /2mcav <C ilp. Most of 
the experiments that will be discussed in the next sec- 
tions only involve the lowest states of the lower polariton 
branch for small wavevectors and a description in terms 
of a single polariton branch with parabolic dispersion of 
mass ruLP is accurate. More details on the basic proper- 



ties of microcavity polaritons can be found in ( 


Deveaud 


2007 


Deveaud et al. 


2003 


2009 


Kavokin et al. 


2008 


) 



where the ajjp,cr{k) and dLp.o-(k) operators correspond to 
hybridized excitations resulting from the linear superpo- 
sition of exciton and cavity photon modes, the so-called 
exciton-polaritons. Thanks to the translational symme- 
try along the cavity plane, the k wavevector and the 
spin a remain good quantum numbers for the new quasi- 
particles. 

The dispersion of the upper and lower polariton 
branches uijp^ai)^) and ujLP.a^}^) has the typical anti- 



3. Polariton-polariton interactions 

As it is typically done in the many-body physics, a 
widespread procedure is to describe the system in terms 
of an effective model Hamiltonian that is able to repro- 
duce the exciton-exciton interactions without invoking its 
elementary constituents: the idea is to replace the com- 
plex Coulomb interactions between the electron and holes 
by a simple two-body interaction potential involving the 
exciton as a whole. In the simplest version, low-energy 
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scattering at a relative wave vector much smaller than 
the exciton Bohr radius can be accurately described by 
a contact two-body interaction potential term, 

era' 

x*;f,,,(r)#^^^(r), (28) 

where the spin indices cr, a' run over the circular polar- 
ization basis (7, a' ~ a±. Rotational invariance for a 
contact interaction potential imposes that total exciton 
spin is conserved and that V^^a-+ — 



(7 - ,(7- 



V^^ and 



V. 



XX 



XX 



- i^s ■ 



A first estimation of V'^'^ in terms of a microscopic 
quantum well electron-hole model of the exciton was pro- 



vided in (Ciuti et at. 1998 Tassone and Yamamoto 
1999 ) within the Born approximation, giving 



Ge^as/e where e is the dielectric constant of the material 
and as is the exciton Bohr radius. This positive and fea- 
tureless value of the triplet interaction constant appears 
to be in reasonable agreement with experimental mea- 



surements (Amo et al. 2009a Ferrier et al. 2011) and 



provides the most relevant contribution when the micro- 
cavity is excited by circularly polarized light. 

For other polarization configurations, the much richer 
physics of the singlet channel needs to be considered: dif- 



ferent experiments (Amo et al. 2010b Kavokin et al. 



2005b) have reported values of the singlet interaction 



constant with opposite signs; first systematic studies of 
the singlet interaction constant as a function of the micro- 



cavity parameters were recently performed in (Paraiso 
20101 Vladimirova et al. 2010). In analogy to the Fes- 



hbach resonance phenomenon in atom-atom scattering 
when the energy of the incident atomic pair is in the 



vicinity of a molecular intermediate state ( Chin et al. 
201 Oh, several authors (Carusotto et al. 2010[ Savasta 



et al.\ |2003| |1999| |Wouters[ |2007D have suggested the 
importance of the biexciton intermediate state (i.e. a 
two-electron, two-hole bound complex in the singlet scat- 
tering channel: in the vicinity of the scattering reso- 
nance, the effective interaction constant Vg''^ in the sin- 
glet channel is expected to be strongly enhanced and to 
change sign. More details on the theory of biexciton Fes- 



hbach effects will be presented in Sec IX. B. 3 

In principle, exciton-exciton interactions may also 
transform a pair of a± bright excitons into a pair of spin 
±2 dark excitons (Ciuti et al. 1998). While these pro- 



cesses are important in isolated quantum wells, they no 
longer conserve energy in microcavities: because of the 
Rabi coupling between photons and bright excitons, the 
dark spin ±2 excitons are at much higher energy than 
the lower polariton branch and can only play a role as 
non-resonant intermediate states of high-order processes. 

An additional interaction channel originates from sat- 
uration of the exciton oscillator strength( Glazov et al. 



2009 Rochat et al. 2000 Tassone and Yamamotol 1999 



Pauli exclusion principle for electrons and holes forbids 



that another exciton be created at a distance shorter than 
the Bohr radius from an existing exciton. At the lowest 
order in the exciton density and retaining only the bright 
exciton states, this can be modeled as an effective quartic 
Hamiltonian term of the form 



^sat — 



1 



(r)*^^^,(r)x 



X *;f,,,(r)*c,a(r)+li.c. (29) 

with a saturation potential V|°J, = S„ ,y/ Mlf^/risat. As 
it is discussed in detail in the recent paper by IGlazov| 
et al.\ [20091 the exact value of the saturation density 



Usat cx T/os depends on the specific shape of the wave- 
function of the relative motion. For typical values of 
semiconductor microcavities, the saturation term is gen- 
erally significantly smaller than the exciton-exciton con- 
tribution. 

As previously mentioned, in many relevant experimen- 
tal conditions, one can restrict to the bottom of the 
lower polariton branch and approximate the dispersion as 
parabolic. In this regime, we can rewrite the interaction 



Hamiltonians ( 28 1 and ( 29 ) in terms of polariton opera- 
tors and keep only the terms involving the lower polari- 
ton branch. This finally leads to an effective polariton- 
polariton contact interaction 

-Hlp-lp = \ J rf'r ^ Vi-;, *lp,.(r)*lp^,,(r)x 

x^'LP,a'{r)^lpJr) (30) 



with 
V, 



LP 
a .a' 



^l^' + 2 u^/ u^/ V^'. (31) 



Recent experiments ( Cristofolini et al. 2012) have 



demonstrated strong coupling of a cavity photon with 
a hybrid exciton in a double quantum well geometry, 
where the direct and indirect exciton states are mixed 
by coherent electron tunneling events across the barrier. 
As a result, the polariton acquires a finite electric dipole 
moment, that is expected to strongly reinforce polariton- 
polariton interactions and, possibly, induce new effects 
due to long-range dipole form of the interaction poten- 
tial ( Astrakharchik et al. 2007 Boning et al. 2011). 



D. External potentials affecting the in-plane motion of 
cavity excitations 

In analogy to magnetic and optical traps for atomic 
gases (IPitaevskii and Stringaril 120041), several strategies 



have been explored to implement spatially- and spin- 
dependent external potentials for cavity photons and ex- 
citons, creating Hamiltonian terms like 

^Pot- Jd'rJ2 [<7^(r)*L,.(r)^cav,.'(r)+ 

a,cr' 

+ V;(7'(r)*t (r)#exc..'(r)l. (32) 
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Provided the amplitudes of the potentials and 

^(exc) 

are much smaller than the Rabi energy Mln and 
for a sufficiently smooth spatial variation, the resulting 
potential acting on lower polaritons can be written as 



overgrowth technique able to create a position-dependent 
thickness of the cavity spacer (see top panel of Fig. [5|. 
This technique provides a quite flexible tool to design ex- 
ternal potentials with arbitrary shapes on characteristic 
spatial lengths in the few fini range. 

A more brutal way of confining polaritons is to de- 
sign a micropillar structure by etching away all the layers 
forming the top mirror and the cavity layer (and possi- 
bly also the lower mirror down to the substrate): in this 
case, light is confined in the in-plane directions by the 
large refractive index mismatch at the air-semiconductor 
sation experiment of lBalih et am2007| was created by interface (IBurstein and Weisbuchl [19951). As with the 



ej^(r) = \un' V^->{r) + \u--\-^ V;^7^(r). (33) 

We will start by considering several experimental proto- 
cols to generate scalar potentials of the form y^'^^''^^'^^ — 

y(cav,oxc) 

The polariton trap used in the Bose-Einstein conden- 



applying a mechanical stress, resulting in an energy red- 
shift of the exciton state ( Negoita et al. 1999 ) and thus 
an attractive potential for excitons {V^'^^'^' < 0) with an 
approximately harmonic profile at the bottom. 



(b) 




Bragg 
mirror 



GaAs 
AlAs 

I quantum well 
/.=835 nm 



Pulsed 





Potential landscape 



FIG. 5 Top panel: figure from |E1 DaiT et aL\ [20061 An ef- 
fective lateral confinement potential for cavity photons (and 
hence for cavity polaritons) can be obtained with a position- 
dependent cavity thickness (achievable in semiconductor mi- 
crocavities via growth, etching and regrowth). Bottom panel: 
figure from |Sanvitto et al.\ |2011| An effective lateral confine- 
ment potential for cavity photons can be obtained by using 
control laser spots with different polarizations as first achieved 



in (Amo et al. 2010b I. 



Another concept of polariton trap exploits the depen- 
dence of the cavity mode frequency on the thickness of 
the cavity layer. This feature is usually exploited to ob- 
tain a full scan of the exciton-photon detuning across 
the resonance on a single microcavity sample by grow- 
ing a wedge-shaped wafer (Weisbuch et al. 1992). In 



the experiment of (Sermage et al. 2001 ), a polariton ac 



celeration effect was observed as a consequence of the 
cavity wedge: within a parabolic band approximation, 
the polariton motion can be described by a Newton-like 
equation rtiLpd^r/dt^ = —VrfujLp{i^ = 0). 

A more sophistica ted development of this same idea 
was implemented in (|E1 Dai'f et al.\ |2006| |Kaitouni et al. 



2006 Lu et al. 2005 



regrowth technique, arbitrary geometries can be realized 



with full three-dimensional confinement ( Gerard et al. 



or with free motion along photonic wires (Kuther et al. 



1996[ [Ohnesorge et al.\ |1997[ [Reithmaier et al.\ |1997l 



1998[|Wertz et aL[ |2010[ [Zhang et al||1995[ ). The direct 



processing of micropillar devices may come at the price 
of larger lateral losses at the cavity/air interface. The 
recent realization of double pillar configurations show- 
ing a sizable tunnel coupling of the photonic wavefunc- 



tions (de Vasconcellos et al. 2011) is very promising in 



view of creating more complex structures such as cavity 
arrays. 

Another strategy to apply external potentials to 
exciton-polaritons is to deposit metal films on the sur- 
face of the cavity: depending on the geometry, this can 
either increase the photon energy by squeezing the tails 



of electric field inside the DBR mirror ( Kim et al. 2008 ) 



or even strongly mix the Tamm plasmon with the exciton 
and the cavity photon modes as predicted in ( Kaliteevski 



et al. 2009). Experimental evidence of the strong cou- 



pling of a Tamm plasmon with quantum well excitons in 
a single DBR mirror-metal/air gap microcavity was re- 
ported in (Grossman^TF^II [2011 Symonds et al. 20091. 
Using metallic layers also allows for a dynamic tuning of 
the external potential by applying DC electric fields to 
the structure: an electric field indeed acts on polaritons 
by reducing the quantum well oscillator strength by spa- 
tially separating the electron/hole pairs and by lowering 
its energy via the quantum confined Stark effect (Kim 



2011) 



et al. . 2008 ) and thermal expansion ( Grossmann et al. 



A technique to generate a periodic potential using sur- 
face acoustic waves (SAWs) propagating in the semicon- 



ductor structure was developed in ( de Lima et al. . 2006 ) 



and then applied to fragment a polariton condensate into 
an ar ray of strongly elongated gas es with reduced coher- 
ence ( Cerda-Mendez et al. 2010). Microscopically, the 



SAW strain field simultaneously changes the thickness 
and the refractive index of the cavity layer, modulating 
the exciton energy via a deformation potential. Each of 
these effects is proportional to the real part of the local 
amplitude of the sound wave. As a result the effective 
potential moves in space at a speed equal to the SAW 
phase velocity, on the o rder of the sound speed in the 
material (3 x 10'^ m/s in Cerda-Mendez et al. 2010), i.e. 



to create polariton boxes with an much smaller than all other characteristic speeds of the 
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polariton system. 

An even more flexible all-optical technique was devel- 
oped in (Amo et al. 2010b): a strong ct_ laser field is 



used to inject cr_ po aritons into the cavity with a given 
spatial profile. Assuming that one can neglect spin-flip 
processes and that a large number of (t_ polaritons is in- 
jected, polariton-polariton interactions between counter- 
polarized polaritons create an effective potential for cr_|_ 
polaritons, whose geometry and time-dependence can be 
easily controlled via the a_ laser fleld. As it is shown in 
the bottom panel of Fig. [5] quite localized potentials can 
be gener ated by this techniqu e. For instance, they were 
used in ( Sanvitto et aZ.[[2011[ ) to study vortex dynamics 
in polariton superfluids. With respect to SAWs, optical 
potentials have the advantage that they can be modu- 
lated on extremely short time scales and designed with 
almost any desired spatial shape down to the /im scale. 

Even though the external potential generated by all 
these methods is mostly a spin-independent one, in some 
cases it is essential to take into account in the model its 
spin-dependent component: the different reflection am- 
plitude of TE/TM polarization states at the cavity/air 
interface in a micropillar device is responsible for a split- 
ting of the two linear polarized states ( Kuther et al.\ 
[1998); an analogous effect may also arise from the cou- 
pling of the exciton with mechanical stress in a laterally 
patterned device (Dasbach et al. 2002). In some pro- 



classical evolution equations for the expectation values 
of the quantum field operators fo-(r) = {^c,cr{'^)) and 
Vai'i^) = (^'x,CT(r)): the former is (approximately) pro- 
portional to the expectation value of the in-cavity electric 
field, the latter is proportional to the matter polariza- 
tion due to the exciton transition. The mean-field equa- 
tions of motion for 5o-(r) and T^o-lr) are obtained from the 
Heisenberg equations of motion for '^c..u{t^) and '^c,a{Y) 
by replacing every instance of an operator with the cor- 
responding expectation value. 

Of course, this approach is exact as long as we are 
restricting ourselves to terms in the Hamiltonian that 
involve at most two operators, e.g. the free dynamics of 
the coupled cavity photons and excitons, their mutual in- 
terconversion, as well as the one-body pumping and loss 
terms discussed in the previous section. An approxima- 
tion is instead made on the interactions terms, for which 
one is assuming, 

(*x,.W^x,.(r)*x,.(r))« 

^{^^x,M){^xA^)){^xA^))- (35) 

In nonlinear optics, a mean-field approximation of this 
kind is implicitly made whenever one writes the polariza- 
tion of a medium as the product of a nonlinear suscepti- 
bility x'") times the n-th power of the classical electric 



posed applications such as (lUmucahlar and Carusottol field (|Boyd| |2008[ [Butcher and Cotter[ |2008[ ): this ap- 



2011), it is essential to be able to impose a specific form 
of spin-dependent potential to the polaritons. The Zee- 
man shift under a static magnetic field can be used to 
split the a± components of t he exciton. A (possibly spa- 



tially modulated) strain field (van Doom et al. . 1996 ) or a 



lateral patterning with submicrometer periodicity ( Flan- 
ders 1983[ ) can be used to generate a linear birefringence 



proach provides accurate results in most optical media as 
nonlinear effects require the presence of a large number of 
photons. If one is interested in the interplay of the opti- 
cal nonlinearity with quantum fluctuations, a linearized 
treatment is generally enough to account for the rela- 



tively weak quantum fluctuations (Walls and Milburn 
2006). An alternative, semiclassical approach is based 



in the cavity material, which would then result in a spin- 
dependent photonic potential of the form 

Fj:t^)(r) = A(r) [cos(20(r))ai;Usin(2^(r))a(^j, 

(34) 

where the angle 0(r) fixes the orientation of the bire- 
fringence axis and the basis of circularly polarized states 
a = a± la used. 



III. THE DRIVEN-DISSIPATIVE GROSS-PITAEVSKII 
EQUATION 

A. The mean-field approximation 

In the previous section we have introduced a quan- 
tum mechanical model for the coupled quantum field 
dynamics of cavity photons and quantum well excitons: 
the interaction terms describing for exciton-exciton col- 
lisions make the dynamics non trivial and are responsi- 
ble for a number of nonlinear and quantum phenomena. 
A standard approximation to attack this kind of prob- 
lems is the so-called mean field approximation, based on 



on the truncated- Wigner representation of the quantum 
field, as discussed in the next subsection. The quantita- 
tive validity of these approaches has been tested by com- 
paring with exact calculations for the simplest model of 
single-mode nonlinear cavity (Carusotto 2001 1 as well as 
with multiniodc Wigner simulations ( Verger et al. . 2007 1 . 



Of course, this approximation breaks down completely as 



soon as one enters the photon blockade regime ( Imamoglu 



et al. . 19971 where a single photon is able to substantially 



modify the response of a device and the discrete, quan- 
tum nature of photons starts playing a crucial role. 

In the context of quantum gases of material parti- 
cles, a classical partial differential equation for the su- 
perfiuid order parameter was first written by |Gross]|1961| 
and IPitaevskli] |1961| to describe quantum vortices in 
liquid Helium. Starting from modern formulations of 
the Bogoliubov theory of the dilute Bose gas, a Gross- 
Pitaevskii equation (GPE) of the form 



2m 



I^-P* (36) 



can be derived from first principles to describe the 
dynamics of the condensate wavefunction ^'(r,i): the 
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macroscopic fraction of particles that populate the con- 
densate mode behave in a collective way and the quantum 
atomic matter field ^{r,t) behaves as a classical field 



^'(r,i) (Lcggett 2001 Pitaevskii and Stringari 2004 1. 



Generally speaking, the GPE is quantitatively accurate 
as long as the occupation of modes other than the con- 
densate one is small: in three dimensions this requires 
that the temperature is much lower than the transition 
temperature for Bose-Einstein condensation and that 
the gas is dilute, i.e. the atom-atom scattering length 
a is much shorter than the mean interparticle spacing 
na^ < 1. 

From a physical standpoint, the GPE is then for the 
matter field what Maxwell equations are for quantum 
electrodynamics in nonlinear media. An important dif- 
ference is however worth emphasizing: while the global 
phase of the electromagnetic field has a direct and ob- 
servable physical meaning, conservation of the total num- 
ber of particles implies that the expectation value of the 
matter field of massive particles exactly vanishes. Even if 
treatments based on a spontaneous breaking of the phase 
symmetry are popular in the literature and have a deep 
physical foundation in terms of the BEC phase transi- 
tion ( Gunton and Buckingham 1968j ) , some doubts may 
remain on the consistence of descriptions based on a clas- 
sical atomic field, especially for finite systems. The pic- 
ture has been reconciled by recent theoretical works that 
have developed particle-number conserving versions of 
the Bogoliubov theory (I Cast in and Dum 1998 Gardiner 



[T997) and have investi gated in full detail the meaning of 



the condensate phase I 


Castin and Dalibard 


1997 


Cirac 


et al. 1996 


Javanainen and Yoo 


1996 


Leggett 


1995 


). 



B. Generalized GPE under coherent driving 

In many experimental circumstances, it is not neces- 
sary to work with the pair of equations of motions for the 
photonic £o-(r) and excitonic 'P<j(i") fields and one can re- 
strict to a single classical field describing the lower polari- 
ton field 5'ip(r,t) in a single spin state. This simplified 
description is generally legitimate provided the Rabi fre- 
quency ^Ifi is much larger than all other energy scales of 
the problem, namely the kinetic and interaction energies, 
the pump detuning from the bottom of the lower polari- 
ton, and the loss rates. The resulting Gross-Pitaevskii 
equation for polaritons has the form 



iat*Lp(r,t) 



2mLp 

+ 5Lp|*Lp(r,i)P^'Lp(r,t) 



V^ + VLp{r) *Lp(r,t)+ 



^Lp{r,t)+irjLpE'''^r,t). (37) 



With respect to the Gross-Pitaevsku equation (36) de- 
scribing atomic Bose-Einstein condensates, Eq.(37) in- 
cludes additional terms to account for the driven- 
dissipative nature of the polariton gas, namely a loss rate 



proportional to 7lp and the coherent pumping propor- 
tional to the incident field i?'"'^(r, t). 



The parameters appearing in the polariton GPE (37) 
are defined in terms of the microscopic description of 
polaritons discussed in the previous section: niLP — 
TOcav/lw^^P is the lower polariton mass, w£p is the fre- 
quency of the bottom of the lower polariton branch, 
l^Lp(r) is the external potential felt by polaritons as de- 



fined in (33), gLP is the polariton-polariton interaction 



constant introduced in (31 ) and the t/lp parameter quan- 
tifies the coupling of the polariton to incident radiation, 
Vlp — Va ■ Of course, an identical equation holds 
for planar cavities embedding a nonlinear medium: in 
this case, the polariton reduces to a pure photon and the 
interaction constant gi^p is related to the nonlinear sus- 
ceptibility x'^^ of the medium. In the nonlinear optics 
literature, this equation often goes under the name of 
Lugiato-Lefever equation (Lugiato and Lefever 1987). 



C. Quantum fluctuations and Wigner representation 



Going beyond the Gross-Pitaevskii equation (37) and 
fully include the quantum and thermal fluctuations of the 
quantum field 5'lp is a very difficult task in the general 
case. In the course of the years, a number of different 
methods have been developed to attack this problem in 
different regimes, for instance Path Integral Monte Carlo 



for the thermal equilibrium state (Ceperley 
agrammatic techniques ( [Houcke et al. 2010 



1995), di 



Prokof 'ev 



and Svistunov 2008), density matrix renormalization 
group techniques fo r the both the gro und state and the 
temporal dynamics ( Schollwock[ 2005). 

In this section we will then limit ourselves to a short 
review of phase space techniques originally developed in 
the quantum optical context and recently applied with 
success to the study of quantum fluids of atoms and 
photons. The basic idea is to represent the state of 
the quantum field as a (non-necessarily positive) prob- 
ability distribution function on a suitable classical phase 
space and describe the time-evolution of the field in terms 
of a Fokker-Planck-like partial differential equation or a 
stochastic partial differential equation. A general intro- 
duction to the basic concepts of phase space representa- 



tions of a quantum field can be found in (Gardiner and 



Zoller 2004) 



A crucial feature in view of efficient numerical simu- 
lations is that the probability distribution be positive- 
valued and the time-evolution be described by a Fokker- 
Planck partial differential equation (Risken 1996) 



dP{x,t) 
dt 



M 



d 



z— 1 ^ 



with a drift force F and a positive-definite diffusion ma- 
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trix D. The probability distribution P(x, t) is here de- 
fined on a M-dimensional space with real coordinates 
X;, but probability distributions for complex quantities 
are straightforwardly included by considering the real 
and the imaginary parts as independent real variables. 
Provided the diffusion matrix D is positive-definite, the 
Fokker-Planck equation ([38| can be mapped onto a sys- 



tem of M Ito stochastic differential equations ( Gardiner 



2009 ) of the form 



dxi = Fi{-K,t) dt + dWr 



(39) 



the incoherent pumping case requires some modeling of 
the relaxation mechanisms in the device: first attempts 



in this direction were recently reported in (Wouters and 



Savona 20091. Typical implementations of the Wigner 



representation require discretizing the quantum field on 
a d dimensional discrete lattice of N'^ points enclosed in 
a finite integration box of side L. In this geometry, the 
Wigner distribution is a function of the M'^ complex am- 
plitudes ipi = '0(ri) of the field at the lattice positions r^. 
Its time-evolution is described by the Fokker-Planck-like 
equation 



with a Wiener noise satisfying 



dW^dWj = D,j{yi,t)dt, 



(40) 



that can be efficiently simulated on a computer by taking 
the statistical average over many different realizations of 
the Brownian motion. 

Unfortunately, in most cases of actual interest this 
is not possible as either the diffusion matrix D is not 
positive-definite or additional terms with higher order 
partial derivatives of P are present in the right-hand side 
of (38): for the system of interacting bosons under in- 



vestigation here, non-positive diffusion terms appear in 
the time-evolution of the Glauber P and Husimi Q rep- 
resentations, while third order derivative terms appear in 



the time-evolution of the Wigner W (Vogel and Risken 



19881 ). 

A possible way out based on the so-called Positive- 



P representation was proposed in (Drummond and Gar- 
|dincr^ 1980^,Drummond and Walls[|T980) and transferred 



to the realm of atomic condensates in ( Steel et al. 1998 1 . 



A remarkable application of this representation to multi- 
mode optical systems addressed squeezing of quantum 
solitons propagating in nonlinear optical fibers (Carter 
[1987| ). 



et al. 



Unfortunately, in many other relevant 
cases the stochastic differential equations in the dou- 
bled phase space show strongly divergent trajectories 
that intrinsically undermine the numerical stability of 

The use of improved 



the method (Gilchrist et al. 



schemes based on the so-ca 



1997) 



was proposed in (Deuar and Drummond 2002 Drum 



led gauge-P representation 



mond et al. 2004 ) and has provided interesting results 



on one-dimensional interacting Bose gases (Deuar et al. 



2009). Another exact stochastic approach to the many- 



body problm in (Carusotto and Castin 



2001 Carusotto et al. 2001 ) using a Hartree ansatz in the 



canonical ensemble: while the method could be proven 
to be mathematically consistent and to have a finite sta- 
tistical error at all times, useful results could only be 
obtained in simple geometries with a limited number of 
particles ^Carusotto and Castin 2003). 



In the following of the section we shall focus our at- 
tention on the Wigner W representation that has turned 
out to be most useful in practical calculations. For the 
sake of simplicity, we restrict to the coherent pumping 
case for which a self-contained Hamiltonian description 
is available (Carusotto and Ciuti 2005). Extension to 



dW 

~dt 



+ 



2 AV d^* d^p. 



4AF2 dlpi 



d 



d-ipi 



(41) 



where AV — (L/Af)'^ is the volume of the elementary cell 
of the discrete lattice. The drift force term on the site 
Fi{tp} = F{-0}(r = Ti) involves a deterministic evolution 
of the field very similar to the right-hand side of (37), 



F{^}(r) - 



■"LP 



Hlp 



V2 + Fip(r) 



+ ffip(|V(r,i)|2-— ) V(r,t)+ 



*77Lpi?"^(r,t). (42) 



The derivatives with respect to the complex variable 
ip are defined as 

d 



_§_ 



d 



d 

dip* 



dRc[jp] + *Olm 

is always positive and can be straightforwardly mapped 
onto a noise term with local correlations in space. On 
the other hand, the third derivative terms can not be 
included in a standard stochastic differential equation 
of the form ( 39 1 . Recent works ( Plimak et al. . 2001 ) 



The second-order derivative term 



have tried to solve this issues by generalizing the concept 
of stochastic differential equation to the case of stochas- 
tic difference equations with discrete time steps of finite 
size At: provided one does not attempt to take the limit 
At — > 0, a stochastic process can be found whose average 
over the noise recovers all terms of the Fokker-Planck- 



like equation (41 ). Given the large statistical error of the 



method, the application of this Positive- Wigner method 
has been so far limited to very simple, few-mode models. 
Another strategy was proposed in (Polkovnikov |2003 ) 
in terms of a perturbative expansion in the number of 
quantum scattering events accounting for the third-order 
derivative terms 

The main practical interest of the Wigner represen- 
tation follows from the different scaling of the various 
terms of (41 1 in the dilute gas limit — > oo, g^p 
at a constant interaction energy gLp\'ip\'^. In this limit. 
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the noise is responsible for a statistical fluctuation of 
the field ip around its mean-field value on the order of 
Atp (X AV^^^^. Correspondingly, the characteristic mag- 
nitude of the third-order derivative term is roughly esti- 
mated on the order of g^p/AV to be compared to the 
diffusion term which is of the order of 'Jlp- Provided 



finite temperature can be found in (Davis et al. 20121; 



an interesting application of the stochastic GPE (Gar- 



diner and Davis 2003) to the formation dynamics of an 



atomic condensate was discussed in ( Weiler et al. , 2008 1 . 



7LP > gLp/AV, 



(43) 



one may then expect that accurate results for a non- 
equilibrium gas of photons can be obtained by means 
of truncated Wigner calculations where the third-order 
derivative term is completely neglected and the stochastic 
partial differential equation has the form; 



dV(r,t) = F{iP}{r)dt- 



ILP 

AAV 



dW{T,t) (44) 



with a zero-mean, complex Gaussian noise term dW sat- 
isfying 

dW{Y,t)dW{r' ,t) = (45) 
dW*{Y,t)dW{r' ,t) = 2dtSr,r'- (46) 



D. Incoherent and saturable pumping 

The discussion in the previous subsections addressed 
the case where the microcavity is driven by a coherent, 
quasi-resonant pump: in this case, the microscopic de- 
tails of the system are under control and one can de- 
velop an ab initio description of the system. Many re- 
cent experiments were performed under different pump- 
ing schemes where the polariton gas does not inherit any 
coherence from the pump and/or coherence is quickly lost 
during the relaxation process towards the bottom of the 
lower polariton branch. This is the case of the incoherent 
and/or far blue-detuned laser pump used in (Kasprzak 
Richard et al. 2005a|b ), the large angle op- 



et all 



2006 



tical drive of (Deng et al. 2003), as well as, more re- 
cently, the electrical injection of polaritons of (Bajoni 



Under the condition (43), 
that are introduced by the third-order derivative term 
are quickly washed away by the losses and the classical 
noise associated to jlp- in particular, the magnitude of 
the classical noise is the "right" one to simulate in the 
most accurate way the quantum dynamics on a classical 
computer. In (Sinatra et al. 2002), it was shown how 



the non-classical correlations etal. .2008 Khalifa ei aZ. 2008 Tsintzos aL 20091. 



A detailed study of the microscopic mechanisms in- 
volved in the kinetics of incoherent pumping schemes can 
be found in (Porras et al. 2002 ): Typically, all these tech- 



the truncated Wigner approach is able to capture the 
quantum fluctuations at least at the level of Bogoliubov 
theory. 

For typical experimental parameters of dilute photon 
gases in planar devices, the condition imposed by (43) 
on the lattice spacing appears to be fully compatible 
with the characteristic range of wavevectors that are in- 
volved in the physics of the system: a first application of 
the truncated Wigner method to the critical fluctuations 
across the parametric oscillation threshold of a planar 



niques end up accumulating a quite significant density of 
incoherent polaritons in the so-called bottleneck region in 
k-space in the vicinity of the inflection point of the lower 
polariton dispersion: further relaxation to the bottom of 
the lower polariton band via phonon-polariton scattering 
is in fact slowed down by the reduced density of states of 
the final states (the use of a finite density of electrons in 
the quantu m well to overcome this i ssue was theoretically 
proposed by Malpuech et al. 2002 ) . At high enough po- 



semiconductor microcavity was reported in (Carusotto 



and Ciuti 20051 and is reviewed in Sec V.C.2 More re- 
cently, the truncated Wigner method has provided the 
first quantitative evidence of analog Hawking radiation 
from acoustic black holes in flowing polariton fluids ( Ger- 



ace and Carusotto, 2012). Remarkably, a careful appli- 



cation of the truncated Wigner method provides useful 
results also in the conservative case of atomic gases for 
which 7lp = and the (sufficient) validity condition (43) 
is not fulfilled. In this case, the truncated Wigner equa- 
tion is a purely deterministic one and the quantum fluc- 
tuations enter the model via the randomness of the field 
Tp at the initial time of the simulation. A detailed char- 
acterization of the power and the difficulties of the trun- 
cated Wigner method to study the thermal equilibrium 
state and th e dynamical properti es of atomic gases can 
be found in (Sinatra et al. 2002). A general review of 



lariton densities (i.e. for strong enough pumping intensf 
ties), another relaxation mechanism based on polariton- 
polariton collisions becomes active: a pair of polaritons 
in the bottleneck region collide and are respectively scat- 
tered to the bottom of the LP branch and to the large 
wavevector region where polaritons have a mostly exci- 
tonic nature with a very large density of states. Given 
the bosonic statistics of polaritons, this relaxation pro- 
cess turns out to be stimulated as soon as the phase-space 
density of polaritons already present at the bottom of the 
lower polariton starts being of order one: when stimula- 
tion overcomes losses, a macroscopic coherent population 
of polaritons accumulates in the final state and a conden- 
sate appears. On the other hand, the phase-space den- 
sity on the excitonic branch always remains much smaller 
than one, which guarantees effective irreversibility of the 
scattering process. 

A phenomenological description of the condensate dy- 
namics under an incoherent pumping inspired to the 
semi-classical theory of laser (Lamb 19641 was intro- 
duced in (Wouters and Carusotto 2007a I: stimulated 



a wider class of C-field methods for simulating the non- 
equilibrium dynamics of degenerate Bose gases at zero or 



scattering into the condensate is described by an ampli- 
fication term in the field equation of motion, with an 
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amplification rate rate which is a monotonically grow- 
ing function R[n}i\ of the local density of the polariton 
reservoir nfi{r) in the bottleneck region. In practice, this 
corresponds to adding one more term to the polariton 
GPE (ISTl) of the form 



et al. 



1998 ) in the completely different context of the the- 



idt-^Lpir, t) = ... + -R[nR{r, t)] ^Lp{r, t) (47) 

and describing the polariton reservoir density via a rate 
equation of the form 

dtnniv, t) = P- i?K(r, t)] |*Lp(r, i)P - 7i?. nR{v, t). 

(48) 

The intensity of the incoherent pumping is phenomeno- 
logically described by the pump term P{y). In the typi- 
cal case where the characteristic relaxation rate ^r of the 
reservoir is much faster than all other scales ^, one can 
adiabatically eliminate ur. In the simplest case where 
R[nR\ = RriR, the amplification term then reduces to 
the form 



■"LP 



2mLp 



V^ + VLp{r) *ip(r,t) 



llLP 



+ 9LP |*Lp(r, t)\' *Lp(r, t)-^ *Lp(r, t)+ 
iP(r,t)/2 ^ , , 
7_R + i?|1'Lp(r, i)r 

that closely resembles the complex Ginzburg-Landau 
equation currently used in the theory of pattern for- 
mation in nonlinear dynamical systems (Aranson and 



Kramer[ |2002[ [Cross and Hohenberg[ |1993[ ). For this 



reason, some authors prefer to directly call it the com 
plex Ginzburg-Landau equation (cGLE) for the polari- 
ton condensate. A slightly different, but almost equiva- 
lent form of (|49]) was independently introduced in ( Keel- 
ing and Berloff, ,2008, ) . A generalization of the ( 49 ) in- 
cluding the energy-dependence of the amplification term 
was proposed in (Wouters and Carusotto 2010) to ac- 



count for the stronger scattering into the lowest polariton 



2009) 



states (Sarchi and Savona 2008 Wouters and Savona 



The phenomenological nature of this model allows for 
its direct application to a wide class of incoherent pump- 
ing schemes: the details of the specific configuration enter 
via the functional form of the scattering rate R and the 
other parameters of the theory: the only requirement is 
the presence of an irreversible process that is able to in- 
ject extra polaritons in the lowest energy states and that 
can be stimulated by pre-existing polaritons. It is worth 
mentioning that a very similar model was used in ( Kneer 



^ Note that the relaxation rate 'yji does not need to correspond 
to the actual decay rate of reservoir polaritons, but may simply 
account for their redistribution among the different states of the 
reservoir. 



ory of atom lasers, still-to-come devices that should pro- 



duce a coherent atomic matter wave as an output ( Bloch 



et aL[|1999[|Chikkatur et aL[ |2002[ [Lahaye et~al]\200^ 



Extension of this model beyond the mean-field approx- 
imation to include quantum fluctuations of the polariton 
field requires a more sophisticated model of the ampli- 
fication mechanism where the reservoir degrees of free- 
dom are consistently taken into account. A first at- 
tempt in this direction based on the Wigner representa- 
tion of the polariton field was reported in (Wouters and 



Savona 2009). A completely different approach based on 



a Keldysh diagram solution to a non-equilibrium Fermi- 
Bose model was proposed in ( Szymahska et al. 2006 ) and 



provided interesting predictions for the coherence prop- 
erties of a polariton condensate. 



IV. POLARITON CONDENSATES UNDER 
QUASI-RESONANT EXCITATION 

In the previous section we have introduced the general- 
ized Gross-Pitaevskii equation that describes the dynam- 
ics of the polariton condensate at the mean-field level. 
Now we shall proceed to review its application to derive 
the stationary state and the spectrum of the elementary 
excitations of the condensate under the different pump- 
ing schemes: The present section is devoted to the case of 
a coherent pumping at the bottom of the lower polariton 
branch, the next section is devoted to the so-called para- 
metric regime with a coherent pump in the vicinity of the 
inflection point of the lower polariton branch, while the 
following section is devoted to the case of a generic inco- 
herent pumping. To facilitate the reader in appreciating 
the novel features that stem from the non-equilibrium, 
driven-dissipative nature of the polariton condensate, we 
shall start with a short survey of the basic properties of 
standard condensates at equilibrium, e.g. in ultracold 
atomic gases. 



A. Equilibrium condensates 

1. The condensate wavefunction 

The textbook definition of Bose-Einstein condensation 
in an ideal Bose gas is that a macroscopic fraction of the 
particles occupy the lowest single-particle energy state: 
in free space, this is the k = plane wave, while in a 
trap it is the lowest eigenstate of the trap potential. In 
presence of weak particle-particle interactions, the con- 
densate wavefunction 4'o(j) is obtained as the lowest fi 
eigensolution of the time-independent Gross-Pitaevskii 
equation (Pitaevskii and Stringari 2004) 



2m 



V^Mr) + ^trap(r) Mr)+9\Mr)\^^oir), 

(50) 
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where m is the mass of the particles, Vtrap(r) is the trap 
potential, and g = Airfi^a/m is the interaction constant in 
terms of the two-body scattering length a. A few general 
features of this equation are worth mentioning. 

In agreement with gauge symmetry, the global phase of 
the wavefunction 0o(r) is not fixed. As it has no physical 
consequence, it can be arbitrarily chosen. The eigenvalue 
fi gives the oscillation frequency of the condensate phase 
under the time-dependent GPE ([36]). In a thermodynam- 
ical context, it corresponds to the chemical potential of 
the condensate, i.e. the energy that is required to add 
one more particle to the gas. In free space (Vtrap — 0), 
the equation of state relating the chemical potential and 
the (spatially uniform) density |0oP has the simple form 



(51) 



To the best of our knowledge, the local phase of 0(r) is 
constant throughout space in all known cases even in the 
presence of a trap potential Vtrap, which means that 0o(r) 
can be reduced to a purely real function ^. Physically, 
having a purely real (/'o(r) means that no macroscopic 
current is flowing across the condensate in its ground 
state. 



1 

2mi 



[0:(r)V0o-li.c.] = 



(52) 



As we shall see in what follows, polariton condensate 
show very different features: the system is far from a 
thermodynamical equilibrium state and the macroscopic 
occupation of a single state ^'Lp(r) can not be derived 
from a free energy minimization argument. As a result, 
the shape of ^'Lp(r) is determined by a complex dynami- 
cal balance of pumping and losses: depending on the geo- 
metrical configuration, stationary states with a spatially 
dependent phase can be found, which physically corre- 
sponds to a non- vanishing polariton current J^p along 



the cavity plane. For instance, polariton fluids with a 
well-defined and non-zero in-plane momentum k can be 
created using a coherent pump with a non- vanishing inci- 
dence angle 0i„c • A radial flow in a polariton condensate 
under a spatially localized incoherent pump was invoked 
in ( Wouters et al. . 2008 1 to explain the ring-shaped mo- 



mentum distribution observed in ( Richard et al. 2005b I ; 



a similar mechanism due to polariton-polariton interac- 
tions underlies the macroscopic flow in a one-dimensional 
geometry in photonic wires (Wertz et al. 2010). Quan- 



tized vortices are spontaneously present in the steady 
state of a polariton condensate under a incoherent pump 



2008). 



in the presence of disorder (Lagoudakis et al. 
2. Bogoliubov dispersion of elementary excitations 



The dynamics of weak perturbations on top of a di- 
lute Bose-Einstein condensate can be described within 
the Bogoliubov theory (Castin 2001 Pitaevskii and 



Stringari, 2004): weakly excited states of the condensate 



are characterized by bosonic excitation modes, whose fre- 
quency and spatial profile can be obtained by lineariz- 
ing the Gross-Pitaevskii equation around the equilibrium 
state. 

A weak modulation of the condensate wave- 
function can be represented as 0(r, t) = 
[<j)oir,t) + S(t){r,t)] e-^t'*/''' with small (5(/)(r,<). In- 
serting this ansatz into the GPE (36) an linearizing 



around the stationary solution leads to a pair of evolution 
equations 



ihdt 



Sr{r,t) 



= C 



Bog 



<50(r,t) 
Sr{r,t) 



(53) 



with the Bogoliubov operator defined by 



C 



Bog 



2m 



-gn 



2gn — fj. 



2m 



gn 
Kxt(r) 



2gn + fi 



(54) 



In the simplest case of a spatially homogeneous system, 
the equilibrium condensate wavefunction is constant in 
space (j)o{j^,t) = \fn and the Bogoliubov modes can be 
chosen in a plane-wave form of wavevector k. In this 
case, the Bogoliubov operator reduces to a matrix and 
its eigenvalues give the so-called Bogoliubov dispersion 
of excitations. 



fi^Bog(k) = ±^— (— +2^j, (55) 
the + (— ) sign referring to the so-called positive (neg- 



ative) Bogoliubov branch, i.e. the modes with positive 
(negative) Bogoliubov norm. As it is discussed in full de- 
tail in textbooks on BEG, the energy of the physical ele- 
mentary excitations is determined by the positive branch 
only (IGastinl 120011 IPitaevskii and Stringaril 12004 1. 



Depending on the value of the wavevect or k as com- 
pared to the so-called healing length ^ = yJWJm^,, two 
regimes can be identified in the Bogoliubov dispersions: 
for large momenta fc^ ^ 1, the Bogoliubov dispersion 
recovers the parabolic shape of single particles with a 
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global energy shift due to the Hartree interaction energy, 



a;Bog(k) 



2m 



(56) 



On the other hand, small wavevector A:^ ^ 1 excitations 
have instead a phononic nature with a sonic dispersion 



WBog(k) ~ Csfc 



(57) 



with a sound speed = ^ ^i/m. The presence of a sharp 
corner at k = crucially depends on the relative value of 
the non-diagonal elements ign of Z^Bog and the Hartree 
shift ±{2gn — fx) = ±gn in the diagonal ones: physically, 
it can be seen as a remarkable consequence of the dy- 
namical stability of the system combined with the global 
gauge symmetry under rotations of the global phase of 
the wavefunction. Experimental studies of the elemen- 
tary excitation spectrum of liquid Helium showing both 
the phonon and the roton branches can be found via 



neutron scattering experiments (Palevsky et al. 1958 



Yarnell et al. 



19591. Recent studies of the Bogoliubov 



dispersion in dilute atomic conden sates using two-ph oton 
Bragg scattering are reviewed in (Ozeri et al. 2005). 

As we shall see in Sec |IV.C| breaking this condition in 
the driven-dissipative case leads to several novel phenom- 
ena in the excitation spectrum of polariton condensates. 



B. Stationary state under coherent pumping 

In this subsection, we review the steady-state 
properties of polariton condensates under a coherent 
continuous-wave pump at frequency Wjnc with different 
spot geometries and, for the sake of simplicity, consider- 
ing regimes where the form (371 of the generalized po- 
lariton GPE can be applied. 



1. Optical limiting and optical bistability 

When the cavity is driven by a coherent, continuous- 
wave pump with a plane-wave spatial profile, 



£;'""(r,t) = 



mc g^kii 



(58) 



we can 



LP' 



look for solutions like ^Lp(r,t) — 
c-»i^inct Inserting this ansatz into the 



generalized GPE (|37[), we obtain the polaritonic analog 



of the equation of state. 



2mLp 



-9LP\^1 



nLP 



LP\ 



^ LP 



iriLpEZ. 
(59) 

relating the polariton field amplitude ^^.p ^'^ pump 
parameters. Comparing this equation with the standard 
equation of state (51), we notice that the gauge sym- 



metry under global rotations of the condensate phase is 
explicitly broken by the coherent pumping term on the 
right-hand side of ( 59 ) . The consequences of this feature 



on the elementary excitation spectrum will be discussed 
in the Sec llVTil 

Furthermore, while the oscillation frequency of the 
matter field in standard equilibrium condensates is re- 
lated to the density via the chemical potential, the oscil- 
lation frequency of the polariton condensate wavefunc- 
tion LP is fixed by the pump frequency Winc- On 
the other hand, the dependence of the polariton den- 
sity riLp = I^'lpP on the pump parameters is strongly 
affected by the resonance condition between the pump 
frequency Winc and the polariton frequency cjip(kinc) — 
uj1p+hkf^^/2mLp, possibly shifted by the interaction en- 
ei'gy 5lpI^lpP- As a result, the behavior of the density 
riLp as a function of the pump intensity Ilp = l-E-incP is 
very different depending on whether ujinc is on the blue 
or on the red side of ojLpj^inc)- This is a well-known fact 
of nonlinear optics ( |Boyd 2008 ) and is illustrated in the 
left panels of Fig|6] 

The Wine < '^Lp(kinc) case is shown in the upper-left 
panel: the density ulp is a continuous, monotonically 
growing bmction of the pump intensity Ilp. As the non- 
linear term tends to shift the polariton frequency further 
away from resonance with the pump, the growth of the 
density is sub-linear, a behavior known under the name 
of optical limiter. 

The Wine > WLp(kine) case is shown in the lower-left 
panel: the density ulp is still a growing function of Ilp 
but shows a hysteretic behavior that goes under the name 
of optical bistability in the nonlinear optics literature. 
The central branch of the hysteresis loop with a nega- 
tive slope is always dynamically unstable. Depending 
on the specific values of pump parameters, also other re- 
gions may be dynamically unstable towards the paramet- 
ric generation of polaritons into other modes at different 
k ^ kinc. These stability issues will be addressed in full 
detail in Sec jlV.C] using the polariton version of the Bo- 
goliubov theory. 



2. Local density approximation and experiments 

Optical bistability behaviors have been observed 
in several experiments with polaritons in microcavi- 
ties (Baas et al. 2004a): while there is good qualita- 



tive agreement with the theory discussed in the previous 
subsection, a quantitative understanding of the peculiar 
geometrical features observed in the experiments requires 
a more refined theoretical treatment accounting for spa- 
tial inhomogeneities. Yet, the equation of state (59) is a 



good starting point to describe the polariton wavefunc- 
tion even in this case. Provided the spatial dependence 
of the pump and of the external potential is slow enough, 
one can approximate the local polariton field at each 
spatial position with the homogeneous system prediction 
( 59 ) using the local value of the pump amplitude and in- 



serting the external potential yLp(r) as a shift of a;£p. 
This procedure of neglecting the kinetic energy term in 
the GPE ( 37 ) provides an accurate approximation of the 
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FIG. 6 Left panel: mean-field energy (meV) for the exciton field versus pump intensity in the optical limiter regime (top) and 
bistable regime (bottom). The dashed line indicates the unstable regions. Right panel: real part (top) and imaginary part 
(bottom) of the excitation frequencies for the non-equilibrium Bogoliubov modes corresponding to the points A,B,C,C',D,E, 
indicated in the left panel. 



time- independent GPE of trapped atomic gases ( 50 ) and 
leads to the so-called Thomas-Fermi profile (Pitaevskii 



and Stringari 20041 



In the polariton case, one has to consider the additional 
features due to the optical bistable behavior: for a given 
value of the pump intensity, several solutions for the po- 
lariton density may be available. An intuitive criterion 



to choose the relevant solution was introduced in (She- 



lykh et aL[|2008[ ) and successufully compared to complete 
simulations of the GPE: while moving in space along a 
streamline (i.e. along the direction of the wavevector 
kinc), one has to continuously follow a given branch of 
the hysteresis loop. Jumps to other branches are only 
allowed when one reaches the end point of the branch. 
Optical multistability effects in the presence of spin de- 



(Gippius et al. 


2007 


(Paraiso et al. 


20101 



Interesting spatio-temporal switching effects when polar- 
ization multistability is combined with a non-trivial spa- 



tial dynamics have been theoretically studied in (Liew 
et al.\ |2008 ) and experimentally demonstrated in ( Amo 
eTon 2010al). 



3. Spatially localized coherent pump 

A completely different regime is realized when the 
monochromatic pump is focused onto a spatially small 
region, with a e.g. Gaussian circular spot of peak ampli- 
tude £'g"'^, radius a and carrier wavevector kjnc = kiac^x- 
In the regime of weak pump intensity where nonlinear in- 
teractions are negligible, analytic calculations show that 
the steady-state polariton field 



(60) 



extends in the direction of the pump wavevector for the 
distance i = Vg/^f that cavity photons ballistically cover 
in a time 7~^. On the other hand, the polariton field 
follows the profile of the pump spot along the orthogonal 
direction y and keeps a narrow width a. The small spot 
assumption is consistent if cr ^, while the wavepacket 
expansion due to the curvature of a;cav(k) is safely negli- 
gible provided a ^ yjhj^hrv^. 

Remarkably, the wavevector k of the ballistically mov- 
ing polaritons in ( 60 1 is not fixed by the incident laser 
wavevector kjnc, rather it is determined by the laser 
frequency lo\^c via the energy conservation condition 
a;Lp(k) — Wine. Correspondingly, the group velocity 
Vg = Vkt^LP has to be evaluated at k = k, which gives 
Vg ~ /ik/mLp. The missing momentum k — kjnc is pro- 
vided by the spatial inhomogeneity of the pump spot. A 
similar physics in the presence of interactions for a inco- 
herent pump with a small spatial spot will be discussed 
in Sec lVI.A.3l 

As we shall review in Sec |VIIl[ the fact that the phase 
of the polariton field in the region outside the pump spot 
is not fixed by the incident laser is crucial to observe 
topological excitations in the polariton fluid past an ob- 
stacle: this feature was theoretically put forward in (Pi- 



geon 


et al. 


et al. 




2011 



2011) and experimentally observed in (Amo 



Nardin ei a/. 2011b Sanvitto ei a/. 20111. 



proposal of analog black holes based on polariton conden- 



sates ( jSolnyshkov et al.\ |2011j), ( [Gerace and Carusotto 
^0T2|. 



C. Elementary excitations and dynamical stability 

Once the steady state polariton field under a contin- 
uous wave pump has been determined, we can proceed 
with the study of small fluctuations around the steady 
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state. This can be done by extending the BogoUubov the- 
ory of dilute Bose gases to the present non-equihbrium 
context of polaritons using the generahzed GPE (ISTI) . A 



pioneering mention of the Bogohubov dispersion of quasi- 



particles on top of a luminous fluid appeared in Chiao and 



|Boyce| [19991 [Tanzini and Sorella[ [19991 where an equr 
librium assumption was implicitly made. A first quan- 
titative study taking into account the driven-dissipative 
nature of the polariton fluid was carried out for the co- 



herent pumping case in ( Carusotto and Ciuti 2004 Ciuti 



and Carusotto 2005). In the following of the discussion 



we shall closely follow the discussion in these works. 

For the sake of simplicity, we shall restrict ourselves 
here to the illustrative case of a spatially homogeneous 
system Vlp{y) = under a coherent pump with kjnc = 0; 
other cases will be reviewed in Sec. IV.C.3I and IVI.BI and 
jVILBj Under these assumptions, the steady state has the 
form 4'Lp(r,t) = ^n^p exp(— iwi„ct) and the linearized 
GPE for the polariton field modulation (S5'Lp(r,t) reads 



idt 



<5*Lp(r) 
<5*lp(r) 



WLp(k) -f 2gLpnLP - Winc - i'lLp/'^ 

-gLpriLP -WLp(k) 



gLpriLP \ f d^Lp{r) 

2gLpnLP + Wine - mp/2 J \ S'i'lpir) 



^ (61) 
being reminiscent of the equilibrium Bogoliubo v equations (53||54). The additional terms are the loss rate jlp and 

pump frequency Winc next to the Hartree energy (Ciuti et al. 2001): in particular, the real part w^p -I- ig^priLp — Wjnc 



is no longer equal to the non-diagonal term gLpn-LP, leading to a wealth of new behaviors (see the plots in the right 
panels of Fig. l6| described by the analytic formula 



a;Bog(k) = ± (' 



^LP 



he 

2rn 



2gLpnLP 



[gLpnLpY 



1/2 



ILP 



(62) 



As expected, the presence of the coherent pump locking 
the phase of the condensate directly reflects in the ab- 
sence of a Goldstone branch whose frequency uj tends to 
zero for k — )■ 0. 



The optical limiter case ajj, 



< UJ 



LP 



is illustrated in 



panel (E): the two Bogoliubov branches are split in fre- 
quency and the sonic behavior at small k ( 57 1 disappears 



The stronger the pump intensity, the wider the gap. The 
imaginary parts of the two branches are flat and equal to 

-ILP- 



In the optical bistability case 



LP-i 



regimes have to be distinguished. On the upper branch of 
the bistability loop [panel (D)], a behavior similar to the 
optical limiter case is recovered as interactions have effec- 
tively shifted the polariton branch ajip(k) + 2gLpnLp to 
frequencies larger than uJiac ■ the Bogoliubov branches are 
split in energy by a gap that increases for growing pump 
intensity. A sonic dispersion of the form (57) is recov- 



ered only at the end-point C of the upper branch, which 
satisfies the effective resonance Wjnc = ^lp + gLP'iT-LP 
condition. Only at this point the oscillation frequency 
Wine happens to fulfill the analogous of the equilibrium 
equation of state (51 ). 



The dynamical instability on intermediate branch (in- 
cluding point C) is the standard Kerr single- mode insta- 
bility of optical bistability (jCiuti and Carusottol |2005 



Drummond and Walls 1980[), signaled by the imaginary 



part becoming positive in a neighborhood of the pump 
wavevector kinc = 0. 

The physics is more interesting for pump values on the 
lower branch: the two Bogoliubov branches cross on a 
ring of wavevectors [panels (A,B,B')]. Because of the anti- 



Hermitian non-diagonal terms in the matrix (61) that 
provide an effective "attraction" between modes (Ciuti 



et al. 2001 Savvidis et al. 2001), the branches stick in 



the vicinity of the crossing point giving rise to a flat re- 
gion where Re[a;Bog] = 0. On the other hand, the imag- 
inary parts Im[a;Bog] are split: as soon as one of them 
becomes positive, the system becomes dynamically un- 
stable (as at point B') and a parametric oscillation sets 
in, as discussed in Sec[V] It is worth noting the qualita- 
tive analogy of this behavior with the soft mode of an 
several open Dicke model (Dimer et al. 2007) when the criti- 



cal point of the Dicke transition is approached (Emary 



and Brandes 


2003 


Nagy et al. 


2010 


2011 



tal investigations of the Dicke p 



lysics using an optically 



driven atomic condensate embedded in an optical cavity 



have recently appeared in (Baumann et al. 2010 2011). 

In Sec jVII.B we shall discuss in detail the consequences 
of the form ( 62 ) of the Bogoliubov dispersion onto the su- 
perfluidity properties of the polariton fluid and we shall 
review the experiments that have been performed in the 
different regimes. A direct study of the Bogoliubov dis- 



persion was performed in (Kohnle et al. 2011 ) by means 



of an angle-resolved four-wave mixing technique as pro- 
posed in ( jWouters and Carusotto 2009) . In this way. 



the transition from the parabolic single-particle disper- 
sion at low densities to the sonic dispersion at high densi- 
ties was observed: under the pulsed pump configuration 
used in this experiment the polariton Bogoliubov disper- 
sion recovers in fact the standard one (55) for equilib- 
rium condensates, the only difference being the global 
—ilLp/2 decay rate. In contrast to previous studies of 
the Bogoliubov dispersion with photoluminescence teen- 
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hiques ( Utsunomiya et al. 



2008 1 , the improved sensitiv- 
owed to detect a signal also 



ity of four-wave mixing al 
from the negative Bogoliubov branch. 



V. PARAMETRIC OSCILLATION REGIME 

Historically, the first configuration that was used to ex- 
perimentally study the spontaneous onset of macroscopic 
coherence in the polariton gas and the collective dynam- 
ics of coherent polariton fluids in planar microcavity sys- 
tems was in the so-called optical parametric oscillation 
(OPO) regime. In this section, we shall review the main 
properties of the OPO condensate from the point of view 
of hydrodynamics. An early review of this physics can 



be found in ( Ciuti et al. 2003 1 



A schematic sketch of the pumping scheme is shown in 
the upper panel of Fig[7j a coherent pump is shined on 
the microcavity at a finite incidence angle such that the 
wavevector kp lies in the vicinity of the inflection point 
of the lower polariton dispersion. Because of collisions, 
polaritons are then scattered into a pair of other modes 
at different wavevectors ks and k^. In nonlinear optical 
terms, the modes at k^ and k^ represent the signal and 
zdter modes of an optical parametric amplifier/oscillator, 
whereas the kp mode is the pump mode: the p+p — > s + i 
parametric scattering is mediated by the strong x^'^^ op- 
tical nonlinearity resulting from polariton-polariton in- 
teractions as shown in the lower panel. Differently from 
most other parametric devices, the peculiar shape of the 
lower polariton dispersion and the choice of the pump 
wavevector kp in the vicinity of the inflection point al- 
lows for the parametric process to occur in a remarkable 
triply-resonant way with conservation of both total en- 
ergy and momentum. 



A. Parametric amplifier and parametric luminescence 

For low pump intensities, the parametric scattering 
process takes place in a spontaneous way and produces an 
incoherent luminescence into the signal and idler modes, 
characterized by a broad spectral distribution in both 
energy and wavevector (i.e. angle, see middle panels of 
Figj?]). If a additional seed beam is incident on the mi- 
crocavity around the signal (idler) wavevectors, this gets 
coherently amplified ( ,Savvidis et al.\ 2000ap (see FigjsJ 
and in turn generates a four-wave mixed beam at the 
idler (signal) wavevector (Ciuti et al. 2000). Following 



the first observations, many authors have experimentally 
unveiled different aspects of this sort of parametric am- 



plifier physics ( 


Erland et al. 




2001 




Messin et al. 2001 


Saba et al. 


2001 


), with a special attention to its coherent 



and quantum aspects ( 


Huynh et al. 


2003 


Kundermann 


et al. 


2003 


Savasta et al. 


2005 
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FIG. 7 (Color online) Top panel: sketch of a polariton- 
polariton scattering process conserving both total energy and 
momentum. The first experimental observation of such pro- 
cess was reported by Sav vidis et al.\ |2000a| Middle: contour 
plots of angle- and frequency-resolved polariton parametric 
scattering in the spontaneous (top) and stimulated (bottom) 
regime. Figure from |Savvidis et al.\ |2000b| Bottom panel: 
in the coherent regime, such scattering process can be seen 



as a polariton four-wave-mixing process (Ciuti et al. 2000 



20011 due to the nonlinearity of the excitonic component, in 
particular via the Coulomb interaction. Figures are courtesy 
of P. G. Savvidis. 



B. Optical parametric oscillator 

As usual in bosonic systems, scattering processes can 
be stimulated by an existing population in the final 
states, here the signal and idler modes. When the pump 
intensity is strong enough for the stimulated scattering 
rate to overcome losses, a new stationary regime with a 
macroscopic occupation of single signal and idler modes 
is reached and a pair of coherent signal and idler beams 
are emitted with a narrow distribution in the energy and 



in-plane wavevector (Baumberg et al. 2000 Stevenson 



et al. 2000). Other general aspects of the parametric os- 
cillation in microcavity polariton systems were reported 
in ( Butte et al. 2003 Houdre et al. 2000b Savvidis et al. 
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FIG. 8 Left: amplification spectra due to stimulated polari- 
ton scattering (Savvidis et al. 2000a I. If the pump is cir- 



cularly polarized, due to the scattering selection rules, the 
amplification (responsible for the very sharp peak in the red 
curve) occurs for co-circularly polarized seeding, while no ef- 
fect is visible for counter-circular polarization. Right: power- 
dependance of the emission under continuous wave pump ex- 



citation. Figures from ( Baumberg et al. 2000 1 



2000b Tartakovskii et al. 2002 ) 



From a fundamental point of view, the paramet- 
ric oscillation in spatially extended geometries such as 
planar microcavities is an interesting example of non- 

The sym- 



1975). 



equilibrium phase transition (Haken 
metry that is spontaneously broken at the critical point 
is the U{1) symmetry corresponding to the simultane- 
ous and opposite rotation of the global signal and idler 
phases by the arbitrary angle tp, 



C. The signal/idler condensate 

Early theoretical modelling of the polariton parametric 
effects in planar microcavities was based on a three-mode 
approximation and consisted in a set of coupled nonlinear 
differential equations for the time-evolution of the lower 
polariton field C-number amplitudes in the pum p kp, sig- 



nal ks, and idler k^ modes ( Ciuti et al. 2000 Gippius 



et al. 2004 Whittaker 2001 



2005a). This corresponds 



to making a three-mode ansatz for the lower polariton 
field: 



^LP{K t) = S{t) (5k,k, e-'^^' + P{t) ,5k,k, e-^" 

+ I{t) 5, 



ik,k,e-*"'*, (65) 



where the signal/idler frequencies and wave vectors are 
related by = 2ujp — ujs and k^ = 2kp — kg. Below 
the pump oscillation threshold, the pump-only solution 
with 5* = / = is dynamically stable. Above threshold, 
the pump-only solution becomes dynamically unstable 
(as shown in the B' panel of Fig|6| and a new steady- 
state is reached with finite signal and idler amplitudes 
that spontaneously break the U{1) phase symmetry. De- 
pending on the specific pump parameters, a variety of 
bistable and reentrant behaviors as a function of pump 
intensity were anticipated; in particular, the paramet- 
ric threshold can be a second-order one with a continu- 
ous growth of the signal/idler intensity past the critical 
point, or a fi rst-order one with a sudden jump to the 
OPO r egime (Whittaker 2005a| [Wouters and Carusotto 



2007c|) . Experimental investigations of this rich nonlinear 
lies h£ 
2009| 



dynamics ha ve appeared in ( [Baas et al.\ |2004b| [Demenev 

et al. 



1. Spatial patterns 



I ^ le""^ 



(63) 



A main drawback of the three-mode approximation 



even t hough the mean-field equation s are (approxi- 
mately (Wouters and Carusotto 2007b[ )) invariant under 
the U{1) symmetry, a specific value of the signal and idler 
phases is chosen at every instance of the experiment. The 
strong analogies between the parametric oscillation op- 
eration and a Bose-Einstein condensation transition are 
apparent in the Penrose-Onsager criterion 



lim 



(Et(r)E,(r'))^0 



(64) 



based on the ansatz ( 65 ) is that it is limited to spatially 
homogeneous geometries and does not provide a criterion 
to determine the specific modes k^^^ into which OPO op- 
eration will take place. Naively, one could consider the 
kg i modes for which parametric gain is the strongest and 
then perform a local density approximation. However, 
kg j themselves strongly depend on the pump intensity 
and the theory has to be solved self-consistently. Some 
insight can be obtained by applying techniques borrowed 
from pattern formation theories in nonlinear dynamical 
systems (Cross and Hohenberg 1993). 



Alternatively, the spatial pattern can obtained by 
numerically integrating the non-equilibrium Gross- 

Pitaevskii equation. This was first done by [Whittakerl 

that was used in the experiment ( Baas et al. 2006) and in |2005bt A discussion along these hues recently appeared 



the numerical simulations of ( Carusotto and Ciuti 2005 ) 
to assess the appearance of long-range order in the signal 
(or equivalently the idler) emission. In both theory and 
experiments, the signal field Es(r) was isolated from the 
pump and idler ones by an angular filtering in k-space. 



in (Marchetti and Szymaiiska 2011 ) for a flat-top pump 
spot geometry (see Fig|9j). As a function of pump in- 
tensity, the parametric oscillation operation first appears 
in the center of the spot and then extends to the whole 
pumped region (top-right panels). More details on the 
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FIG. 9 Figures from ( [Marchetti and Szymariska| [20TT| ) . Top 
panels: calculated evolution of the signal, idler and pump 
state intensities as a function of the renormalized pump inten- 
sity (left) and space profiles of the filtered signal at different 
values of the pump intensity in the presence of some photonic 
disorder potential (right). Middle paneh photonic density 
spectra in momentum-frequency space. Bottom panel: full 
emission (first panel) and filtered emission of signal (second 
panel), pump (third panel) and idler (fourth panel). The ar- 
rows in the second panel indicate the current pattern that is 
present in the steady state of the signal field. 



spectral shape of the polariton wavefunction are shown 
in the middle panels of FigjOj as we are dealing with 
a steady-state, the pump, signal and idler emissions are 
monochromatic at uip_s,i- On the other hand, the finite 
spatial size of the emission is visible in their significant 
broadening in k space around their central wavevectors 
The weaker satellite emissions that are visible at 
frequencies ujs2 = Zws — Wp, uJi2 = 2uii—uip are due to sec- 
ondary parametric scattering processes; experimentally, 



they were first observed in (Savvidis et al. 2001 Tar- 



takovskii et al. 



2002D . 

The real-space shape of the polariton wavefunction in 
the parametric regime is illustrated in the lower left panel 
of the same figure. As it was mentioned in (Carusotto 



and Ciuti 2005), interference of the signal, idler and 
pump beams is responsible for moving fringes in the po- 
lariton density pattern: from this point of view, the OPO 
phase transition can be interpreted as the spontaneous 
formation of a stripe pattern in the intensity profile of 



optics literature ( 


Lugiato et al. 


1997 Oppo et al. 


1994 


Staliunas and Morcillo 


2003 


Vaupel et al. 


1999 


I. 



In the finite pump spot geometry under considera- 
tion here, the interplay of nonlinearities, parametric gain 
and losses is responsible for the appearance of a com- 
plex structure of macroscopic currents also in the the 
filtered signal emission (second panel from the left). 
In (Marchetti et al. 2010), it was shown that the behav- 



ior can be even richer when the pump beam is spatially 
narrow and the configuration with spatially smooth sig- 
nal and idler can become dynamically unstable towards 
the spontaneous appearance of quantized vortices. Once 
the vortex has reached its equilibrium position, the sys- 
tem remains stably in the new steady-state with a macro- 
scopic current flowing around the cloud. 



As originally pointed out in (Whittaker 2007), vortices 
can also be forced into the system by seeding the signal 
(or idler) modes with a temporally short pulse with a 
Laguerre-Gauss spot profile. Such a triggered paramet- 
ric device is often indicated by the acronym TOPO. For 
suitably chosen pump conditions, the orbital angular mo- 
mentum imprinted by the seed can remain in the signal 
polaritons for macroscopically long times until the pump 
beam is switched off. As we shall see better in Sec VILE 
the presence of several metastable states with non-trivial 
macrocopic current patterns is one characteristic signa- 
ture of superfluid behavior. 



2. Coherence of the signal/idler condensate 

The Gross-Pitaevskii description reviewed in the pre- 
vious section is able to capture a good deal of the mean- 
field physics of the signal/idler emission well above the 
threshold. As usual in statistical mechanics, fluctuations 
play a crucial role in the critical region in the neighbor- 
hood of the transition point, which calls for more sophis- 
ticated theoretical approaches to sort out the analogies 
and differences between the parametric critical point and 
the standard Bose-Einstein condensation phase transi- 
tion in equilibrium statistical mechanics. 

A numerical study of the coherence properties of the 
signal emission across the OPO threshold was performed 



in (Carusotto and Ciuti 



Carlo technique of Sec I 



2005 1 using the Wigner Monte 
.C to include quantum fluctua- 



tions. Below threshold, it results from numerical calcula- 
tions that the signal emission has a incoherent, thermal 
nature with a short-ranged first-order coherence and a 
significant spatial bunching on a distance scale set by 
the coherence length. As the threshold is approached, 
the coherence length becomes of the order of the system 
size. Related features in the correlation functions were 



addressed in the nonlinear optical literature in (Gatti 



and Lugiato 1995 Zambrini et al. 20001 and often go 



under the name of "quantum images" (Lugiato et al. 
1997). A recent work on the universal critical fluctua- 



tions around the optical parametric threshold appeared 
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in (Drummond and Dechoum 2005). Above threshold, 



first-order coherence extends to the whole system and 
intensity fluctuations are almost completely suppressed 
(^(2) I These numerical observations are in agree- 
ment with the usual picture of second order phase tran- 
sitions and with the experimental measurement of coher- 
ence by |Baas et cd] |2006| 

As the polariton system is naturally a two-dimensional 
one, one may expect that fluctuations have a major 
impact on the long-distance coherence of the signal 
emission. At equilibrium, the Bose-Einstein condensa- 
tion phase transition in two-dimensions is replaced by a 
Berezinski-Kosterlitz-Thouless (BKT) transition to a su- 



perfluid but still non Bose-condensed state (Minnhagen 



1987) and true long-range order of a Bose condensate 
is only observable at strictly zero temperature. So far, 
the spatial size of the systems considered in numerical 
calculations of two-dimensional OPO systems was not 
large enough to numerically solve this issue. Indirect 



evidence based on analytical calculations of ( Szymanska 



et al. 2006) for a different but related model suggests 
that also at non-equilibrium the transition is of the BKT 
type. 

Quasi-condensation behavior is instead clearly visible 
in calculations for a one-dimensional photonic wire ge- 
ometry: even well above threshold, coherence has an 
exponential decay in spa ce ( Carusotto and Ciuti[ [2005 



Wouters and Carusotto 20061. At the threshold, the 



coherence length shows a clear increase and grows to 
macrocopic values, still it does not diverge to infinity. 
On the other hand, the suppression of intensity fluctua- 
tions is almost unaffected by the reduced dimensionality. 
Pioneering experimental observation of suppressed coher- 
ence in one-dimensional parametric oscillators appeared 



in ( Cerda-Mendez et al. 2010) and then in fSpano et al. 
2012). An experimental study of the BKT-like power- 



law decrease of coherence in a two-dimensional polari- 
ton gas under incoherent pumping has recently appeared 

Following the theoretical pro- 



in (R oumpos et al 
posal in (Small et 



20121 



al. 2011), pioneering experimental 



studies of the BKT transition of classical nonlinear waves 



have recently been presented (Situ and Fleischer 2012) 



using an optical configuration similar to the kinetic BEC 



experiment of (Sun et al. 2012) 



Experimental investigations of the temporal coherence 



of the signal emission were reported in ( Krizhanovskii 



et al. 2006 ). On one hand, the interplay of the OPO op- 
eration with spatial inhomogeneities due to sample disor- 
der (mainly of photonic origin) and the pump beam pro- 
file was responsible for the simultaneous presence of sev- 
eral coherent signal beams at slightly different frequen- 
cies with different spatial distributions. On the other 
hand, the temporal coherence of each frequency compo- 
nent showed a coherence time in the few 100 ps range, 
almost a couple of orders of magnitude longer than the 
polariton lifetime. Even l onger coherence tim es in the 
ns range was reported in (Spano et al. 2011) using a 
high-quality sample where the OPO emission enjoys a re- 
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FIG. 10 (Color online) Figure from ( Wouters and Carusotto| 
|2007b[ ). Top panel: imaginary part of the energy of the ele- 
mentary excitations without (left) and with (right) symmetry- 
breaking probe. Bottom panel: same but for the real parts. 
The letter G indicates the Goldstone mode. 



markably good spatial uniformity and spatial coherence: 
the experimentally observed trends suggest the combined 
effect of number fluctuations and inter-particle interac- 
tions in the condensate as the main source of decoherence 



according to the single mode model of (Whittaker and 
Easthaml 120091). 



3. Elementary excitations: the Goldstone mode 

As usual in the statistical mechanics of phase transi- 
tions, the spontaneously broken continuous U{1) symme- 
try is responsible for the appearance of a soft Goldstone 
mode in the dispersion of the elementary excitations of 
the system; in the long- wavelength limit fc — > 0, its fre- 
quency dispersion wg(^) tends to zero in both the real 
and imaginary parts. Physically, the Goldstone mode 
corresponds to a spatially slow twist of the order param- 
eter, which feels a vanishingly small restoring force in the 
long- wavelength limit. 

Among the most celebrated examples of Goldstone 
modes in condensed matter physics we can mention the 
zero-sound mode of superfluid Helium 4 or dilute Bose- 
Einstein condensates and the magnon excitations in fer- 
romagnets: zero-sound is related to the spontaneous 
breaking of the U{\) gauge symmetry of the quantum 
Bose field below the Bose-Einstein condensation temper 
ature ( Forster[ 



1995[ Huang 1987 Pines and Nozieres 



1998), while the magnon branch is related to the spon- 
taneous breaking of the rotational symmetry of the 
magnetic moment orientation below the Curie temper- 



ature (Landau and Lifshitz 1980). 



Calculations of the elementary excitation spectrum 
on top of an OPO coherent state were performed 



in (|Wouters and Carusotto 
a t 



2006 



2007b I starting from 
iree-mode model of OPO operation. Small spatial 
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fluctuations of the mode amplitudes are included in the 
form 



S 



5 + u,e'('^'=^-"*) 



-i{Ak x-Auj t) 



(66) 



P P + Mpe*('^'=^""*) +w;e-*('^'=^~^"*) (67) 



/ ^ I + u^e'^'^'''^""*^ +<e-*('^''^-'^"*\ (68) 

and the dispersion Aa;(A/c) of the elementary excita- 
tions is obtained by plugging this ansatz into the full 
wave equation (37). A pioneering experimental inves- 
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tigation of these excitations was reported in (Savvidis 
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et al. 2001|), including the remarkable observation of a 
significant off-branch luminescence from the negative Bo- 
goliubov modes. 

An example of theoretical dispersion taken 



from (Wouters and Carusotto 2007b[ ) is shown in 
the left panels of Fig|10[ note in particular the Gold- 
stone mode (indicated by the " G" label) whose frequency 
tends to zero (in both real and imaginary parts) as 
Afc — )■ 0. Its functional dependence on Afc is however 



FIG. 11 (Color online) Experimental observations of polari- 
ton condensation under incoherent pumping with a large spot 
of diameter 25 ^m. The different panels show the emission 
pattern in the {k,E) plane for different values of the pump 
intensity. The condensation threshold is situated at a pump 
intensity value between the second and the third panel. Fig- 
ure from (IRichardl 12004 1 . 



very different from the usual one (|55|) of equilibrium 
Bose systems which starts as uj{k) 2 
with a singularity at k — 0. 

Here, no singularity appears in the dispersion rela- 
tion ujci^k) of the Goldstone mode around Afc = 0. 
The continuous and non-vanishing slope of the real part 
Re[w(Afc)] cx Afc is due to the flow of the pump polari- 
tons which are injected with a finite wave vector kp and 
are then able to drag the elementary excitations. On the 
other hand, the imaginary part oi u)Q[l^k) has the low- 
Afc form Im[w(Afc)] k, —a (Afc)^ with a positive a > 0. 
As a result, the Goldstone mode of a planar OPO consists 
of a spatially slowly varying twist of the signal and idler 
phases: a localized perturbation of the signal/idler phase 
will not propagate as a sound wave, but rather relax back 
to the equilibrium state while being dragged by the pump 
polariton flow. Almost simultaneously, a similar diffusive 
behavior of the Goldstone mode was obtained in a com- 
pletely different model of polariton condensation under 
a incoherent pump in ( Szymahska et al. 2006 ) ; a review 
of these results will be given in Sec |VI[ 

The soft nature of the Goldstone mode is destroyed 
in the presence of external flelds that fix the direction 
of the order parameter. While this can appear as quite 
artificial for Bose-Einstein condensation of material par- 
ticles ( [Gunton and Buckingham 1968 Huang, 1987), it 
naturally happens in ferromagnets in the presence of an 
external magnetic field that imposes the direction of mag- 
netization. As it was discussed in (Wouters and Caru- 
pottOj ^20075] ), the signal/idler phases of an OPO can be 
fixed by an additional, weak laser field that resonantly 
drives the signal (or the idler) mode. The effect of this 
additional weak laser field on the elementary excitation 
spectrum is illustrated in the right panels of FigjlOj now 
the imaginary part of the dispersion does not tend to 
zero for Afc — >■ 0, which signals the presence of a restor- 
ing force than tends to bring the signal/idler phase back 
to their steady-state determined. 



|fc| for smafl fc's VI. CONDENSATES UNDER INCOHERENT PUMPING 



In Sec. III.D we have introduced the generalized Gross- 
Pitaevskii equation (49) describing the mean-field dy- 



namics of a driven-dissipative polariton condensate un- 
der an incoherent pumping. In this section we shall re- 
view its application to specific configurations of actual 
experimental interest and we shall compare the theoreti- 
cal predictions with observations. Even though our focus 
will be concentrated on microcavity polariton systems for 
which the analogy with BEC of material particles is the 
clearest, most of the features are expected a much wider 
application range and can be observed to any spatially 
extended laser oscillator, for instance vertical cavity sur- 
face emitting lasers, macroscopic laser cavities at large 
Fresnel number, or even laser oscillation in disordered 



media, the so-called random lasing (see e.g. Tureci et al. 
20081). 



A. The condensate shape 

The first quantity to consider is the shape of the con- 
densate in both real space and in momentum space: an 
unexpected feature of the first, pioneering experiments 
of polariton condensation under an incoherent pumping 
was in fact the strong dependence of the momentum dis- 
tribution of the coherent polaritons on the size of the 
pump spot (Richard et al. 2005a|b[ ). For a large pump 
spot (cr « 20 /im), the momentum distribution was wide 
but centered at k = as shown in the two rightmost 
panels FigfTT] On the other hand, for a small pump spot 
(cr w 3 /xm) the momentum distribution was centered on 
a narrow ring in k-space, as shown in Fig |12[ d). 
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FIG. 12 Experimental observations ol polariton condensa- 
tion under incoherent pumping with a small spot of diameter 
3 fim. The first and second rows are obtained below and above 
the stimulation threshold, respectively. Left column: far-field 
emission. The middle column: emission pattern in the {k, E) 
plane. Right panels: intensity profile along the azimuth angle 
on the 6 = 17° ring indicated as a dashed line in (a,d). Figure 
from (iRichard et a/.| |2005b|). 



the second erm on the right-hand side accounts for the in- 
teractions between the condensate polaritons with an in- 
teraction constant g^p already discussed in (31 1. In con- 



trast, an explicit expression for the last term Awi.cs(-P) 
describing the (generally repulsive) interactions with the 
reservoir polaritons requires a microscopic modeling that 
goes beyond the scope of the present review. In recent ex- 



periments (Wertz et al. 20101, it appears that this latter 



term is the dominating one in the pumped region. The 
fact that in such a model a stable condensate can appear 
in any k mode is of course an artifact of the assumption of 
a spatially homogeneous pump P(r) = P and disappears 
in more sophisticated calculations including in the gen- 
eralized GPE ( 49 ) the s patial profile of the (u navoidably 
finite-size) pump spot (Wouters et al. 2008]): the drift 



of polaritons outside the pumped region reduces the ef- 
fective amplification of k 7^ moving condensates and 
provides a criterion to select the mode into which con- 
densation does occur. The next two subsections describe 
the steady-state condensate shape for simple geometries. 



1. Spatially homogeneous system 



2. Condensation under a large pump spot 



A physical explanation of these observations can be 
obtained in terms of the generalized GPE that in the 
present case of an incoherent pumping has the complex 
Ginzburg-Landau equation form (49). As a first step. 



let us consider the spatially homogeneous case with a 
spatially constant pumping P(r) — P and without any 



trapping potential (Wouters and Carusotto 



2007a) 



For low pumping P < Pc = JrJlp, the only dy- 
namically stable mean-field solution is = 0. For 
a stronger pumping P > P^ the stimulated scattering 
overcomes losses and the zero mean-field solution be- 
comes dynamically unstable. Depending on the initial 
fluctuation that seeds the condensation process, a stable 
condensate can appear in any momentum mode. 



^Lpiv,t)^^>lpe 



.-ikr ^—iujt 



with a condensate density 

nLP = \^Ip? = 



IR 

R 



- 1 



(69) 



(70) 



Given the U{1) symmetry of the GPE equation under 
global rotations of the condensate phase, this latter is 
arbitrarily chosen at every instance of the experiment 
according to the same spontaneous symmetry break- 
ing mechanism as in the equilibrium BEG phase transi- 
tion (Gunton and Buckingham 1968). A completely dif- 



ferent approach to non-equilibrium condensation based 



on a Bose- Fermi model was proposed in ( Szymahska 



et al. 2006 ) , giving the same qualitative features for the 



phase transition. 

The condensate frequency is determined by the pump- 
ing intensity via the equation of state 



a;Lp(k) + gLP |*lpP + At^rcs(^') 



(71) 



The condensate shape in the case of a large pump spot 
with a smooth spatial profile was theoretically studied 



within the local density approximation in (jWouters et al. 
2008|). While the local density nLp(r) 



Whplj 



univocally determined by the local pumping intensity via 
(70), more care has to be paid to the phase: in contrast 



to the standard GPE ( 50 ), breaking of time reversal sym- 



metry by the pumping and loss terms of ( 49 ) allows for 



the steady-state wavefunction to have a spatially vary- 
ing phase, corresponding to macroscopic currents flowing 
across the condensate cloud. 

As the frequency u; of a condensate is constant 
throughout the cloud, the spatial variation of the pump- 
ing intensity -P(r) and of the condensate density nip(r) 
has to be compensated in ( 71 1 by a non-trivial spatial de- 



pendence of the local wavevector k(r) = V[Arg[\E'Lp(r)]. 
The geometric details of the resulting polariton current 
are easily understood in the simplest case of cylindrically 
symmetric pump with a monotonically decreasing pump 



intensity along the radial direction, see Fig 13^a-d); an 



example of the complex current pattern that appears in 
a polariton condensate under a OPO pumping for was 
shown on the bottom row in Fig[9j 

Under the reasonable assumption that rotational in- 
variance is not spontaneously broken, the wavevector 
at the center of the cloud k(r = 0) = and k(r) 
is everywhere in the outward radial direction, k(r) — 
k{r)er. From the explicit forms of ( [7T| ) and ( [70| , the 
flow wavevector fc(r) is found to be an increasing func- 
tion of the radial coordinate r and gets to its maximum 
value /cmax at the external edge of the condensate cloud 
where the condensate density vanishes (?T.(7'max) = 0). 
Within the local density approximation, the radius r„iax 
of the condensate cloud coincides with the point where 
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FIG. 13 (Color online) Numerical results of generalized GPE 
simulations in the absence of disorder for respectively a large 
CTp = 20 fim (a-d) and a small = 2 fim (e-h) circular excita- 
tion pump spot. Panels (a,e) give the (fc, E) emission pattern, 
(b,f) the polariton distribution in momentum space, (c,g) the 
polariton distribution in real space and (d,h) the local wave 
vector fcc(r). The wave vector kco of a free polariton at Uc is 
indicated by the dotted lines in (a,b,e,f); the dashed line in 
(g) is the analytical approximations to the density tail, and 
the dashed line in (d,h) are the LDA predictions to the local 
wave vector. All quantities in real (momentum) space de- 
pend only on the radial coordinate r — |rj (fc = |k|). Figure 
from JWouters et al.\ 120081) . 



the local pump intensity equals the critical pump value, 
P{fmaic) — Pc- Most remarkably, the maximum momen- 
tum fciiiax corresponds to the value of wavevector at which 
the free polariton dispersion equals the condensate fre- 
quency WLp(^max) = 

Even though no local measurement of the condensate 
phase has yet directly confirmed this prediction for the 
spatial dependence of k(r) in wide condensate clouds. 



several other experimental observations (Richard et al. 



2005a) are in agreement with this theory: in particu- 
lar the emission pattern in the (fc, E) plane consists of 
horizontal segments at a constant w, whose extremes co- 
incide with the intersection with the free polariton dis- 
persion. Another indirect evidence of the presence of 
a non-trivial current pattern in the condensate is vis- 
ible in the non-centrosymmetric shape of the momen- 
tum distribution shown in the right panel of FigjTT} the 
Fourier transform /(k) of any function /(r) with a spa- 
tially constant phase has in fact a centro-symmetric mod- 
ulus|/(-k)P = |/(k)p. 



3. Condensation under a small pump spot 

The situation is different when the size a of the pump 
spot is smaller than the ballistic propagation distance 



and polaritons are able to travel away from the pump 
spot before decaying. In this geometry, a quantitative 
calculation of the frequency oj and of the peak density 
requires including the kinetic energy terms of ( 49 ) and 



has to be performed numerically. Nonetheless, impor- 
tant qualitative information on the condensate shape in 
both real and momentum space can be obtained under 
the reasonable assumption that rotational invariance is 
not broken. In particular, an analytical solution for 
the condensate wavefunction far from the pump spot 
can be obtained for non-interacting polaritons in the 



form (Wouters et al. 2008) 



A 



-r/2£ 



(72) 



with A a normalization factor, fc^ax is the above-defined 
maximum wavevector and £ = v-axax/lLP is the spatial 
decay rate of the ballistically flowing polaritons at speed 
Vmax = Vkw(fcmax)- The additional \l\fr dependence is 
a consequence of particle conservation for a radial flow in 
two dimensions. The non-interacting polariton assump- 
tion is generally valid as soon as the density n(r) is low 
enough for gLpn{r) < ujLpiKj,^^). 

Under the condition £ ^ cr, we expect that the to- 
tal number of polaritons ballistically expanding at /cmax 
is larger than the number of polaritons under the pump 
spot and that the momentum distribution will be approx- 
imately given by the Fourier transform of ( 72 1 . As a re- 



sult, the condensate shape in momentum space consists 
of a ring in k space of radius fcmax and width propor- 
tional to This analytic guess was numerically ver- 
ifled by [Wouters et al.\ |2008[ the result for the steady- 
state condensate wavefunction \I'Lp(r) is illustrated in 
Figjl3je-h) and provides a theoretical interpretation to 
the observations by [Richard et al.\ |2005b| reported in 
FigjT2] This picture of ballistic motion of condensed 
polaritons away from the pumped region has been con- 



firmed by later experiments (Christmann et al. 2012 



Wertz et al. 2010 1 measuring the local wavevector of the 
expanding polaritons at different positions. 

For the sake of completeness, it is worth mentioning 
that an alternative explanation of the observed ring in- 
volves fragmentation of the condensate into several sub- 
condensates at different k's. The conceptual difference 
between these two states is formally illustrated as 
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v/2 V 1 



N 



1 



|vac) ^ 7^(4) {aiyi |vac). 

\ 2 >■ 

(73) 

at equilibrium, energy argu- 



For condensates 
ments 



Leggett 




2001 


Nozieres 




1995 



favor 



the 

coherent superposition state in the LHS of ( |73[ ), but 
these arguments no longer apply in the driven-dissipative 
polariton case. Mutual coherence of the emission in the 



±k directions was experimentally measured in (Richard 



2004[ Richard et al. 2005b ) with a Billet interferometer, 
which confirmed that one is indeed dealing with a single, 
non-fragmented condensate with a spatially complex 
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wavefunction ^Lp{r). A direct observation of the 
mutual coherence of polaritons bahistically moving in 
opposite directions away from the small pump spot was 



reported in (Wertz et al. 2010) 



4. Effect of disorder 

A complete account of the experimental findings re- 
quires including in the model the effect of the disorder 
that is unavoidably present in the microcavity samples, 
mostly of photonic origin. The main effect of disorder is 
to introduce a strong modulation of the condensate den- 
sity, the highest density being concentrated around the 
bottom of the disorder potential wells. It is then natural 
to wonder about the mutu al coherence of the diff erent 
regions: in the experiment (Kasprzak et al. 2006), the 



spatial coherence of the different regions was assessed by 
measuring the contrast of fringes when light from differ- 
ent spatial positions is made to interfere. This observa- 
tion was a crucial step in establishing the Bose-condcnsed 
nature of the polariton system in contrast to the Bose- 
glass phase described in (Malpuech et al. 20071). The 



temporal dynamics of the onset of long-range coherence 
was experimentally investigated in ( jNardin et al. , 2009[ ) . 
Remarkably, the problem of the spatial and temporal co- 
herence of the laser emission in disordered systems is one 
of the key issues of present-day research in random lasers, 
i.e. devices that exploit disorder to facilitate laser oscil- 
lation (|Wiersma||2008|). 



This physics was further investigated in (Baas et al. 
[2008) : when the potential barrier separating the high 
density spots is too large for tunneling to lock the phase, 
a synchronized state with a single condensate is re- 
placed by a desynchronized state with multiple indepen- 
dent condensates. A theoretical model of phase locking 
phenomena in non-equilibrium polariton condensation in 



multiple- well geometries was put forward in (Eastham 



2008 Wouters, 2008). 



While all these features have a counterpart in equilib- 
rium BEC in ultracold atomic gases or liquid Helium, 
novel features stem from the interplay of disorder with 
the possibility of macroscopic currents in the steady-state 
wavefunction. Numerical examples of such patterns for 



the case of a wide pump spot are presented in (Wouters 



et al. 2008 1 : while the short-distance profile of the cur- 



rent pattern is strongly affected by the disorder, its global 
shape with a radially outgoing current is robust against 
disorder. 

From a closer look at the real-space experimental im- 
ages of the non-trivial condensate phase pattern, it was 
possible (Lagoudakis et al. 20081 to ascertain the pres- 



ence of quantized vortices in the steady-state conden- 
sate. On one hand, the number and the position of 
the vortices strongly depends on the pumping condi- 
tions, e.g. the position of the pump spot with respect 
to the disorder potential and its intensity. On the other 
hand, their fixed position for given pumping conditions 





FIG. 14 (Color online) Figure from f Lagoudakis et a/. [[2008 1 



(a) Interferogram with a vortex: the fork-like dislocation can 
be seen in the red circle, (b) same interferogram, but this time 
the vortex is overlapped with a different region of the conden- 
sate and a different fringe orientation is selected. The vortex 
appears to be pinned in the same position, (c) Real-space 
phase profile calculated from the interferogram in (a). The 
red circle encloses the vortex (same real-space area as in a,b). 
(d) Phase as a function of the azimuthal angle for a range of 
different radii as shown in the inset of (d) (magnification of 
(c)). 



pinpoints a deterministic mechanism for vortex nucle- 
ation and excludes an interpretation in terms of ther- 
mal vortices as in the Berezinskii-Kosterlitz-Thouless 
transition to a non-superfluid state in finite-temperature 

In the 



two-dimensional Bose gases (Minnhagen 



1987) 



present experiment, vortices spontaneously appear at the 
interface of spatial regions with macroscopic currents in 
different directions. With a careful choice of the pump 
beam profile, vortex-antivortex pairs were experimen- 



tally observed (Roumpos et al. 2011). The spontaneous 
appearance of half-vortices in spinor condensates via a 



related mechanism was observed in (Lagoudakis et al. 
2009). The observation of topological defects when sev- 



eral atomic condensates are merged together was re- 



ported in ( Scherer et al. 2007 ) 



A striking self-oscillation effect was observed in the 



time-resolved data reported in (Lagoudakis et al. 2010) 



for a polaritonic Josephson junction. With a careful 
choice of the position on the sample allowed to identify a 
spatial region where the spatial dependence of the disor- 
der potential exhibits a sort of double-well configuration 
and the condensate density is concentrated at the bottom 
of the two traps. Above a certain value of the pump in- 
tensity, time-resolved measurements of the relative pop- 
ulation and of the relative phase showed an oscillatory 
behavior that bears analogies with the ac Josephson ef- 
fect. 
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5. Condensation in a trap potential 

Differently from tlie experiments discussed in the pre- 
vious section wliere tlie spatial size of tlie condensate 
was mostly determined by tlie size of the pump spot, the 
work by |Balili et al.\ |2007| used a mechanical stress to 
confine polaritons in a harmonic potential as discussed in 
Sec jlLDl This strategy took inspiration from experiments 
with atomic condensates and the experimental findings 
were similar: the condensate appears at the bottom of the 
harmonic potential, and its size is determined by an in- 
terplay of harmonic trapping and repulsive interactions. 
A similar observation was reported in the photon BEC 
experiments by [Klaers et al.[ |2010a| 

Some care has to be paid if one wishes to use the gener- 
alized GPE ( 49 1 to model a harmonically trapped polari- 



et al. 2012 ) suggests that condensation is occurring into a 



ton condensate under a spatially homogeneous incoher- 
ent pumping: as shown in ( Keeling and Berloff 2008 ) , a 



naive application of the theory would in fact predict that 
the rotationally symmetric steady state is dynamically 
unstable towards the exponential growth of high angu- 
lar momentum modes at the surface of the condensate 
and a new stationary state is eventually reached in the 
form of a vortex lattice. The discrepancy of this theoreti- 
cal prediction from the experimental observations can be 
traced back to the finite pump laser spot used in (Balili 



et al. 2007 1 . It is worth mentioning that a regular lattice 



of vortices was instead observed in (Scheuer and Oren- 



stem 



1999) using a VCSEL system electrically driven 



time-dependent polariton wavepacket that bounces back 
and forth within the trap potential and gets amplified 
each time it hits the pumped region. 



7. IVlultiply connected geometries 

An interesting physics can be observed when the ge- 
ometry of the condensate is multiply connected. In the 
simplest case of a ring, the pump intensity is maximum 
along a circle and vanishes at the center. As a result, 
several metastable steady states exist for the condensate 
wavefunction. 



*Lp(r) = f{r)e 



ihe 



(74) 



labeled by the value of the quantized angular momen- 
tum L (ILeggettl 119991 IMueller et aLl 119981). Differently 



well above threshold. 



from simply-connected finite-spot geometries where the 
drift of polaritons outside the pump spot suppresses the 
effective amplification rate and makes a moving conden- 
sate to be quickly replaced by a condensate at rest, fi- 
nite super-current states in multiply-connected geome- 
tries can live for macroscopically long times: unwinding 
the condensate phase requires in fact creating a node in 
the condensate wavefunction, which is a quite unlikely 
event. The strict conceptual link between metastabil- 
ity of super-currents and superfluidity effects will be the 
subject of Sec |VII.E[ 



6. IVIore complex trap potentials 



B. Stability and the elementary excitation spectrum 



Investigations of condensation in complex geometries 
started from the case of a double well potential: early 
work used a pair of neighboring polariton traps cre- 
ated by disorder. This allowed to observe the tran- 
sition between independent condensates to a synchro- 
nized regime where the condensate phases are locked 



by tunneling across the separating barrier ( |Baas et al. 
[2008). Theoretical m odels accounting for the observed 



After having determined the shape of the condensate, 
the next logical step is to investigate the spectrum of el- 
ementary excitations on top of the condensate. For the 
case of an incoherent pumping, this problem was first 



attacked in ( Szymanska et al. 2006 1 using Keldysh dia- 



gram techniques applied to a non-equilibrium Bose- Fermi 
model of condensation. Here we shall follow a different 



route initiated in (Wouters and Carusotto 2007a) and 



featured appeared in (Eastham 2008 Wouters 2008). based on the non-equilibrium GPE (49), which leads to 



More recently, a similar physics was investigated in an 
artificially fabricated double micropillar system (Galbiati' 
et al. 2011). A time-dependent regime showing sponta- 



neous Josephson oscillations between the two wells were 



qualitatively similar predictions in a perhaps more trans- 
parent way. The main result of these works is the diffu- 
sive behavior of the Bogoliubov dispersion at low k's, 
which recovers the almost simultaneous work (Wouters 



observed in (Lagoudakis et al. 2010). Theoretical stud 



ies of Josephson physics in coherently driven polariton 
systems appeared in ( Sarchi et al. 2008 ) and first exper- 



imental investigations were reported in (Bloch 



Polariton condensation into several trapped states of 



2012) 



complex wavefunction has been observed in ( Wertz et al. 



2010 ) using a combination of the repulsive effect of the 



pump and of the geometrical end of the polariton wire. 
Later measurement of the emission coherence in (Tosi 



and Carusotto 2006 2007b ) for the planar optical para- 
metric oscillator reviewed in Sec lV.C.31 The fact that the 
same behavior is found in several microscopically very 
different models suggests that it must be a generic fea- 
ture of non-equilibrium phase transitions. 

In the simplest case of a spatially homogeneous con- 
densate at rest \l/ip(r, = -Julp e~*"*, the Bogoliubov 
equations describing the dynamics of small fluctuations 
have the form 
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FIG. 15 (Color online) Dispersion of the Bogoliubov modes 
of an incoherently pumped condensate. The ± branches show 



the prediction 1 77 1 for the equilibrium Bogoliubov modes. 



The R branch indicates the fast relaxation of the reservoir at 
rate 7h according to ( |48[ |. The diffusive behavior of (^Bog(fc) 
is clearly visible at low k's. Figure adapted from ( |Wouters| 
and Carusotto 2007a I. 



with an effective damping rate 
PRn 



{jR + Rnf 



1 ■ 



P/Pc - 1 
P/Pc 



(76) 



Within the usual parabolic band approximation, the 
eigenmodes of ( 75 1 give the dispersion 



'^Bog(k) 



(77) 

for the elementary excitations. An example of such a dis- 
persion is plotted in FigJTSj two regimes are clearly iden- 
tified depending whether the argument of the square root 
in (77) is positive or negative. The transition wavevec- 
tor ko corresponds to the point where the equilibrium 
Bogoliubov dispersion crosses the effective damping F: 
as expected, for smaller loss rates 7, the effective decay 
rate T decreases and the system behavior recovers the 
equilibrium case. 

For low values of the wavevector k < kg, the argument 
of the square root is negative and the real parts of the 
two branches stick to the common value Re[a;Bog(k)] = 



0, while the imaginary parts are split. Stability of the 
condensate is guaranteed by the fact that for all k's one 
has Im[a;Bog(k)] < 0. 

As fc 0, one of the two branches tends to a finite, 
purely imaginary value — 2iF: from an explicit calcula- 
tion of the eigenvector, one sees that it corresponds to 
a spatially homogeneous fluctuation of the condensate 
density and does not affect the condensate phase. When 
approaching the critical point from above P — >■ P+, its 



decay rate F tends to zero. In the laser literature (Mandel 



and Wolf 1995 1, such a mode is involved in the intensity 



dynamics of the device at switch-on and determines the 
temporal correlation function of intensity fluctuations 
Of course, a more complex dynamics than simple relax- 
ation is observed as soon as the gain medium shows a sig- 
nificant frequency dependence of the amplification rate. 

A similar amplitude mode is encountered in many con- 
texts of condensed-matter physics when a continuous 
symmetry is spontaneously broken at a thermal or quan- 
tum phase transition. The simplest case is the super- 
fluid to Mott-insulator transition for bosons in a lattice. 
Close to the critical point on the superfluid side, a second 
branch of excitations appears in addition to the Bogoli- 
ubov sound. Physically, this new Higgs-like mode can be 
understood in terms of an oscillating amplitude of the su- 
perfluid order parameter, that is a local interconversion 
of the condensed and non-condensed fractions (Altman 
and AuerbachI [2002{ |Huber et aL[ |2007| [Menotti and 



Trivedi 20081. A summary of physical systems where 



this physics may be observed (atomic gases, antiferro- 
magnets, charg e-density- wave system s, superconductors) 
is reported in (Podolsky et al. 2011): very recently, ex- 



perimental evidence of the amplitude mode was reported 



for ul tracold atoms in an optical lattice (Bissbort et al. 
2011). A related physics was investigated for a Dicke 



model of polariton condensation in (Brierley et al. 2011 1. 



The other mode which tends to in both real and 
imaginary parts as A: — >■ is the Goldstone mode cor- 
responding to the U{1) spontaneous symmetry breaking 
associated to the BEC transition: given the underlying 
U{1) phase rotation symmetry of the problem, the phase 



* In standard lasers, the decay rate F of intensity fluctuations re- 
main finite even at threshold because of the relatively large fluc- 
tuations in a single mode l |Mand cl and Wolf, 1995|. Even though 
such a prediction r(P = Pc) = was obtained within mean-field 
theory, this result is expected to be robust against fluctuations 
provided a spatially extended system is considered. The lin- 
earized fluctuation theory of jSzymanska et aL\\2006^ appears to 
conflrm this expectation. 
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of the condensate is randomly chosen at the phase transi- 
tion and a global rotation of the phase can not experience 
any restoring force. This interpretation is confirmed by a 
calculation of the corresponding k = eigenmode of ( 75 ) 



that indeed describes a pure global phase rotation. Most 
remarkably, the so-called Schawlow-Townes linewidth of 



a standard laser (Schawlow and Townes 1958 Walls and 



Milburn 2006 ) can be seen as arising from the excitation 
of this soft mode by spontaneous emission events. 

For finite < k < kg, the amplitude and phase modes 
are mixed by the kinetic energy term in ( 75 ) but maintain 



while the imaginary part of the dispersion remains com- 
pletely unaffected. This prediction that moving conden- 
sates in ring geometries are dynamically stable up to ar- 
bitrarily high values of the flow speed vq is at odd with 
intuitive expectations that condensation into fast moving 
states should be unfavored. A more refined treatment 
including an energy-dependent amplification mechanism 
was developed in (Wouters and Carusotto 2010): as ex- 



a diffusive nature with a decay rate tending to zero in the 
long wavelength limit. A perturbation imprinted at the 
initial time i = on the condensate phase (p{r) will not 
propagate through the system, rather it will diffusively 
expand in space according to a heat equation 



dt r 



VV(r,i), 



(78) 



Cs = a/ gLpTiLp/fn being defined as the usual sound ve- 
locity at equilibrium. 

For k > ko, the imaginary parts collapse to a single 
value — r/2, while the real part starts from zero propor- 
tionally to VA; — ko- At larger values of fc, the dispersion 
recovers the standard Bogoliubov dispersion of equilib- 
rium systems (black dashed line in Fig 15). As usual 
in Bogoliubov theory, the off-diagonal terms in the ma- 
trix in (75) are responsible for a mixing of the ^5'Lp(r) 
and <5^^p(r) field, i.e. particle and holes operators in 
a quantum framework. As a result, both the positive 
and the negative branches should be observable in lumi- 



pected, it turns out that the region of stability extends 
up to a maximum speed value vq < Ustab- For higher 
speeds vq > Vgtah^ the instability of the condensate is due 
to the birth of another condensate around k = which 
eventually replaces the one at kp. 



VII. SUPERFLUID HYDRODYNAMICS OF THE 
PHOTON FLUID 

Supcrfluid behaviors were first observed in 1938 al- 
most simultaneously by Allen and Miesner in Cam- 

|1938 ) and Kapitsa in 



nescence experiment (Byrnes et al. 


2011 


Keeling et al. 


2005 


Marchetti et al. 


2007 


). A re 


ated prediction was 



put forward for atom laser beams extracted from atomic 



condensates in (Japha et al. 1999). 



bridge (Allen and Misener 

Moscow ( Kapitsa |1938 ) in liquid Helium-4 cooled at tem- 
peratures below the so-called A point at Ta = 2.17 K. 
The first experimental signature was a sudden drop of 
mechanical viscosity when the Helium fluid was flowing 
along a narrow channel. Since then, a number of other 
fascinating features of superfluid Helium-4 have been un- 
veiled (Leggettj 1999 Pitaevskii |1992 ), from metastahil- 
ity of supercurrents (with int eresting applications as a 
gyroscope for rotation sensing ( Sato and Packard] 2012 )), 
to novel excitation branches with out-of-phase oscil a- 
tions of the normal a nd superfluid compon ents (the so- 
called second sound ( Peshkov 1944 1946)), to the re- 
duced moment of inertia ( Andronikashvili and Mama! 
'adze/1966) and the Hess-Fairbank effects (Hess and Fair- 



Even though such negative Bogoliubov branches (also 
called ghost branches) have been observed with polari- 



ibank, ,1967| ) in rotating bucke t exper iments, to the foun- 
tain effect (Allen and Jones 1938), to the frictionless 



motion of impurities through the fluid when the speed 



ton systems in several other pum ping regimes ( |Kohnle| is below some critical value (iPhillips and McClintock 



et aL\ |2011[ [Savvidis et a/.[|2001[ ), no trace ofjt was vis- |1974| |Rayfleld[ |1966|). Even though our conceptual un 



ible in the recent experimental study of the Bogoliubov 
dispersion in an incoherently pumped condensate in a 
planar geometry (Utsunomiya et al. 2008): Most likely. 



this is due to the fact that the emission from this nega- 
tive branch is easily masked by the much stronger back- 
ground of the condensate emission. On the other hand, 
a spectral feature that may be attributed to the negative 
branch was observed in a pillar geometry by |Lagoudakis[ 

mm 

Extension of this theory to the supercurrent states in 
multiply-connected geometries described in Se c|VI.A.7| 
requires more attention. A naive application of ( 75 ) to a 
moving condensate at speed Vq = Kko/m would predict 
(within the parabolic approximation for the lower polari- 
ton branch) a mere Doppler shift of the real part of the 
Bogoliubov dispersion, 



derstanding of superfluidity in liquid Helium is now based 
on London's intuition that superfluidity originates from 
Bose-Einstein condensation of the constituent quantum 
particles and on Tisza and Landau's two fluid hydro- 
dynamics, the experimental evidence of the presence of 
a Bose-Einstein condensate in superfluid Helium is still 
quite elusive and based on indirect neutron scattering 



data ( 


Azuah et al. 


1997 


Glyde et al. 


2000 


2011 Sokol 


1995 


Sosnick et al. 


1990 


). 



The interest in the physics of superfluids was suddently 



revived in 1995 when Cornell and Wieman at JILA (An- 



,(vo) 



(k) = WBog(k) + k • vo, 



(79) 



derson et al. 1995 ) and soon after Ketterle at MIT (Davis 
et al. 1995 ) observed Bose-Einstein condensation in a gas 
of laser- and then evaporatively cooled atoms stored in 
a magnetic trap at nanokelvin temperatures. In contrast 
to liquid Helium where the condensate fraction is limited 
to around 10%, atomic gases have a low enough density 
to be in the dilute regime where the condensate fraction 
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is large and clearly visible in the experimental data. 

In the immediate wake of the first evidence of conden- 
sation, researchers started investigating the superfluid 
properties of the ultra-cold atomic clouds. Thanks to 
the dilute nature of the atomic gas and the flexibility of 
atomic manipulation and diagnostic techniques, a direct 
comparison with microscopic calculation was possible, 
which allowed to get a deeper insight into the basic mech- 
anisms underlying, e.g. the critical velocity for friction- 
less motion of weak impurities across the cloud (Raman 



et at. 1999), the nucleation of quantized vortices at the 



surface of a large defects (Neely et al. 2010), the effect 



of a reduced moment of inertia onto the collective scis- 



sor mode of an anisotropic cloud (Marago et al. 2000), 



the metastability of supercurrents in multiply-connected 
geometries ( |Ryu et al. [2007 ). In spite of these excit- 
ing advances, many subtle features of superfluidity of ul- 
tracold atomic clouds are still waiting for experimental 
confirmation, in particular for what concerns the con- 
ceptual distinction between the superfluid and the con- 
densed fractions, in particular in low-dimensional geome- 
tries (Carusotto 2010). 



In this section we shall give an overview of the re- 
cent theoretical and experimental developments in the 
superfluid hydrodynamics of light. Along the same lines 
of classical studies of superfluidity in condensed-matter 



systems (Leggett 



1999| ), different aspects of superfluid 



behavior will be addressed: in contrast to standard equi- 
librium systems, the different aspects of superfluidity can 
appear separately in driven-dissipative systems, which 
gives rise to a richer variety of effects. 



A. The Landau criterion 

A most celebrated textbook presentation of superflu- 
idity (Pitaevskii and Stringari 2004) is based on the 



so-called Landau criterion for superfluidity, which de- 
termines the maximum speed at which a weak impurity 
can travel across a fluid without experiencing any friction 
force. This criterion was originally proposed by Landau 
to explain superfluidity in liquid Helium and led to the 
prediction of the roton minimum in the elementary exci- 
tation spectrum of the fluid. 

In terms of the dispersion uj{k) of the elementary ex- 
citations in a generic fluid at rest, the Landau critical 
velocity has the form 



Vc — mm 
k 



uj{k) 



(80) 



In its simplest formulation, the idea underlying the Lan- 
dau criterion can be summarized as follows. We con- 
sider a fluid uniformly moving at speed v: as seen from 
the laboratory reference frame, the elementary excitation 
spectrum in the fluid has the Galilean-transformed form 



For slow flow speeds v < Vc^ & finite positive amount of 
free energy is always required to create excitations in the 
fluid at any value of the momentum k. For faster flows 
V > Vc, there exist values of k for which w'(k) < 0, so it 
is thermodynamically favorable to create a large number 
of excitations in the moving fluid. As a result, superflu- 
idity is broken and an effective friction force appears. Of 
course, the validity of the Landau criterion is restricted 
to impurity potentials (due, e.g., to the roughness of the 
container walls) that are weak enough to modify the fluid 
density only in a perturbative way. As we shall see in 
the following of this section, non-perturbative processes 
in the presence of stronger impurities can create bound 
complexes of many elementary excitations (e.g. vortices) 
even at v < Vc', experimentally, this appears to be the 
case in most experiments with macroscopic Helium sam- 
ples. 



1. Bogoliubov formulation of the Landau criterion for weak 
impurities: the equilibrium case 

Landau's thermodynamical argument can be made 
complete and quantitative in the framework of the Bo- 
goliubov theory of dilute gases in the limiting case of a 
weak impurity at rest in a uniformly moving homoge- 
neous condensate ( Astrakharchik and Pitaevskii^ ,2004) 



Carusotto et al. 2006). A unified discussion of the 



wave emission by uniformly moving sources in the dif- 
ferent contexts of superfluid hydrodynamics, Cerenkov 
effect in classical electromagnetism, and classical hydro- 
dynamics of surface waves can be found in (Carusotto 
and Rousseaux 2012 1. 

compact notation 54>{r, t) = 
the weak perturbation created 



Introducing the 

{S<j>{r,t)M*{r,t)r, 
by the impurity in the moving condensate at speed 
V = hko/m can be straightforwardly obtained within 
Bogoliubov theory by the linear equation 



d ^ - - 

ih—Stjj = £v • Si' + Fd, 
at 



(82) 



where the source Fd is proportional to the defect po- 
tential Vd and the unperturbed condensate wavefunction 
(/)o(r) = exp(jko • r), 



'/'o(x) 



(83) 



The Bogoliubov operator for the moving fluid has the 
form: 



(84) 



w'(k) =tj(k) 



^1) 



where the condensate motion is clearly visible as a phase 
factor in the off-diagonal matrix elements. As usual, the 
dispersion of the Bogoliubov modes in the moving fluid 
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FIG. 16 Dispersion of Bogoliubov modes in the flowing con- 
densate (left column), k-space locus E of the resonantly ex- 
cited Bogoliubov modes (center column) and spatial profile of 
the density wake around an impurity at rest (right column). 
The different rows correspond to different values of the (su- 
personic) flow speed v^jc — 1.1, 1.5, 3. The red crosses in the 
left and central columns indicate the condensate wavevector 
ko. Figure adapted from (Carusotto et al. 20061. 



is obtained from the eigenvalues of £v, 



4;)jk) = v.(k-ko 



± 



;i(k-ko)2 /fi(k-ko 



2m 



2m 



2gpo) : (85) 



the finite velocity of the condensate is apparent in the 
Doppler shift term (k — ko) • v. 

The steady state in the presence of the defect potential 
Vd is then obtained from the motion equation (82) as: 



(86) 



where the infinitesimal imaginary term is required to en- 
sure causality of the solution. As it can be seen from 
the resonant denominator, the time-independent defect 
potential is able to effectively excite only the k 7^ kg 
Bogoliubov modes whose energy is Wg^g(k) — 0. The 
shape of the k-space locus E of such modes is shown in 



the central column of Fig 16 for a few values of the flow 
speed. According to the Landau criterion (80 1, two main 
regimes can be identified depending on whether this is 
slower or faster than the speed of sound Cg in the BEC. 

In the subsonic regime < Cs, the intersection lo- 
cus E is empty. As a consequence, no Bogoliubov mode 
can be resonantly excited by the defect and the fluid 
is able to flow in a frictionless way around the defect 
without suffering any dissipation. Nevertheless, the non- 
resonant excitation of the Bogoliubov modes by the im- 
purity is responsible for a sizable density modulation in 



the vicinity of the impurity, that quickly decays to zero 
in space with an exponential law and which results in a 



sizable renormalization of the mass of the impurity ( As- 
[trakharchik and Pitaevskii 2004). Of course, all these 
conclusions are based on a mean-field description of the 
condensate that neglects quantum fluctuations: more so- 
phisticated Bethe ansatz calculations for a strongly in- 
teracting one-dimensional Bose gas ( Astrakharchik and 



Pitaevskii 2004) have anticipated the appearance of a 
finite drag force also at sub-sonic speed. A calculation 
including higher order terms of the Bogoliubov theory led 
the authors of ( Roberts and Pomeau 2005 ) to a similar 



claim for a three-dimensional condensate; this conclusion 
however appears at odd with classical literature on liquid 



Helium (lordanskii and Pitaevskii 1978). 



In the supersonic regime vq > c, the intersection locus 
E is not empty and consists of the closed curve E shown 
in the central column of FigjlB] Resonant excitation of 
these Bogoliubov modes makes the kinetic energy of the 
flow to be partially dissipated into Bogoliubov phonons 
and is associated to a sizable drag force exerted by the 
fluid onto the defect. Some characteristic features can be 
identified in the corresponding density modulation pat- 
tern: 

(i) The emission of low-wavevector Bogoliubov modes 
in the sonic region is responsible for a Mach cone lo- 
cated downstream of the defect with the usual aperture 
9 such that sin 6 — Cs/vq. This emission mechanism is 
closely analogous to the Cerenkov emission of electromag- 
netic waves by charged particles in super-luminal motion 
across a dielectric, v > Vpi^ — Cq/u (Jelley, 19581. 

(ii) A series of curved precursors located upstream 
and laterally to the defect that result from the single- 
particle interference between the incident plane wave 
and the spherical scattered wave. An analytical calcu- 



Gladush et al. 


2007 


); a 


discussed in (Carusotto 



and Rousseaux 


2012 


); their one-dimensional analog 


studied in ( 


Leboeuf and PavlofI 2001 


Pavloff 


2002 



2. Experiments with liquid Helium and atomic condensates 

As originally anticipated by Landau, the critical ve- 
locity (80) is determined in liquid Helium by the ro- 
ton branch of excitations. Experimental studies of the 
drag force as a function of the velocity were reported 
in ([Phillips and McChntock} |1974[ |Rayfield[ |1966[) using 



ions as microscopic impurities traveling across a macro- 
scopic Helium sample. Interpretation of the experimen- 
tal data was made difficult by the number of competing 
phenomena that may be simultaneously taking place, e.g. 
generation of charged vortex rings. 

In recent years, a much more detailed experimental 
access to the microscopic physics underlying the Landau 
criterion has been provided by experiments with ultra- 
cold atomic clouds. Some remarkable results are illus- 
trated in Fig 17 In the upper-left and middle ones, the 
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FIG. 17 Upper left and middle panels: experimental stud- 
ies of weak impurities travelling across dilute atomic conden- 
sates. Left: relative density of scattered atoms as a function 
of the impurity speed vq/cs, showing a marked threshold at 
vo — Cg. Figure from |Chikkatur et al] |2000| Middle: atomic 
momentum distribution after elastic scattering between the 
condensate (lower black dot) and super -sonically moving im- 
purity atoms (upper spot); figure from (Ketterle and Inouye 
2001 1. Upper right panel: real-space density profile of a BEG 



hitting an obstacle at supersonic velocity Wq/cs = 13 in the 
left direction. Figure courtesy of P. Engels and E. A. Gor- 
nell |GornelH |2004[ Lower panels: sequences of images showing 
the dynamics of vortices nucleated by the repulsive potential 
of a blue-detuned laser scanned across a pancake-shaped con- 
densate in the rightwards direction. Figure taken from |Neely| 



impurity consists of atoms of the same species crossing 
the condensate cloud at finite speed after a Bragg scat- 
tering process. The onset of drag as soon as the speed 
of the impurity exceeds the speed of sound is clearly vis- 
ible in the right panel. The ring-shaped momentum dis- 
tribution of scattered atoms shown in the middle panel 
provides information on the Bogoliubov modes that are 
excited by the moving impurity. The right panel shovifs 
the real-space wake created in an expanding condensate 
by the strong repulsive potential of a blue detuned laser 
beam: while the Cerenkov cone in the downstream re- 
gion is strongly modified by a trivial shadow effect, the 
curved precursors upstream of the defect closely resemble 



theoretical prediction in Fig 16 Experimental studies of 
this physics in polariton fiuids were reported in Sec.fAmo 
et al.\ 2009a I and are reviewed in Sec VILE. 2 and Figsj2 
andlioi 



3. Beyond the Landau criterion: strong defects 

In most experiments involving macroscopic bodies 
moving in liquid Helium or strong potentials in atomic 
condensates, the critical velocity for the onset of drag 
turns out to be much lower than the prediction of the 
Landau criterion (80). A possible explanation for this 



of the strong modification of the velocity pattern around 
a large and strong defect: as a result of the spatial com- 
pression of the stream lines, the local velocity on the 
surface of the defect can largely exceed the flow speed at 
infinity. As soon as the critical velocity (80) is locally 



exceeded, quantized vortices start being nucleated at the 
surface of the object and give rise to a drag force. Ac - 
cording to the first calculations in ( [Frisch et aL 19921, 
this happens when the asymptotic speed far from the 
defect exceeds a critical value of a/2/11 = 0.43 times 
the spee d of sound. Subsequent, m o re refined calcula- 

|200r 



tions in (Berloff and Roberts 2000 



Rica 2001) 



predict a slightly different value 0.37 for the critical Mach 
number, in agreement with numerical calculations. Fur- 
ther numerical work for penetrable objects was reported 
in (Winiecki et al. 2000). This dissipation mechanism 



based on vortex nucleation was experimentally confirmed 
in ( Neely et al. 2010|) using a strong blue-detuned beam 
crossing a pancake-shaped atomic condensate: examples 
of experimental snapshots of the density pattern are re- 
produced in the bottom panel of Fig |17| Experimental 
investigations of this physics with polariton condensates 
were reported in ([Nardin et al.\ |2011a[ [Sanvitto et al. 



2011 ) and are reviewed in Sec VIII and Fig 27 



As originally predicted in (El et al. 20061, at fiow 



speeds larger than the speed of sound, vortex nucleation 
is replaced by the appearance of one or several pairs of 
oblique solitons in the wake of the large defect. Stabi- 
lization of the dark soliton against the snake instabil- 
ity and the subsequent decay into a string of vortices is 
prevented by the large flow speed as explained in ( Kam- 



chatnov and Korneev 2011 Kamchatnov and Pitaevskii 



2008 ) . While no evidence of this physics has yet been ob- 



tained yet in liquid Helium nor in atomic clouds, exper- 
imental studies of hydrodynamic generation of o blique 
dark soliton s in polariton flui ds were reported in (Amo 
et al. 2011) and soon after in (jG rosso et aL[[20TT |. This 
physics will be reviewed in Sec | VIII and illustrated in 
FiglU 



B. Coherently pumped polariton condensate flowing 
against a weak defect: superfluidity effects 

The first experiments unambiguously showing super- 
fluidity effects in luminous systems (lAmo et al. 2009a) 



were performed under a coherent pump with a quasi- 
plane- wave spatial profile, along the lines of the proposal 

In this subsection we 



20041. 



behavior was proposed in (Frisch et al. 1992) in terms 



in (Carusotto and Ciuti 
shall review the theoretical concepts underlying the ex- 
periment and then we shall illustrate the experimental 
data. 



1. Theory 

The value of the in-plane pump wavevcctor kjnc and 
of the pump frequency Wjnc are chosen in the parabolic 
region at the bottom of the lower polariton branch, so 
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potential is a time-independent one, the emission of Bo- 
goliubov modes is concentrated on those modes for which 
the real part of the frequency Re[cij(k)] is zero; the effect 
of the finite lifetime of the Bogoliubov excitation on the 
broadening of the emission as well as on the drag force 
was discussed in (Cancellieri et al. 2010). The peculiar 



k (|Lm-i) 



FIG. 18 Left panels: Bogoliubov dispersion of elementary ex- 
citations in a coherently pumped driven-dissipative polariton 
condensate flowing in the rightward direction. Right panels: 
corresponding far-field emission patterns in the presence of 
a weak defect at rest. The pumping regimes considered in 
the different rows arc explained in the text. Figure taken 
from ( |Carusotto and Ciuti, ,2004^ . 



to avoid the parametric instabilities discussed in SecfV] 
and have, in the absence of defects, a spatially uniform 
polariton condensate steadily flowing along the microcav- 
ity plane. However, as we have discussed in Sec IV.B.l 



shape of the k-space locus of excited Bogoliubov modes 
directly reflects on the shape of the k-space density pat- 
tern, which is concentrated on the resonant modes: a few 
examples of such patterns are shown in the right column 
of Fig{T8] 

Panels (a,b) on the top row of Fig 18 illustrate the 
case of low-intensity resonant pump, where polariton- 
polariton interactions are negligible. In this case, the 
locus of resonant modes has a circular shape and the k- 
space emission pattern is characterized by the so-called 



resonant Rayleigh scattering ring ( Freixanet et al 



Houdre et al. 2000a Langbein and Hvam 2002 1 . An 



1999 



the non-equilibrium nature of the polariton system elimi- 
nates the one-to-one relation between the field oscillation 
frequency cjjnc and the density n^p and is responsibile for 
the larger variety of Bogoliubov dispersions discussed in 
Sec irvXl 

The Bogoliubov matrix for a coherently pumped po- 
lariton con densate flowing at speed v = Hkinc/m was first 
derived in dCarusotto and Ciuti 2004 Ciuti and Caru- 



sotto 



analogous feature was illustrated in the left panel of 
FigjT7|in the case of ultracold atomic gas. 

Panels (c,d) on the second row correspond to the case 
of an interacting polariton gas in supersonic motion: the 
pump frequency and intensity is adjusted at the end- 
point C of the bistability loop of Fig[6] so that the low-k 
dispersion has a sonic shape with a sound speed lower 
than the speed of flow. This results into a strong defor- 
mation of the k-space resonant Rayleigh scattering ring 
into an eight-shaped pattern. 

Panels (e,f) illustrate the case of a subsonic motion, 
where the polariton gas behaves as a superfluid: no Bo- 
goliubov mode can be any longer resonanly excited and 
the resonant Rayleigh scattering ring disappears. The 
weak emission that is still visible in panel (f) corresponds 
to non-resonant scattering events favored by the finite 
linewidth of the polariton modes. 

All these features are in close parallel to what is pre- 
dicted for dilute atomic gases. Novel behaviors that stem 
from the non-equilibrium nature of the polariton fiuid are 



2005^ Its form is closely similar to the matrix illustrated in the two lowest rows of Fig 18 Panels (g,h) 



(61) for the polariton condensate at rest except for the 
phase factors e^^'''"" in the off-diagonal terms as in the 



Bogoliubov matrix ( 84 ) for an equilibrium condensate in 
motion. Its eigenvalues give the dispersion of elementary 
excitations. The most relevant examples are shown in 
the left panels of Fig 18 the effect of the finite kjnc is to 
Doppler shift the branches by (k — kjnc) ■ v. As usual, 
the frequency and wavevector of the Bogoliubov modes 
is defined with respect to the oscillation frequency Wjnc 
and the wavevector kjnc of the unperturbed condensate. 

Superfluidity of the polariton fiuid in the spirit of the 
Landau criterion can be investigated by looking at the re- 
sponse of the flowing polariton fluid to a local and time- 
independent impurity potential modeling, e.g., a fabrica- 
tion defect in the microcavity sample or an optically gen- 
erated potential. The perturbation induced in the fluid 
by the weak defect is straightforwardly obtained using 
the formalism review in Sec |VII.A.l| with the appropri- 
ate form (87) of the Bogoliubov operator: as the defect 



correspond to the case where the pump intensity is on 
the lower part of the bistability loop (i.e. in the vicinity 
of point A on Figj6]) and the Bogoliubov branches stick 
within a finite region in k space. As a consequence of the 
reduced linewidth of the Bogoliubov modes, the k-space 
emission is reinforced on a short segment parallel to the 
ky axis corresponding to the sticking region. 

Panels (i,j) correspond to a optical limiter case or to 
the upper branch of the bistability loop where the Bogoli- 
ubov branches are separated in energy. In this case, the 
k locus still has a ring-like shape but no longer intersects 
the pump wavevector. Of course, for larger separation of 
the Bogoliubov branches one may even recover a super- 
fiuid behavior without any resonant scattering channel 
available. 

Examples of real space density patterns for the differ- 
ent cases are shown in Figjl9) The upper-left panel refers 
to the linear regime of non-interacting polaritons whose 
k-space features are shown in Fig{T8la,b). The density 
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FIG. 19 Examples of real-space polariton density patterns 
for coherently pumped polariton condensates flowing in the 
rightward direction against a defect at rest. The defect is 
located at the center of each panel. Upper-left panel is for 
the non-interacting polariton regime of Fig |18[a ,b). Upper- 
right panel is for the superfluid regime of Fig |18[ e,f). Panels 
from (Carusotto and Ciuti 20041. T he lower-left panel is 
for a super-sonic flow regime of Fig 18 c,d). The lower-right 
panel illustrates the Zehra- Cerenkov effect in the vicinity of 
a parametric instability; the corresponding k-space emission 
pattern is shown in the inset. Panels from ( [Ciuti and Caru- 
sotto||2005|). 



pattern consists of a series of parabolic wavefronts ex- 
tending in the upstream direction of the defect and orig- 
inating from the interference of the plane-wave incident 
condensate with the spherical wave of scattered polari- 
tons. Of course, the finite polariton lifetime makes the 
perturbation disappear at large distances from the impu- 
rity. The upper-right panel corresponds to the superfluid 
regime of Fig{l8f^e,f) where no Bogoliubov mode is reso- 
nantly excited and the density perturbation remains lo- 
calized in a very small region in the neighborhood of the 
impurity. The two panels on the lower row correspond 
to a polariton fluid in supersonic motion: the deformed 
shape of the resonant Rayleigh ring shown in Fig|T8|^c,d) 
is responsible for the rich structure of the pattern with 
curved precursors upstream of the defect and a Mach 
cone extending in the downstream direction. 

The lower-right panel of fig{l9] illustrates a novel 
regime of zehra- Cerenkov effect that is only possible far 
from equilibrium when the condensate is close to a para- 
metric instability: the k-space emission is concentrated 
around two symmetric spots where the linewidth of the 
Bogoliubov modes is strongly quenched as a precursor of 
the instability; correspondingly, the real-space image is 
characterized by a horizontal stripe pattern downstream 
of the defect. 
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FIG. 20 Experimental (top) and theoretical (bottom) real- 
space images of a fast moving polariton condensate hitting a 
point defect. The different columns correspond to increasing 
values of the polariton density (from left to right). Figure 



taken from (Amo et al. 2009a| ) 



2. Experiments 



The theoretical predictions reviewed in the previous 



subsection have been experimentally confirmed in (Amo 
et al.\ 2009a[ ). The most significant images are repro- 
duced in Figsj2] and [20j In this experiment, polaritons 
are coherently injected by a laser beam with a wide spot 
incident on the microcavity sample at a finite angle so to 
generate a polariton directed in the downwards direction. 
The static impurity consists of a localized fabrication de- 
fect naturally present in the microcavity structure. 

The first and second rows of Figj2] show the experi- 
mental images for the real and k-space emission patterns. 
The different columns refer to increasing polariton den- 
sities from left to right. The left column correspond to a 
non-interacting polariton fluid: the parabolic wavefronts 
in real space and the resonant Rayleigh scattering ring in 
k-space are clearly visible. The weaker spot that is visi- 
ble at the upmost point on panel (iv) can be attributed 
to a coherent multiple scattering events on cavity dis or- 
der ( |Langbein et aL[ |2002| |Miiller and Delande| |201lD . 

The right column correspond to a superfluid regime 
where the k-space emission is quenched and the dis- 
turbance created by the impurity remains spatially lo- 
calized in its close vicinity. Together with the experi- 



ment (|Amo et al 
in Sec 



2009b ) that will be reviewed later on 
VII. D[ these observations provide the first exper- 



imental evidence of superfluid behaviors in a luminous 
fluid. Theoretical images obtained by solving the non- 
equilibrium Gross-Pitaevskii equation are shown in the 
lower rows. Exception made for the additional micro- 
cavity disorder that was not implemented in the calcula- 
tions, the agreement with experimental data is very good 
and confirms that the main source of broadening of the 
k space images is the finite spatial size of the polariton 
fiuid. 
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Fig 20 illustrates the case of a coherently pumped po- 
lariton fluid moving along the plane at a faster speed as 
a result of a larger incident wavevector kjnc- In this case, 
the available laser power was not sufficient to penetrate 
a superfluid regime. The supersonic flow regime is vis- 
ible in the real-space patterns shown in the central and 
right columns: the precursors in the upstream direction 
are always visible and the aperture of the Cerenkov cone 
increases for increasing polariton density (i.e. decreasing 
Mach number). 

Soon after, a similar experiment was performed 
in ( Amo et al.\ 2010b I using another laser beam to cre- 
ate the defect potential by injecting control polaritons 
with counter-polarized tT_ polarization with a suitably 
designed spatial distribution The all-optical nature of 
this method appears extremely promising in view of gen- 
erating more complex potential geometries where to trap 
polaritons: linear-shaped barriers were explored in the 
same work ( Amo et al. 2010b|), while o ther configura- 



tions were considered in (Sanvitto et al. 2011) 



C. Ballistically moving polariton fluid 

The main drawback of the coherent pumping configu- 
ration is that the condensate phase is locked by the pump 
laser, apparently restricting the range of hydrodynamic 
phenomena that can occur. In particular, topological 
features such as solitons and vortices are forbidden, as 
they involve the spontaneous development of non-trivial 
phase patterns in the fluid. This issue was first pointed 



out in (Bolda et al. 2001), who proposed to use a tem- 
porally pulsed coherent pump to inject polaritons and 
then follow in time their free evolution when hitting a 
defect. Such a pump configuration was adopted in ( Ser- 



mage et al. 2001 1 to study the in-plane propagation of 
polariton bullets in presence of a finite acceleration. 

Another strategy to study hydrodynamical effects 
in fluids of polaritons is inspired to the experiment 
of ("Richard et al.\ |2005b|) that we have reviewed m 



Sec |IV.BT3 the condensate of polaritons is generated in a 
spatially restricted region from which it propagates away 
in a ballistic way before decaying on a length scale set 
by £ = Vh^t /i- As theoretically explained in (Wouters 
2008), the group velocity Wgr is determined by the 



et al. 



repulsive potential felt by polaritons in the pumped re- 
gion, while 7 is the inverse polariton lifetime. A more 
complete study of the density and flow profile of such 
a ballistically propagating polariton condensate can be 
found in (Kamchatnov and Korneev 2011). 



Focussing our attention onto a small enough spatial 
region where the condensate density n^p and local mo- 
mentum k can be taken as approximately constant, the 
Bogoliubov matrix has the form 
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FIG. 21 Theoretical images of the ballistic propagation of a 
polariton condensate against a large defect. Panel (a): spatial 
profile of the pump spot localized in the top half-plane, just 
above a point defect. Panels (b-f) show the different flow 
behaviors: (b) superfluid propagation around the defect; (c) 
turbulent creation of vortex-antivortex pairs (unlike the other 
images, this one is not stationary in time); (d-f) generation 



of a pair of dark oblique solitons. Figure from ( Pigeon et al. 



exception made for the linewidth of the modes, the 
resulting spectrum is then the standard sonic one of equi- 
librium condensates (84). Note how a similar spectrum 



is obtained for condensates generated by a spatially ex- 
tended, but temporally pulsed coherent pump; a n exper- 
iment addressing this last case was reported in (Kohnle 



et al. 2011) 



A direct proof of coherent ballistic polariton propaga- 



tion was provided by the experiment in (Wertz et al. 



2010) by comparing the momentum distribution in the 
different spatial regions with theoretical predictions and 
then by measuring the interference pattern between po- 
laritons located on opposite sides of the pumped region. 
These observations complements the k-space Billet inter- 



ferometer experiment of (Richard 2004 Richard et al. 
12005b'). Intri guing observation of back-scattering sup- 
pression in supersonically flowing condensates have been 



recently reported in (Tanese et al. 2011 ): the experimen- 



tally observed superfluid flow appears at odd with the 
expectation based on the naive Landau criterion used so 
far. In the same work, a novel possible mechanism for 
creating an effectively superfluid state via parametric in- 
stabilities is proposed. 

The potential of ballistic propagation schemes for 
studying the hydrodynamic nucleation of topological ex- 
citations in a condensate was pointed out in (Pigeon 
using a configuration with a coherent pump 



in a half-plane geometry as shown in Fig 21 a). Panels 
(b-f) show snapshots of the polariton density for increas- 
ing values of the Mach number v/cs. For intermediate 
values of the Mach number, turbulent vortex nucleation 
is observed and then a pair of oblique solitons are formed. 
Qualitatively, this phenomenology is reminiscent of the 
theoretical works in the context of liquid Helium ( Frisch 



(87) 



\t-al. 



2006 Kam 



1992) and ultrac c )ld at oms ( El et al._ 

;hatnov and Pitaevskii [2(508). Experiments along these 
lines will be reviewed in Sec lVIIll 
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D. Spontaneous symmetry breaking and non-equilibrium 
polariton superfluidity 

Historically, the first experiment that has addressed 
superfluidity properties of polariton fluids was reported 
by |Amo et al.\ |2009b[ The chosen configuration is very 
different from the one of the experiments discussed so far, 
and raises many intriguing questions about its theoretical 
understanding. 

A spatially wide, continuous- wave pump is incident on 
the cavity at a finite k- vector close to the inflection point 
of the polariton dispersion with an intensity slightly be- 
low the parametric threshold discussed in Sec|V] A sec- 
ond probe pulse excites the cavity with a small wavevec- 
tor, creating a slowly moving polariton bullet that prop- 
agates within the pumped region. Parametric scattering 
processes from the pump strongly enhance the lifetime of 
the signal polariton bullet fBallarini et al. 2009 1 . The re- 
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FIG. 22 Figure taken from (Amo et al. 2009b I. Top panels 



markable experimental observation is that this polariton 
bullet does not spread in space during propagation and 
is not destroyed by the collision against a defect: snap- 
shots of real and momentum space distributions illustrat- 
ing this superfluid-like behavior are shown in Fig{22} 

Since the experimental results were published, several 
mechanisms have been invoked to explain the observa- 
tions. As -to the best of our knowledge- no consensus 
has been yet reached in the community, we limit our- 
selves here to an overview of some most significant fea- 
tures and open questions. An alternative point of view 



real space images of a "bullet" of signal polaritons hitting a 
natural defect of the microcavity. The bullet of signal polari- 
tons has a very long lifetime thanks to the parametric gain 
induced by the continuous-wave pump. Bottom panel: corre- 
sponding k-space space images. 



hitting the defect, see the right panel of Figj22j The 
spatial modulation of the density profile that is visible in 
the two central panels is attributed to the Cerenkov wake 
that is imprinted by the defect onto the pump beam. A 
naive reasoning based on the Landau criterion applied to 
the elementary excitations spectrum of an OPO conden- 



on this experiment can be found in the review ( Marchetti 
and Szymariska 2011). 



A main difficulty of the theoretical description is that 
the experiment is carried out in the vicinity of the 
OPO critical point, where nonlinear effects play a domi- 
nant role in determining the so-called TOPO (triggered 
OPO) spatio-temporal dynamics. With respect to the 
already complex discussion of the OPO equation of state 



sate shown in Fig 10 does not provide an explaination of 
this experimental observation: because of the diffusive 
nature of the Goldstone mode, a number of k-modes are 
in fact available into which polaritons can be scattered 
by the defect. On this basis, one would rather expect the 
signal polariton bullet to be immediately destroyed after 
hitting the defect. 

The unexpected outcome of numerical simulations of 
the Gross-Pitaevskii equation ( 49 1 for polariton conden- 



sate under i ncohe rent pumping in ( [Wouters and Caru- 



sotto 



2008 



in (Wouters and Carusotto 2007c), the temporal and 



2010) has suggested a possible path to rec- 
oncile the picture. A naive application of the Landau cri- 
terion ( 80 ) to the real part of the diffusive spectrum of el- 



spatial evolution of the signal/idler wavepackets has now ementary excitations ([77|) shown in Fig |15| would predict 
to be fully taken into account. 

Slowly decaying and non-spreading polariton 
wavepacket propagation has been numerically stud- 



that the critical speed vanishes and an impurity is able 
to emit phonons in the fluid independently of the value of 
the flow speed. The numerical results shown in the upper 



led in (Szymahska et al. 2010) by solving the full OPO 



wave equations: the polariton bullet consists of spatially 
superposed signal and idler beams; parametric scattering 
from the pump is able to compensate losses and keep 
the spatial shape constant. For a simplified model with 
a spatially infinite pump beam, non-spreading solitonic 
solution have been numerically found and characterized 



in (Egorov et al. 2010 2011). A dedicated experimental 
study of bright polariton solitons in microcavities 
under a cw-pump in a OPO-like configuration recently 
appeared in (Sich et al. 2011). We expect that these 



studies of soliton propagation may provide useful insight 
in the mechanism underlying the polariton bullets in the 



panels of Fig 23 strongly disagree with this expectation 
and indicate that at low speeds the incoherently pumped 
condensate indeed behaves as a superfluid and is almost 
unaffected by the defect. The critical speed for the on- 
set of the usual wake is of the order of the (equilibrium) 
speed of sound = \/gLpnLp/rnLp. These numerical 
observations are confirmed by the clear threshold that is 
visible close to Cg in the velocity dependence of the fric- 
tion force shown in the lower left panel: the weaker the 
effective loss rate P, the sharper the threshold. 

Taking inspiration from classical work in electrody- 
namics of absorbing media (Tait 1972), an analytical 



2009b). 



TOPO configuration of (Amo et al. 

Another, even more challenging feature of this exper- 
iment is the robust shape of the polariton bullet after 



understanding of these results can be obtained by noting 
that in a stationary state both the condensate wavefunc- 
tion and the Bogoliubov modes oscillate at the single 
frequency w, while the imaginary part of the dispersion 
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FIG. 23 Upper panels: Generalized GPE simulation of an in- 
coherently pumped condensate hitting a weak and stationary 
defect at three different speeds v/cs = 1.5,1,0.4 across the 
critical value for superfluidity (from left to right). Lower left 
panel: force exerted on the defect by the moving condensate 
as a function of the speed v. The different (thin black solid, 
blue solid, green dashed, red dotted lines) curves correspond 
to growing values of the non-equilibrium parameter. Lower 
center and right panels: Real and imaginary parts of the com- 
plex wavevector k of the zero-frequency Bogoliubov mode as 
a function of v. Figure taken from ( Wouters and Carusotto 



20101 



(77) has to be reabsorbed into the complex wavevector 
k. Restricting our attention to zero-frequency WBog — 
Bogoliubov modes propagating along the velocity axis, 
the real and imaginary parts of the complex wavevector 
k are plotted in the bottom central and left panels of 
Figj23]as a function of the speed v: the generalized Lan- 
dau critical velocity Vc corresponds to the threshold for 
the appearance of a non- vanishing real part. The weaker 
the effective decay rate T, the closer Vc to the equilibrium 
prediction c^. 

As we have reviewed in Sec|V] polariton condensates in 
the OPO regime share the same spontaneous C/(l) sym- 
metry breaking mechanism and the same diffusive Gold- 
stone mode as the incoherently pumped ones. On this 
basis, we then expect that the mechanism of superfluid- 
ity illustrated in Fig |15| for incoherently pumped conden- 
sates may provide a physical explanation of the intriguing 
superfluidity observations in the OPO regime reported 



in (Amo et al. 2009b). 



E. Metastability of supercurrents 



A striking manifestation of superfluidity is the 



metastability of supercurrents ( Leggett 1999). For sim- 
plicity let's consider a multiply connected geometry in 
the form of a torus: single-valuedness of the conden- 
sate wavefunction ^(r) imposes the so-called Onsager- 
Feynman quantization condition to the supercurrent 



dl 



2Trn 



-N„ 



where Vs is the superfluid flow velocity and the integer 
N.^; is the so-called winding number indicating the num- 
ber of times the phase of the wavefunction winds up in a 
loop around the torus. 

While in classical hydrodynamics any friction process 
is able to continuously slow down the flow around the 



torus, the quantization condition (88) makes the life- 
time of finite 7^ supercurrents extremely long in 
repulsively interacting superfluids with no extra internal 
degree of freedom even in the presence of a significant 
wall roughness. In order for the winding number 
to vary, a node has in fact to appear in the condensate 
wavefunction, which (exception made for the critical re- 
gion just below the superfiuid critical temperature Ta) 
requires surmounting a very high free-energy barrier. As 
a result, a number of different mctastablc supercurrent 
states exist for the condensate, labelled by the winding 
number A^^,. 

This p hysics was firs t investigated in superfiuid liq 



uid He (Reppy 1965 Reppy and Depatie 1964 



nen 



Vi- 



1961 ) and, more recently, in ultra-cold atom con- 



densates (Ramanathan et al. 2011 Ryu et al. 2007): in 



this latter case, it was explicitly shown that the toroidal 
geometry is essential to guarantees metastability of the 
supercurrent; a vortex in a harmonically trapped con- 
densate can in fact escape from the cloud with no energy 
cost (Fetter and Svidzinsky 2001). 
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FIG. 24 Snapshots of the time evolution of the phase (upper 
row) and of the density (bottom row) of an OPO polariton 
condensate triggered by a pulsed probe carrying a m = 1 vor- 
tex. To improve visibility of the vortex dynamics, the contri- 
bution of the unperturbed polariton signal (in the absence of 
the probe pulse) is subtracted from all data. Figure adapted 
from | |Sanvitto efai] |2010[ ) . 



An experiment along these lines was carried o ut in po- 
lariton condensates in a OPO configuration in (Sanvitto 



et al. 2010): a pump beam with a wide spot profile is 
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shone on the microcavity with an intensity sHghtly above 
the threshold for OPO operation. A temporally short 
and spatially narrow probe pulse in a Laguerre-Gauss 
state with a finite orbital angular momentum is used to 
force OPO operation to occur into a rotating state: after 
the extra polaritons injected by probe pulse have dis- 
appeared, a vortex remains imprinted into the conden- 
sate wavefunction in the signal/idler modes and lasts for 
macroscopically long times as shown in Fig 24 Analo- 



This definition of the normal fraction is widely used in nu- 



merical quantum Monte Carlo calculations ( jPoUock and 
Ceperleyl (19871) to assess the transition to a superfluid 



gously to t he case of harmon ically trapped atomic con- 
densates of ( Ryu et al. 2007 ) , the fate of the vortex is to 



state and then quantitatively evaluate the superfluid frac- 
tion: a most interesting feature in view of studies of ex- 
otic superconductor and superfluid states is that it does 
not require a ■priori knowledge of the microscopic nature 
of the order parameter nor of its symmetry. 

To get a more intuitive physical understanding of the 
definition ( 92 1 , one may note that a purely transverse 



eventually drift out of the condensate: the time-scale of 
this decay process is however orders of magnitude longer 
than the radiative lifetime of a signal polariton. 

Further experimental work has investigated the size of 
the vortex core under a continuous-wave probe beam in- 



vector potential A(r) naturally appears in the description 
of a mechanical system in a rotating reference frame at 
angular velocity 17 ( Cohen- Tannoudji 2001-21, 



jecting angular momentum ( Krizhanovskii et al. 2010) 



Arot(i") = r2 X r 



(93) 



and the appearance of vortex-antivortex pairs when the 
spatial size of the probe pulse is much smaller than the 
pumped area (Tosi et al.. 2011). The dynamical stabil- 



ity of multiply- charged vortices was experimentally ad- 



dressed and numerically confirmed in (Sanvitto et al. tial shape proportiona 



the current pattern induced by the vector potential ( 93 ) 
via the Hamiltonian (189 1 then has a rigid-body spa- 
As 



to the normal fraction /„ 



2010|): while in trapped atomic condensates with simply in the rotating bucket gedanken experiment (iLeggett 



connected geometries, these objects are generally dynam- 
ically unstable towards the splitting into several singly- 
charged vortices (Fetter and Svidzinsky 2001), stable 



1999), the superfluid fraction stays at rest in the fixed 



m — 2 vortices could be observed in the polariton case. 
A related theoretical study of metastable vortices and 
supercurrents in polariton condensates under an incoher- 



star reference frame, while the normal fraction is dragged 
into rotation by the vessel. The Andronikashvili exper- 



iment ( Andronikashvili and Mamaladze 1966) showing 



ent p umping scheme appeared in (Wouters and Savona 
20101). 



F. Response to transverse vector potentials 

From a formal point of view, the most precise and 
quantitative definition of superfluid fs and normal /„ 
fractions involves the response of the quantum gas to a 
weak tran sverse vector potential (Hohenberg and Mar- 
tin 1965). The interaction Hamiltonian describing the 



superfluidity of liquid Helium samples was based on an 
implementation of the rotating bucket idea using a tor- 
sional oscillator. 

Recent developments in the manipulation of ultracold 
atomic clouds have demonstrated the possibility of gener- 
ating artiflcial gauge flelds coupl ed to the motion degrees 
of freedom of neutr al atoms (jPalibard et al. 2011 Dum 
and Olshanii 1996 Lin et al. 2009). The measurement 



of the response of the atomic cloud to a su itably tay 
lored artificial gauge fields was proposed in (Carusotto 
and Castml |2011[ jCooper and Hadzibabic[ |2010 John 



coupling of the particles to a vector potential has the 
form 



et al. 2011) as a way to measure the normal and super- 



V^- Jdh i{r) -Air) 



(89) 



j(r) 



0^(r)V<^(r) -h.i 



(90) 



with the current operator defined as usual as 

h 

2im 

For a spatially homogeneous system of density n, the 
susceptibility tensor relating the average current to the 
applied vector field can be written in Fourier space as 



fluid fractions of an ultracold atomic cloud independently 
from Bosc-Einstein condensation. 

A related theoretical proposal for the case of non- 
equilibrium polariton condensates appeared in ( jKeeling 
2011). A many-body calculation using the Schwinger- 



Keldysh diagrammatic approach for non-equilibrium sys- 
tems has shown that the superfluid fraction of the polari- 
ton gas according to the definition ( 92 ) remains finite in 



(j)(q,t^) = x(q,t^) A(q,w). 



(91) 



and the normal fraction /„ is defined as the low- 
momentum, low-frequency limit of the transverse suscep- 
tibility 



/„ = lim lim — XT(q,a;). 



(92) 



spite of the vanishing Landau critical velocity ( 80 ) and 
the diffusive nature of the Goldstone mode. In contrast to 
the atomic and liquid Helium cases, the non-equilibrium 
nature of the polariton condensate makes the normal 
fraction to remain finite even at T = 0. The artificial 
gauge field for polaritons can be generated by applying 
a real magnetic field to the microcavity in a suitable im- 
balanced anti-Helmholtz configuration. Other proposals 
to implement artificial gauge fields for polaritons will be 
reviewed in Sec lIX.El 
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FIG. 25 Figure taken from ( Amo et al.\\2QV\\ . Experimental observation of hydrodynamically generated oblique dark solitons 
in a polariton superfluid (left). The phase interferogram (middle panel) clearly reveals a phase jump across the soliton lines. 
With bigger defects, a quadruplet of dark solitons is observed (right panel). 



VIII. HYDRODYNAMIC FORMATION OF DARK 
SOLITONS AND VORTICES 

As we have reviewed in the introductory section of 
this article, a great deal of the early literature on su- 
perfluid photon hydrodynamics dealt with the physics of 
quantized vortices in the coherent photon fluid. How- 



ever, in spite of a number of theoretical proposals ( Bolda 



file for a large flow speed well above the speed of sound 
Cs are shown in Fig 25 In addition to the curved pre- 
cursors located upstream and laterally to the defect, a 
pair of oblique dark solitons is apparent in the wake of 
the defect. Their solitonic nature is confirmed by the 
large phase jump ~ tt across the dark region that is 
visible in the phase pattern. In contrast to the moving 
dark solitons that appear past a strong obstacle in a one- 



et al. 2001 Staliunas 1993) and experimental devel- dimensional flow (Engels and Atherton 2007 Hakim 



opments using photorefractive oscillators (Vaupel et al. 1997 Wan et al. 2010), the tip of the oblique solitons 



[1996) and nonlinear optic al crystals in a propagating ge- 
ometry ( ,Wan et al. 2008), no complete study of the fun- 
damental hydrodynamic processes of soliton and vortex 
nucleation at the surface of a large defect was reported. 
Even more surprisingly, even though this mechanism was 
predicted to play a crucial role in determining the actual 
critical speed of superfluid liquid Helium (iFrisch et al. 



shown in Fig|25|is pinned at the defect position. 

The aperture angle a can be related to the Mach num- 
ber and to the phase ju mp across the soliton us i ng th e 
theory of dark solitons (Pitaevskii and Stringari 2004). 
This predicts a soliton speed equal to Vso\ = Cs cos(0/2), 
where the phase jump 6 is related to the relative value 



[l992| ) , its first direct observation in a condensed mat 
ter context appeared only very recently using an atomic 



of the minimum density rimin by cos (0/2) 



i/n. In 



condensate in (Neely et al. 2010) 



order to be at rest in the laboratory frame, the soliton 
speed has to be equal to Wgoi — v sin a, where a is the 
aperture angle of the soliton as indicated in the left panel 



In this section we shall review the recent experimen- 
tal investigations of the physics of a superfluid of light 
hitting a large defect. Thanks to their flexibility, these 
experiments using microcavity polaritons have provided 
the first complete characterization of the different behav- 
iors as a function of the fiow parameters. 



A. Oblique dark solitons 

The first experimental investigation of a flowing polarl 
ton condensate hitting a large defect appeared in (Amo 



of Fig 25 Combining these relations, one obtains for the 
aperture a the expression 



sma = — cos - 
V \2 



(94) 



a direct consequence of the cos(6'/2) factor in (94) is that 



the aperture of the oblique soliton pair is always narrower 



than the one of the Mach cone of linear waves ( El et al. 



2006 Gladush et al. 2007) 



Remarkably, the same cos(6'/2) factor (94) predicts 
that solitons are kinematically allowed also for sub- 
sonic speeds v < Cg] however, theoretical work for loss- 



et al.. 2011): with respect to the earlier demonstration less systems(Brand and Reinhardt 2002 Zakharov and 



by the same group of polariton superfluidity ( Amo e< a/. Rubenchik 1973) has shown that such sub-sonic soh 



2009a), a different fabrication defect in the same micro- 



cavity sample was chosen, with a suitably larger spatial 
extension. Another crucial difference was the implemen- 
tation of the proposal by Pigeon et al.[ |2011| to use a 
pump spot spatially restricted to the upper half space 
upstream of the defect. Without this judicious trick, the 
phase of the condensate would be locked to the pump 
laser and no topological excitation could appear in the 
fluid. Experimental images of the density and phase pro- 



tons should be dynamically unstable against the decay 
into vortices via the so-called snake instability (Ander- 
son et al.\ |2001 ) unless the flow is supersonic (Kamchat 



nov and 



itaevskii |2008) . The experiment (Amo et al. 



2011) appears in partial disagreement with this predic- 



tion, as stable solitons are observed even for sub-sonic 
speeds: according to [Kamchatnov and Korneev] |2011[ 
the discrepancy may be explained as a result of the fi- 
nite lifetime of polaritons. The quite different physics 
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FIG. 26 Different regimes of flow in the presence of a large 
fabrication defect in the microcavity. From left to right, the 
different columns correspond to growing values of the Mach 
number vjcs- For each value, the top panels show the real 
space density profile. The middle panels show the interfero- 
grams giving the phase. The bottom panels show the first- 
order coherence function gr^^' of the emission. Figure taken 
from (lAmo et aZ.1 120111). 



of solitons under a pulsed excitation was investigated 



in (Grosso et al. 2011): within a few picoseconds, soli- 
tons decay into vortex streets. Numerical simulations in 
the same work point out the crucial role of cavity disorder 
in triggering this instability mechanism. The rightmost 
panel in Figj25] confirms the prediction by |E1 et al.\ |2006| 
that for very large defect sizes, the single oblique solitons 
should be replaced by a fan structure of many oblique 
solitons. Very recently, the study of oblique solitons was 
extended to the case of half-solitons in spinor conden- 
sates, where the interplay of the flow and the spin-orbit 
coupling effects discussed in Sec |II.A.2| are responsible 
for the appearance of peculiar spin structures around the 



Flayac et al. 


2011 


Hivet et al. 


2012 



B. Turbulent behavior 



The outcome of this same experiment for different val- 



ues of the Mach number v /cg is summarized in Fig 26 
The left panels correspond to a higher polariton density 
and lower v/cs- a superfluid behavior is visible where 
the perturbation induced by the defect remains spatially 
localized (upper row) and the condensate phase keeps a 
plane wave form (central row). The right panel shows 



again the fast flow case discussed in the previous sub- 
section, where a pair of dark lines are visible in the den- 
sity pattern and are associated to a sudden jump in the 
phase pattern. In both these cases, the coherence func- 
tion g^^'^ remains close to unity at all positions, which 
signals that the condensate is fully coherent. 

As it was expected from the theoretical calculations 
bylFrisch et a/.[ [19921 pigeon et a/.l[20lH|Winiecki et al.\ 



|2000 the situation is completely different for the interme- 
diate value of the speed that is illustrated in the central 
column: while the density profile exhibits a significant de- 
crease behind the defect, the phase pattern shows strong 
irregularities in the same region with a number of vortex 
singularities distributed with no apparent order. An in- 
terpretation of these observations in term of a turbulent 
behavior is confirmed by the spatial behavior of the g^^-* 
coherence function: the regions of suppressed coherence 
in the wake of the defect signal the presence of moving 
vortices with a complex nucleation and drift process. 

Before proceeding, it is interesting to note that the 
corresponding regime in a one-dimensional geometry con- 
sists of a time-dependent nucleati on of solitons past the 
obstacle, as originally predicted in (Hakim 19971 and ex- 



perimentally observed with atomic condensates in (En- 
gels and Atherton |2007 ). Some first e vidence of this ef 



feet in a fluid of light was reported in (Wan et al. 2010) 



using a nonlinear crystal in a propagation geometry. 



C. Vortices 

The experiments reviewed in the previous subsection 
were performed under a spatially localized, but tempo- 
rally continuous-wave coherent pump. This pumping 
regime is ideal to study the stationary state of the po- 
lariton fluid, but has a limited access to time-dependent 
phenomena such as the vortex nucleation process: as this 
effect relies on a dynamical instability of the stationary 
flow at the surface of the defect, the exact nucleation time 
is strongly affected by any experimental disturbance or 
system imperfection. To experimentally observe quan- 
tized vortices in the wake of the defect as in the theoret- 
ical image of Figj2l|c) , one then needs either a very high 
temporal resolution of the imaging apparatus or the pos- 
sibility of taking a single-shot image of the fluid. While 
the former issue can be solved using streak-camera de- 
vices with almost ps temporal resolution, the latter one 
still appears very challenging with present-day technol- 
ogy- ^ 

In the experiments of (Nardin et al. |2011b ) and soon 
after of (Sanvitto et al. 2011), this problem was circun- 



vented using a train of identical short pump pulses. The 
spacing of the pulses is chosen much longer than the life- 
time of excitations in the sample, so that each realization 
of the experiment can be considered as independent. In 
particular, the exact nucleation time of the vortices is 
deterministically pinned to the temporal proflle of the 
pump pulse, which allows to repeat the experiment a 
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FIG. 27 Temporal sequence of experimental images of the 
scattering of a leftward moving polariton cloud on a struc- 
tural defect. The three rows show the polariton density (a), 
the fringes of the measured interferogram (b) and the polari- 
ton phase (c). The defect position is indicated by the green 
circles. Vortices start being visible in the third column and are 
indicated by white markers (x for vortex, -|- for anti- vortex) 
on the density plot and are circled in red on the fringes and 
phase plots. Dotted circles indicate short-lived vortices. On 
the later columns, the previous motion of the long-lived vortex 
pair is indicated by white dots on the density plots. Figure 
taken from (Nardin et al. 2011b|. 



huge number of times and accumulate good statistics on 
the polariton density and phase pattern. As suggested 
in (Bolda et al. , 2001 ), short pump pulse durations guar- 



antee that the condensate phase is free to evolve even 
for pump pulses with a plane wave spatial shape. Fig 27 
shows a series of snapshots at different times before and 
after the arrival of the coherent pump pulse (of duration 
3 ps). The polariton flow is from the right to the left 
and hits a structural defect in the microcavity sample: 
the nucleation of pairs of vortices with opposite circula- 
tion is apparent in the snapshots for t = 4.7 ps shown 
in the third column; after being nucleated, the vortices 
are dragged downstream by the polariton flow. Images 
from an analogous experiment carried out using a gas of 
ultracold atoms by Neely et al.\ |2010| are shown in the 
lower panel of Fig. [17{ 



A similar experiment was reported in (Sanvitto et al. 



2011 ) using an optically generated defect potential ( Amo 



et al. 2010b ) instead of a structural defect of the micro- 
cavity. This apparently minor difference was responsible 
for some qualitative differences in the process of vortex 
nucleation, in particular because of the soft nature of 
the optical potential. This same work also reported im- 



plementation of a proposal in (Pigeon et al. 20111 to 
permanently trap and store the quantum vortices that 
are hydrodynamically generated in the wake of a natural 
defect. The most significant results are summarized in 
Fig 28 the plane-wave, continuous-wave coherent pump 
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FIG. 28 Vortex storage in a triangular trap. Real-space 
images (top row) and corresponding interferograms (bottom 
row) for polaritons injected by a continuous-wave laser beam 
and flowing in the rightward direction against a natural defect 
in the microcavity (marked by rectangles in a). The different 
columns correspond to the no-mask case (a,g) and then to 
different positions of the triangular mask, i ndicated by the 
blue dotted triangle. Figure adapted from (Sanvitto et al. 
20TT|. 



is blocked by a triangular-shaped mask that suppresses 
the coherent pumping in the region inside the dashed 
triangle. If the mask is located right behind the defect 
(b-d), the vortices that are hydrodynamically generated 
by the defect remain permanently stored within the un- 
pumped region where the condensate phase is free. Their 
escape is indeed forbidden by the pump beam that fixes 
the phase in the external region. If the mask is too far 
away (e) or is not present (a), no vortex can be created 
at all because the condensate phase in the neighborhood 
of the defect is fixed by the pump laser. As these figures 
are obtained under a continuous-wave pump with a long 
integration time, the only vortices that are visible in the 
figure are the ones whose position is stationary in time. 



IX. STRONGLY CORRELATED PHOTONS 

Most of the physics reviewed so far originates from 
the collective behavior of a large number of interacting 
photons or polaritons coherently sharing the same wave- 
function. As a result, the system can be described in 
terms of a mean-field theory based on a generalized non- 
equilibrium Gross-Pitaevskii equation, where the driven- 
dissipative nature of the photon fluid is accounted for. 
On the other hand, the mean-field theory is inaccurate 
when classical and/or quantum fiuctuations become im- 
portant and the state of the many-particle system is no 
longer described by a single macroscopic wavefunction. A 
first example of situation where fluctuations play a cen- 
tral role around the critical threshold point (see Secjv] 
for the case of the optical parametric oscillation). As 
long as interactions between polaritons are weak enough, 
a theoretical description of the coherence properties of 
the strongly fluctuation polariton field across the OPO 



48 



threshold in (Carusotto and Ciuti 2005) could be ob- 



tained in a semiclassical way in terms of the stochastic 
Gross-Pitaevskii equation within the truncated Wigner 
representation reviewed in Sec. |IILC| 

The physics is much richer when interactions between 
polaritons are strong and induce sizable quantum cor- 
relations in the fluid. A simple illustrative example of 
this phenomenon is provided by the so-called superfluid 
to Mott-insulator transition predicted for Bose particles 



trapped at the minima of a periodic potential (Fisher 



et al. 1989 Jaksch et at. 1998 ) and experimentally ob- 



served with ultra-cold atoms in optical lattices ( Greiner 



et al. 2002 1 . Depending on the ratio of the hopping en- 



ergy J between neighboring sites and the on-site interac- 
tion energy U , the ground state of the A^-particle system 
at integer filling v = N/M ranges from a coherent state 
describing the superfluid (for J /U ^ 1) 



I*sf) = 




vac) 



to a Mott insulator state (for J/U ^ 1) 



I* 



M 

M.)-n(^i)' 



(95) 



(96) 



While the former state is well captured by a mean-field 
theory based on the GPE, the Wigner representation of 
the latter would involve phase space regions where the 
Wigner function is negative, signaling the importance of 
quantum correlations and preventing the use of semiclas- 
sical stochastic approaches. 

Another illustrative example of strongly correlated 
state that appears in the theory of strongly interact- 
ing Bose fluids is the so-called Tonks- Girardeau gas of 



impenetrable bosons in a one-dimensional geometry ( 01 



shanii 1998). In stark contrast with mean-field states 



where the many-body wavefunction is the product of 
single particle states, the many-body wavefunction of 
a Tonks-Girardeau gas is characterized by zeros when- 
ever two particles approach to each other. A remark- 
ably elegant exact solution to this problem was obtained 
in (Girardeau 1960[ ) using a rigorous mapping between 
the one-dimensional gas of impenetrable bosons and a 
gas of non-interacting spinless fermions: the energy spec- 
tra of the two systems are identical, as are all configu- 
rational probability distributions in real space. On the 
other hand, the signature of the strong correlations ex- 
isting between the bosons are clearly visible in the mo- 
mentum distributions, which remain instead quite dif- 
ferent from the ferniionic ones. Experimental evidence 
of Tonks-Girardeau gases of ultracold atoms trapped in 
one-dimensional geometries was reported in ( ^Kinoshit^ 
et al. 2004 Paredes et al. 2004 1 : the strongly-correlated. 



fermionized nature of the strongly interacting Bose gas 
was assessed from macroscopic observables such as the 
cloud size, its internal energy and the momentum distri- 
bution of the atoms. We refer tolCazalilla et aLl|2011|for 



a recent review of the physics of one-dimensional Bose 
gases in both the condensed-matter and the atomic con- 
texts. 

In this section we will review the recent theoretical 
advances towards the realization of strongly correlated 
gases of photons or polaritons. A complementary point 



of view on this same physics can be found in ( Houck 



et al. 



callec 



et al 



2012). The basic block in this direction is the so- 



photon blockade effect first proposed in ( Imamoglu 



19971: in the presence of a strong enough opti- 



cal nonlinearity, photons in single-mode cavities behave 



as effectively impenetrable particles ( Birnbaum et al. 
20051 [Faraon et al.\ [2008| [Lang et al.\ |2011[ [Reinhajd 



et al. 2011). When a system of many cavities arranged 



in a lattice structure is considered, photon-photon in- 
teractions are expected to drive the photon gas from a 
coherent, superfluid state to a correlated Mott-insulator 
state. The observability of this phase transition in real 
istic devices was first theoretically investigated in ( An- 
gelakis et al.[ |2007[ [Greentree et al.\ |2006[ |Hartmann 



et al. 2006 2008) following the analogy with material 



particles, namely assuming a quasi-equilibrium condition 
and neglecting the driven-dissipative nature of the pho- 
ton fluid. A related theoretical proposal to generate a 
Tonks-Girardeau gas of fermionized photons in an optical 
fiber appeared in ( Chang et al. . 2008 ) . The next crucial 



step was to include in the model photon losses and the 
pump mechanism that is used to replenish the photon 
gas. While these issues were well known in the litera- 
ture on photon hydrodynamics reviewed in the previous 
sections, the first works addressing and possibly taking 
advantage of the non-equilibrium nature in the context 



of str ongly correlated photon gases we re ([Gerace et al. 
20091) and soon after (ICarusotto et al 



2009) 



A first experimental realization of large optical nonlin- 
earities at the single photon level leading to a strongly 
correlated stream of photons was reported very recently 



in (Peyronel et al. 2012) using the dramatically rein- 



forced optical nonlinearity of coherently dressed atoms in 



a Rydberg EIT configuration ( Sevingli et al. 2011 Dudin 



Fleischhauer[ |2012j [Shahmoon et a/.[|2011 | 



and Kuzmich[ 2012; Gor shkov et aL[|2011 Petrosyan and 



A. Non-equilibrium Bose-Hubbard-like systems 

Many of the recent studies of the physics of strongly in- 
teracting photon gases are based on lattice models where 
photons are trapped in an array of optical cavities. A few 
examples of such configurations in a two-dimensional ge- 
ometry are schematically illustrated in Fig jS] using an ar- 
ray of photonic crystal cavities (central panel) or an array 
of superconducting transmission line resonators (lower 
panel). At a simplest level of description, the system 
can be modelled by a generalized driven-dissipative Bose- 
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Hubbard Hamiltonian of the form 



resulting Jaynes- Cummings- Hubbard lisLmiltonian reads 



H — Ti-BH + T~Lf + T~Lenv — 



[F^ (t) 4 + F* (t) C,] + Henv , (97) 



where each cavity is assumed to support a single photon 
mode whose destruction (creation) operator is Ci (c|^). 
Photon tunneling through the non-perfectly reflecting 
cavity mirrors are responsible for the hopping amplitude 
between neighboring cavities with an amplitude J. The 
photon-photon interactions are included by a two-particle 
term of proportional to the on-site energy HU propor- 
tional to the x^ ^^ third-order optical nonlinearity of the 
cavity medium ( Drummond and Walls 1980 Ferretti and 



Gerace 20121, 



U 



4eoVcff £2 



(98) 



where VcS is the effective volume of the cavity mode, 
X^'^^ and e are the spatially averaged relative dielectric 
constant and optical nonlinearity of the cavity material. 

The driven-dissipative nature of the photon system is 
accounted for by the terms on the third line of ( 97 1 : co- 



herent excitation by a pump laser field is described by the 
site- and time-dependent external field Fi{t), while cou- 
pling with the environment is described by Henv Most 
works have considered an environment consisting of a 
zero-thermal radiative bath into which photons from each 
site can escape at a rate T. Once the environment degrees 
of freedom have been traced out, the corresponding term 
in the master equation for the reduced density matrix of 
the system has the Lindblad form 



dp 
'dt 



1 

ih 



(99) 



A different but related form of Lindblad term in the mas- 



ter equation (99) was used in (Hoffman et al. 2011) to 
describe the driving of the cavity array by a broadband 
incident radiation. Of course, incoherent pumping of the 



cavities by an amplifying bath as discussed in Sec III.D 
may be another interesting option, but so far it has not 
been investigated yet in the context of strongly interact- 
ing photons. 



The driven-dissipative Bose-Hubbard model of (97) is 
just an example of a wide class of models describing sys- 
tems of strongly interacting photons. Another promi- 
nent model is obtained by replacing on each site the 
two-particle interaction term with a Jaynes-Cummings 
interaction term with a two-level emitter of transition 
frequency ojog strongly coupled to the cavity mode. The 



+ Sa;eg|e),(e|, + hg (^c^\e),{g\i + cj|g),(e|,) (100) 

With |g)j and |e)i, we indicate the ground and excited 
states of the emitter on site i. The coupling between 
the photon mode and the two-level system is quantified 
by the vacuum Rabi frequency g. Hamiltonians of this 
form can be realized in the physical systems shown in 
FigjS] by coupling a semiconductor quantum dot to a 
photonic crystal cavity (middle panel) or a superconduct- 
ing Josephson artificial two-level atom to the microwave 
transmission line resonator (lower panel). Overall, most 
of the predictions of the Jaynes-Cummings Hamiltonian 
( 100 ) are in close qualitative agreement with the ones of 



the simpler Bose-Hubbard one ( 97 ) ; some remarkable dif 



ferences concerning the many-body physics of the photon 



and in (Grujic et al. 2012 



L( 


Schmidt and Blatter 


2009 


2010 


2012 





Another variation of the generalized Bose-Hubbard 
model involves a cavity mode strongly coupled to a non- 
linear matter excitation, 



^ + elk ) + 



+ubrbibA+j2f'j4cj- (101) 

{id} 

where Swexc is the energy of the matter excitation created 
by the bosonic operator b] and flu is the vacuum Rabi 
frequency of the light-matter coupling. This polariton- 
Bose-Hubbard model naturally arises in the description 
e.g. of a semiconductor micropillar cavity mode strongly 
coupled to an quantum well exciton transition: as dis- 
cussed in Sec |II.C.3l strong interactions arise from the 
exciton-exciton interactions. In the spatially confined 
geometry, the ir strength can be q uantified by the single 
parameter U (Verger et al. [2006 ). 



Provided the energy scale set by the vacuum Rabi 
frequency is larger than the nonlinear interaction, the 
pumping and the losses, one can restrict the descrip- 
tion to a single polariton mode and recover a Bose- 
Hubbard Hamiltonian of the form ( 97 ) with rescaled hop- 
ping and nonlinear interaction parameters, UcS = [/ 
and Jeff = |C|2 J. 



B. Photon blockade effects 

The first example of nonlinear optical process at the 
single photon le vel is the so-called photon blockade effect, 
anticipated in (Carmichael 1985[ ) and fully elaborated 
in (Imamoglu et al. 19971 in analogy with the Coulomb 



blockade effect of electron transport through mesoscopic 
devices (iFulton and Dolanl 119871 iKastnerl 119921). 
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The basic idea is to consider a single-mode cavity (i.e. 



a single site of the Bose-Hubbard Hamiltonian (97)) with 
a strong enough optical nonlinearity to be in the J7 ^ 7 
regime. In this case, the two-photon state is spectrally 
shifted of an amount 2U much larger than its linewidth 
27 and a laser field resonant with the linear cavity fre- 
quency loq will not be able to inject into the cavity more 
than one photon at a time. In physical terms, the pres- 
ence of a single photon in the cavity is able to effectively 
block the entrance of a second one: as a result, the trans- 
mitted light across the cavity show a strong antibunching, 
closely analogous to the one observed in the resonance 



fluorescence by a single two-level atom (Kimble et al. 



1977). The figure of merit of this antibunching is the ra- 
tio U/T, with the zero-delay photon correlation function 
5(2)(0) = {c^c''cc)/{c^cf = l/(l-K4(C//7)2) quantifying 
the ratio of having two photons on a single site. General- 
ization of the photon blockade to nonlinear atom optics 



using matter waves was proposed in ( |Carusotto 2001 



Kolomeisky et al. 2004). 



1. Jaynes-Cummings nonlinearities 

The main ingredient to observe the traditional pho- 
ton blockade effect is a large difference between the 
two-excitation state energy En=2 and 2En=i, i.e., twice 
the energy of the single-excitation state, regardless of 
the physical origin of the nonlinearity. In particular, 
\En=2 — 2-E'n=i I ^ fi.7 is a sufficient photon blockade con- 
dition, being 7 the loss rate. Indeed, the photon blockade 
effect turned out to be experimentally easier to observe in 
Jaynes-Cummings systems where the optical nonlinear- 
ity is due to a single two-level emitter strongly coupled 
to the cavity mode. 



As it is shown in the left panel of Fig 29 the level 
scheme of the Jaynes-Cummings model consists of an in- 
finite series of doublets above a non-degenerate ground 
state. When the emitter is exactly on resonance with the 
cavity mode, the splitting of the n"^ doublet is equal to 
g being the emitter-cavity coupling energy. For a 
pump field on resonance with, e.g., the lower state of the 
n = I one-excitation doublet, the closest two-excitation 
state is then detuned by a frequency amount (2 — ^/2) g. 
As a result, a sufficient condition to observe significant 
antibunching in the transmitted light is that the vacuum 
Rabi frequency g is much larger than the loss rates of 
both the cavity mode and the emitter. 

In experiments, a sizeble antibunching was observed 
in the optical domain using either a single atom in a 



macroscopic optical cavity (Birnbaum et al. 2005), or 



to a slightly less extent- using a single quantum dot in 



a photonic crystal nanocavity (Faraon et al., 2008 Rein- 



hard et al. 2011). This latter case is illustrated in the 



right panels of Fig 29 The most spectacular results were 
obtaine d in the microwave domain of circuit QED exper- 
iments (Lang et al. 2011), where the large g/T ^ 100 
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FIG. 29 Left panel: Sketch of the level structure of the 
Jaynes-Cummings model with a laser pump on resonance with 
a one-photon transition. Figure from ( Birnbaum e£ aZ. 20051. 
Right panels: second-order correlation g^'^\Q) for a single 
quantum dot in a photonic crystal microcavity. Upper panel: 
Theoretical calculations for (/'^'(O) as a function of cavity- 
exciton and laser detuning. Blue (red) regions correspond 
to sub- (super-) Poissonian statistics. Lower panels: Experi- 
mental results (points) compared with theoretical curves for 
specific values of the cavity detuning. Figure from |Reinhard| 
\eral\^ffn\ 



the right panels of Fig|36j An experimentally different 
strategy to assess photon blockade effects was adopted 
in (Hoffman et al. 2011) using a broadband microwave 



excitation instead of a coherent drive. 

An interesting new perspective on photon blockade ef- 
fects was put forward by |Shen and Fan |2007| instead 
of the usual quantum optics approach of studying the 
optical response of a coherently driven two-level emitter, 
they applied Bethe ansatz techniques to solve the quan- 
tum dynamics of a two-photon state in terms of the two- 
photon wavefunction. In particular, they extracted from 
the scattering S matrix interesting predictions for bunch- 
ing and antibunching features in the scattered light. 



2. Resonant optical nonlinearities in multi-level systems 

It is common wisdom in nonlinear optics that the 
strongest optical nonlinearities are found in the vicinity 
of optical resonances, where unfortunately also absorp- 



leads to the very pronounced antibunching illustrated in 



tion is ge nerally large (Boyd 2008 Butcher and Cot- 
ter 20081. The devices in the strong light-matter cou- 
pling regime that we have mentioned so far overcome 
this difficulty using very narrow transitions with weak 
non-radiative losses. Another, different strategy is to fol- 
low the pioneering proposal in (Harris et al. 1990| based 
on Electromagnetically Induced transparency (ElT) me- 
dia (Arimondo 1996 Fleischhauer e< a/. 2005[ ): by suit- 
ably dressing an otherwise strongly absorbing three-level 
medium with a coherent dressing field, a subtle destruc- 
tive interference effect turns out to suppress one-photon 
resonant absorption, while keeping a large matrix ele- 
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the atomic context where strong correlations are induced 
in cold molecular gases by a strong dissipation (Diirr 



et al. 2009 Syassen et al. 2008), it was realized that 



gi3 y wi 



FIG. 30 EIT scheme to enhance the strength of optical non- 
linearity in a coherently dressed four level atom. The cavity 
mode at frequency ujc is coupled to the 1 — >■ 3 and 2 4 
atomic transitions, the coherent field dressing the atom on 
the 2 — >■ 3 transition has a Rabi frequency Q. . The parame- 
teres g\s and 1724 quantify the respective dipole couplings and 
5, A and e are the detunings. Figure from (Hartmann et al. 



an effective blockade effect without significant losses can 
also be observed via a large Im[?7]: the underlying mech- 
anism is based on the so-called Zeno effect, which effec- 
tively forbids a second particle from entering a site if the 
two-photon loss rate exceeds tunneling. A possible imple- 
mentation of this dissipative blockade idea may involve 
the strong dissipative nonlinearity ( |102[ ) of EIT systems 
in the A = regime, as discussed in (Harris and Ya- 



mamoto 



19981. An alternative scheme to reinforce the 



photonic nonlinearity was proposed in (Brandao et al. 
20081). 



2006) 



ment for the two-photon transitions that are responsible 
for the optical nonlinearity. 

The idea of using a single-mode cavity filled with 
an EIT medium was at the heart of the first pro- 



Even though EIT allows to suppress absorption losses, 
the resulting value of the resonant optical nonlinear- 
ity remains of the same order as the one of bare two- 
level atoms. A strategy to further increase the ef- 
fective photon-photon interactions using EIT with a 
pair of counter-propagating dressing fields was proposed 
in ( Andre et al. 2005 1 . As it was experimentally demon- 



strated in (Bajcsy et al.. 2003), this configuration allows 



posal for photon blockade (Imamoglu et al. 1997 (see 



also Grangier et al. 1998 and Imamoglu et al. 119981). 
Under suitable assumptions on the light-matter cou- 
pling parameters indicated in level diagram of Fig |30[ 
a large on-site nonlinear interaction energy can be found 



Imamoglu 19991, 



one to take advantage of Bragg scattering processes to 
freeze light as a stationary excitation within the EIT 
medium while still keeping a sizable photon fraction in 
the polariton. 

Another strategy to enhance photon-photon interac- 
tions involves the strong dipole-dipole interactions of Ry- 
dberg atoms, i.e. alkali-like atoms with a single electron 



in the form ([Hartmann and Plenio[ |2007[ [Werner and| promoted into a highly excited orbital ( Saffman et al. 



EIT 



(102) 



The real part of ( 102 ) is responsible for the detuning of 
the two-photon state: the blockade effect sets in as soon 
as the on-site interaction energy U exceeds the effective 
decay rate "/bit of the filled cavity. As it was experi- 
mentally demonstrated in (Laupretre et al. 2011), the 



2010): when the third state of the EIT lambda configu- 
ration consists of a Rydberg state, the presence of a sin- 
gle photon in the system is able to drive away from the 
Rydberg-EIT resonance condition all atoms contained 
in a surrounding volume of mesoscopic size. For large 
enough optical densities, a second photon traveling across 



dramatically reduced group velocity of EIT media is re- 
sponsible for a corresponding suppression of the cavity 
decay rate by the same factor^ 



fer a 


sizable phase shift ([Sevingli et al. 2011 


Gorshkov 


et al. 


20111 Petrosyan and Fleischhauer 2012 


Shahmoon 


et al. 


2011[), which can be exploited to obtain an effective 



blockade effect. 

All these schemes to enhance the optical nonlinearity 
are however not limited to discrete geometries involv- 



lEIT 



7 



1 + 



(103) 



From (102), it is immediate to see that a small value 
of A is favorable to reinforce Ueit'- a naive upper limit 
to Ueit is set by t he tw o-photon losses described by the 
imaginary part of (102). This statement was partially 



ing cavities: recent theoretical (Chang et al. 2008) and 
experimental works have investigated fiber geometries 
with a gas of atoms coupled to strongly confined opti- 
cal modes propagating along the fiber with mode waists 
as small as 100 nm: to this purpose, hollow-core pho- 



tonic crystal fibers have been used (Bajcsy et al. 2009 



Venkataraman et al. 2011), as well as standard optical 



overturned by the recent works in (Hafezi et al. 2012 



fiber strongly tapered by mechanical stretching (Vetsch 



Kiffner and Hartmann 2010 ): inspired to related work in 



^ Note that the use of detuned EIT media to increase the group 
velocity above c has a much more complicate effect on the cavity 



decay dynamics (Laupretre et al. 2012 1 



et al. 2010 1. In both cases, the nonlinear optical medium 
consists of atoms that are injected in the hollow core of 
the fiber or are optically trapped at its surface. The quite 
high optical depths achieved in this way have led to ef- 
ficient all-optical switching at low powers using a EIT 
sche me ([Bajcsy et aL[|2009| or a two-photon absorption 
one (Venkataraman et al. 2011): the transmission of a 
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FIG. 31 (a) Dispersion of the polariton frequency versus in- 
plane momentum; the exciton-cavity detuning is chosen such 
that the lower polariton is on (Feshbach) resonance with 
the biexciton energy (horizontal dot ted line). Panel (b). 
Schematic diagram of the states belonging to the Ntot = 0, 1, 2 
manifolds exactly on Feshbach resonance 2Elp -\- V = Eb- 
Panel (c). k — anticrossing of the lower and upper po- 
lariton branches as a function of the bare cavity frequency 
Eq. The horizontal dashed line indicates half the biexciton 
energy. Panel (d). Position of the different states belong- 
ing to the Ntot ~ 2 manifold as a function of the bare cav- 
ity frequency E^- Dashed line: bare \pp) and \B) states. 
Solid lines: states originating from the mixing of the biex- 
citon state \B) and the \pp) state. The vertical dotted line 
indicates the Feshbach resonance point. Rightmost panel (e). 
T = value of the second-order optical correlation function 
gf'^>(T) = {p\t)p\t + T)p{t + r)p{t)) as a function of the 
exciton-biexciton coupling G for a system on Feshbach res- 
onance. Parameters: Ex = 1.4 meV; HUr = 1.5 meV; 
G = 0.03 meV. Figure from ([Carusotto eFoTlpOlOl 



atoms is resonant with a long-lived quasi-bound molec- 
ular state. In typical experiments, an external magnetic 
field is used to tune the energy of the quasi-bound molec- 
ular state close to the energy zero of scattering states. At 
this point, the low-energy scattering amplitude diverges 
and a dilute ultra-cold atomic gas has the strongest pos- 
sible interaction. Remarkable recent experiments with 
atomic Fermi gases have exploited such Feshbach reso- 
nances to demonstrate superfluid behavior at tempera- 



tures comparable to the Fermi temperature ( Bloch et al. 



2008 Giorgini et al 2008 ) 



In the wake of the pioneering attempts ( Savasta et al. 



2003 1999), a recent work (Wouters 2007) has pointed 



out the possibility of exploiting an analogous Feshbach 
resonance effect on an intermediate biexciton state to 
enhance the optical nonlinearity of two-dimensional po- 
laritons in planar microcavities. Biexcitons are two 
electron-two hole bound complexes analogous to hydro- 
gen molecules in semiconducting materials and have been 
widely studied durin g the last decades in both bulk and 
confined geometries (Ivanov et al. 1998). As it is de- 



scribed in FigjSlJ a careful tuning of the length of a pla- 
nar microcavity in the strong coupling regime allows to 
bring the energy of a pair of polaritons at the bottom 
of the lower polariton branch close to resonance with 
the biexciton state. This resonance is then responsible 
for a dramatic enhancement of the low-energy scatter- 
ing cross-section for polariton-polariton collisions in the 
singlet channel. An analytic model of the biexciton Fes- 



hbach resonance was put forward in (Carusotto et al. 



2010), where a method inspired by the atom-molecule 
appro ach to strongly interacting degenerate quantum 
gases ( Heinzen et al. 2000 Kokkelmans et al\ [2002 ) is 
used to obtain a approximate closed formula for the po- 
lariton scattering T matrix element in the singlet chan- 
nel, 



probe beam is significantly affected by the presence of 
a control beam containing a number of photons on the 
order of a thousand (for the FIT scheme) and few tens 
(for the TPA scheme). 

Very recently, the large optical nonlinearity in an op- 
tically dressed dense atomic cloud in a Rydberg-EIT 
scheme has allowed for the first experimental claim of 
quantum nonlinear optics at single photon level in a prop- 
agating, cavity-less geometry: the signature of strong 
photon correlation is the conversion of an incident coher- 
ent laser light beam into a non-classical train of single- 
photon pulses (Peyronel et al.. 2012). 



3. Feshbach blockade effect 

A great deal of the recent advances in the field of 
strongly correlated atomic gases have been made possible 
by the so- called Feshbach re sonance effect in atom-atom 
collisions (Chin et al. 2010): the scattering cross section 



\9\ 



E-E% 



"^LP 1^12 ]r.„\Erna, 



(104) 

where E is the collision energy, E'^ is the (bare) biexciton 
energy, g is the polariton-biexciton coupling coefficient, 
and Emax is a UV cut-off energy accounting for the mi- 
croscopic details of the system. On physical grounds, 
the coupling coefficient g is expected to be related to 
the exciton-photon Rabi frequency ^Ir and the exciton 
Bohr radius (or, equivalently, the biexciton size) as by 
q MliiO-B- As usual in many-body theory ( [Petrov 



et al. 2000 1 , the biexciton contribution to the polariton- 
polariton interaction constant in a spatially extended 
two-dimensional geometry is obtained from the T-matrix 
( 104 1 as the amplitude for forward scattering at an en- 



ergy E determined by the oscillation frequency of the 
Bose field. 



is dramatically enhanced when the energy of two colliding 



The same biexciton model ( Carusotto et al. 2010 ) pro- 
vides an estimation of the single-polariton nonlinearity 
coefficient in tightly confined geometries: the predicted 
efficiency of this Feshbach blockade effect to generate a 



53 



non-classical stream of anti-bunched photons from real- 
istic structures is illustrated in Fig 31 'e). Exactly on Fes- 
hbach resonance, one can expect an effective detuning of 
the two-excitation state on the order of Mlj^aB/L (note 
the different scaling on the characteristic spatial size L of 
the cavity mode as compared to (98)). Alternatively, ad- 
ditional decay channels of the biexciton state other than 
the decay into a pair of cavity polaritons can be exploited 
to obtain an efficient dissipative blockade effect. As the 
biexciton mass is not renormalized by the coupling to the 
photon, a possible difficulty is the sizable inhomogeneous 
broadening of the biexciton energy due to short-distance 
fluctuations of the quantum well thickness. 



4. Unconventional photon blockade in photonic molecules 

The discussion of the previous sections pointed out the 
occurrence of effective photon blockade effects as soon as 
the optical nonlinearity U (or the vacuum Rabi frequency 
5' in a Jaynes-Cummings system) largely exceeds the loss 
rate F. This is a sufficient condition to observe strong 
antibunching in the transmitted light, but actually it is 
not strictly necessary. 

This remarkable fact was first investigated in the field 



(a) 



of cavity quantum electrodynamics (Carmichael 1985 



Carmichael et at. 1991): as a result of a peculiar quan- 
tum interference process, a sizable antibunching of trans- 
mitted light is predicted to be observed from cavities 
containing a large number of emitters with moderate val- 
ues of the single emitter-cavity coupling. Experimental 



studies along these lines appeared in (Foster et at. 2000 



Rempe et at. 1991 ). 



Very recently, |Liew and Savona[ |2010| reported a 
numerical study of a driven-dissipative two-site Bose- 
Hubbard system, where a resonantly driven cavity is cou- 
pled by tunneling events to a non-pumped auxiliary cav- 
ity. The remarkable result of this work is that a strong 
antibunching ((/'-^'(O) <C 1) can be obtained with a weak 
Kerr nonlinearity [/ ^ 7. The strongest antibunching is 
achieved for optimal finite values of the on-site repulsion 
energy U and of the detuning between the pump and the 
cavity mode frequency. Independently from the previous 
work in cavity-QED, the origin of the strong antibunch- 
ing was traced back to a sort of quantum interference 
effect. 

The underlying physical mechanism was analytically 
unveiled in a following work (Bamba et al. 2011 ) using a 



wavefunction amplitude approach, which is exact in the 
limit of weak pump limit. Such an approach is based on 
the following ansatz 



1^) = CoolOO) 



(CiollO)-l-CoilOl))- 



'^p*(C2o|20)+Cn|ll)+C7o2|02)) + 



(105) 



to calculate the steady-state of the coupled cavity system, 
where Wp is the pump frequency driving the first cavity. 
Here, \mn) represents the Fock state with m particles in 
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FIG. 32 (a) Optimal nonlinearity (7opt and detuning AiSopt 
for a driven-dissipative two-site Bose-Hubbard system are 
plotted versus the inter-mode coupling strength J normal- 
ized to mode broadening 7 (71 = 72 = 7 and Ei — E2 — E). 
(b) Transition paths leading to the quantum interference re- 



sponsible for the strong antibunching. From (Bamba et al. 
20TT|. 



the first cavity site and n particles in the auxiliary cavity 
(second site). Under weak pumping conditions (Cqo ^ 
C'lo, Cqi ^ ^^20 7 Cii, C02), it is possibl e to determine the 
coefficients C„i„ in a iterative fashion ( Carmichael et al. 



1991) 



Remarkably, there exist values of the parameters for 
which one has C20 = 0, i.e. a vanishing probability of 
having 2 photons in the first cavity: this is indeed what 
is required for an ideal photon blockade. The analytical 
solutions of this condition have a particularly transparent 
form for equal energies and loss rates of the two cavity 
modes {Ei = E2 = E; — j2 = 7) and for J ^ 
7. The condition C20 = is satisfied when the on-site 
interaction U and the pump detuning AE — E — hup 
take the following optimal values: 



pt 



7 
2^3 



and 



2 73 



(106) 



For J ^ 7, it directly implies that Uopt/l ^ 1, i-e., the 
optimal value of the nonlinear interaction energy is much 
smaller than the loss rate. 

It is worth mentioning that this dramatic release of 
the constraint on the U /J parameter to observe photon 
blockade comes at the expense of fast oscillations of the 
photon correlation function ^'■^■'(t) as a function of the 
delay time r: the characteristic scale of these oscillation 
is set by the tunneling rate J ^ F. As a result, photode- 
tectors with a sufficient time resolution are required to 
observe the antibunching effect. 

This unexpected photon blockade effect can be inter- 
preted as arising from the destructive quantum interfer- 
ence between the amplitude of two transition paths hav- 
ing as initial state the vacuum in both cavity sites and 
1 20) as final state, as shown in Fig 32 The interference 
is between the following two paths! (a) the direct exci- 

tation from |10) — > |20) (solid arrow), meaning that two 
photons are sequentially injected by the resonant driv- 
ing field in the first cavity, (b) tunnel-coupling-mediated 
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transition |10) A |01) (|11) A |02)) A |20) (dot- 
ted arrows). For an optimal value of the nonlinearity 
and an optimal value of the pump frequency detuning, 
the two above mentioned transition paths have the same 
amplitude and opposite phases, leading to the destructive 
interference. 

It is worth noting that the nonlinearity Ui of the 
pumped cavity mode is irrelevant for the antibunching. 
This implies that a finite, but weak nonlinearity is re- 
quired only in the auxiliary (undriven) photonic mode to 
achieve the perfect destructive quantum interference that 
leads to C20 — 0. In the same paper (Bamba et al. 2011 ) 
it was shown that an analogous effect also occurs for a 
two-site Jaynes-Cummings-Hubbard system: in contrast 
to the usual single-site photon blockade effect where a 



large vacuum Rabi frequency is needed ( Birnbaum et al. 



2005 Lang et al. 2011 ), the two-site geometry allows for 
a strong antibunching 17*^^) (0) ~ even when (7/7 ~ 1. 
A related proposal has recently appeared in ( [Majumdar 



et al. 2012) 



This paradigm of photon blockade with weak nonlin- 
earities has been extended to the case of cavities with 



polarization splittings in (Bamba and Ciuti 2011 1. The 



additional degree of freedom provides a way to tailor such 
quantum interference photon blockade effects. In partic- 
ular, using a polarized pump driving the first cavity, it is 
possible to get photon blockade for the counter-polarized 
mode in the auxiliary cavity, thus providing an efficient 
way to filter spurious pump photons through spatial and 
polarization filters. 

Given the novel physical mechanism underlying anti- 
bunching in weakly non-linear photonic molecules, a pri- 
ori it is not clear whether the predicted antibunching in 
a single cavity geometry may produce strong quantum 
correlations in an array of coupled photonic molecules. 



This important question was addressed in (Bamba et al. 
2011): a geometry consisting of ring of three photonic 



molecules whose the driven cavities are coupled with each 
other by a tunnel coupling. For a suitable choice of pa- 
rameters, a nearly perfect antibunching could be revealed 
also in this case. 



C. Lossless many-photon gases at (quasi) equilibrium 

The first theoretical investigations of many-photon ef- 
fects in arrays of coupled cavities appeared almost si- 



multaneously in ( Angelakis et al. 


2007 


Greentree et al. 


2006 


Hartmann et al. 


2006). The basic idea of all these 



works closely follows classical approach of many-body 



physics with material particles (Fisher et al. 1989), as 
recentl y anticipated and then observed with ultracold 
atoms ( Greiner ei aZ. 2002 Jaksch aL 1998): depend- 



ing on the ratio between tunneling J and interaction U 



energies in the Bose-Hubbard Hamiltonian T-Lbh in (|97|), 
the ground state of the system exhibits a quantum p 



liase 



transition from a superfluid state ( 95 ) with spatially ex- 



coherence and suppressed number fluctuations on each 
site. 

In the first work ( [Hartmann et al. 2006), the optical 
nonlinearity required for the blockade effect is provided 
by a gas of optically driven four-level atoms in an Elec- 
tromagnetically Induced Transparency scheme inserted 
in t he cavities as illu s trated in FigpO} in the two o th- 
ers (Angelakis et al. 2007 Greentree et al. 2006), a 



Jaynes-Cummings-Hubbard model is considered with a 
simple two-level emitter strongly coupled to each cavity 
mode. The envisaged protocol is to tune the system in 
the deep Mott phase at large U/ J and prepare the sys- 
tem in a state with a single photon per site. Then, U /J is 
adiabatically ramped down into the superfluid phase so 
that the system remains at all times in its ground state. 
A crucial advantage of the photon system over the atomic 
one is that the onset of the phase transition can be as- 
sessed in situ e.g. by measuring the number of photons 
in each cavity and its fluctuations. However, in all these 
works the photon gas is considered as a quasi-closed sys- 
tem, where loss processes are just an hindrance limiting 
the available time for the adiabatic ramp. 



A different geometry was considered in (Chang et al. 



2008): slow photons in a Electromagnetically Induced 
Transparency regime propagate along a one-dimensional 
fiber. By tuning the intensity of the EIT dressing beams, 
one can tune the strength of the effective repulsive in- 
teractions between the photons, as well as reduce their 
propagation speed along the fiber to zero. For sufficiently 
strong interactions, the stored photons within the fiber 
behave as a Tonks-Girardeau gas of impenetrable parti 



cles ( 


Girardeau 


jcsy et al. 2003 



1960). As it was demonstrated in (Ba- 
, the use of a pair of counterpropagating 
EIT beams allows to tune also the photon propagation 
speed along the fiber and to reduce it to zero. 



The experimental protocol that is proposed in ( Chang 



et al. . 2008 ) is then to inject a stream of fast and coher- 
ent photons into the fiber in a weakly interacting state. 
Once photons are loaded, the propagation speed is re- 
duced to zero so that the gas of photons gets stored 
within the fiber. Simultaneously, the interaction con- 
stant is adiabatically ramped up, possibly achieving the 
Tonks-Girardeau regime. As in the previously quoted 
work, the ramps in the interaction constant and/or the 
photonic mass have to be performed slow enough for the 
gas to remain with good fidelity in the ground state. An 
exciting advantage of the fiber geometry is that the in 
situ density correlation function of the photon gas within 
the fiber can be measured by switching the photons again 
into a propagating state and detecting the arrival times 
on an external photodetector. 

Following these pioneering proposals, a number of 
other works have appeared discussing the possibility of 
studying more complex many-body effects in photon 
gases, e. g. the mapping of spi n models onto photonic 
systems (Hartmann et al. 2007), the emergence of Bose 



tended coherence to a Mott insulator state ( 96 ) with no 



glass phases due to the interplay of strong interactions 
with disorder (Rossini and Fazio 2007) (a different mean- 
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field approach to Bose glass phases in weakly interact- 
ing polaritons systems was reported in |Malpuech et al'. 



Ilk 



20071, spin-charge separation effects in two-component 
photon gases (distinguished by their polarization state or 
their frequenci es) (Angelakis et al. 2011 1, a kind of BEC- 



states ( Cho et al. 2008 1 . A recent review of this rich liter- 



BCS crossover ( Huo et al. 201 1 1 , fractional quantum Hall 



ature can be found in ( Angelakis| 2012 Hartmann et al. 
[2008). In all cases, a quasi-equilibrium picture was used, 
where the photon gas is modeled as a closed and conser- 
vative system and the standard concepts of many-body 
physics apply, as in liquid Helium, electronic systems and 
atomic gases. 



D. Non-equilibrium strongly correlated photon gases: 
Josephson interferometers and Tonks-Girardeau gases 

In the last decades, an intense activity has been de- 
voted to the study of the physics of system far from 
thermodynamical equilibrium (Ruelle 2004). In partic- 
ular, a richer variety of many-body behaviors and phase 
transitions was predicted and in some case observed as 
a result of the interplay of driving, dissipation and inter- 
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actions (Mukamel 1999). In this context, a special was 
recently devoted to the engineering of open systems of ul- 
tracold atoms and ions for quantum simulations of equi- 
librium and non-equilibrium quantum phases (Miiller 



et al. 2012) 



While the driven-dissipative nature of optical systems 
is naturally taken into account in the quantum optical 
literature (Carmicheal 1993), most of these works have 
focused their attention on few mode systems: in spite 
of the very restricted spatial dynamics, a rich wealth of 
quantum optical features was already predicted and ob- 
served. Exception are perhaps the pioneering literature 



U/J 



FIG. 33 Upper-left panel: a sketch of an array of 5 nonlinear 
cavities in a ring geometry with periodic boundary conditions. 
The system is driven by a coherent laser and the properties 
of the strongly correlated photon system can be monitored, 
for example, by the spatial-dependent second-order coherence 
function g^. Upper-right panel: steady-state mean number of 
photons as a function of the pump detuning Aujp — ujp ~ lvq 
in the impenetrable boson limit U/J — > -l-oo. The different 
curves correspond to increasing values of the pump amplitude 
Fp/"f — 0.1, 0.3, 1, 2, 3. The vertical dotted lines indicate the 
spectral positions of the peaks predicted by the fermioniza- 
tion procedure with different number of photons. The kets 
|gi....gjv) represent the pseudo-momenta gi...gjv of the oc- 
cupied fermionic orbitals according to Girardeau's the Bose- 



Fermi mapping ([Girardeau 19601. Lower panel: two-photon 



on quantum solitons (Drummond et al. 


1993 


Kartner 


and by different sites (red) when the interaction parameter is 


and Hans 


1993 


Lai and Hans 1989a|b ) and the on-going 


ramped from the weakly- to the strongly-interacting regime. 



research on spatially multi-mode squeezing effects and 
the consequent possibility of sub-shot noise measurement 



of nano-displacements (Treps et al. 2003 2004). 



The framework is even richer when the spatially ex- 
tended geometry is combined with strong photon-photon 
interactions. A pioneering theoretical investigation of 



this physics appeared in ( Gerace et al. 2009 ) for a sim 



plest geometry consisting of a string of a few single-mode 
cavities coupled by tunneling processes: the strong on- 
site interactions result in strong photonic correlations in 
the non-equilibrium steady-state of the system, which are 
visible through the suppression of Josephson-like oscilla- 
tions, as well as through the photon statistics of emitted 
light. 

A different point of view on non-equilibrium strongly 
correlated photon gases was proposed in (Carusotto 



For each value of U /J, the pump frequency is chosen on reso- 
nance with the lowest energy two-particle state. Figure from 
ICarusotto etal\ |2009| 



Tonks-Girardeau nature of the steady state were iden- 
tified in observables as simple as the transmission or ab- 
sorption spectrum of the device: each of the resonant 
peaks that are visible in the spectra in the upper-right 
panel of Fig |33| corresponds indeed to a different strongly 
correlated many-body state of the impenetrable photon 
gas. The labels on the vertical lines indicate the structure 
of each Tonks-Girardeau state in terms of the correspond- 
ing Fermi orbitals according to Girardeau's Bose-Fermi 
mapping (Cazahlla et al. 2011 Girardeau 1960). 



et al. 



2009 1 , where an array of many cavities in the im- 
penetrable photon regime under a coherent, continuous- 
wave pumping was considered, as sketched in the upper- 
left panel of Fig 33 Unambiguous signatures of the An example is shown in the lower panel of Fig 33 for a 



For a pump frequency on resonance with a given many- 
body eigenstate, the microscopic properties of the corre- 
sponding many-body wavefunction can be inferred from 
the coherence functions of the secondary emitted light. 
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pump resonant on the lowest two-body state of a 3 sites 
lattice: depending on the ratio U /J, one moves from a 
weakly interacting regime where the two-body correlation 
functions on the same site and for different sites coincide, 
to a regime where strong on-site anti-bunching coexists 
with strong bunching on distinct sites. In (Umucahlar 



and Carusotto 2012 (, this same approach was applied to 
fractional quantum Hall states of impenetrable photons 
in the presence of an artificial gauge field. 

Different aspects of the steady state of an array of 
nonlinear and lossy optical resonators driven by coherent 
lasers were addressed in ( Hartmannl 2010 ). In particular, 



for weak driving intensity the steady state was shown to 
be dominated by interactions in such a way that photon 
crystallize into dimers localized on neighboring sites and 
anticorrelations appear between distant sites. Much of 
these results were obtained by applying for the first time 



to the optical context the driven-dissipative version (Ver- 



straete et al] |2004[ [Zwolak and Vidal^ 2004, ) of the novel 
numerical techniques based on the time evolving block 
decimation algorithm ( Schollwock 2005). More extended 



studies along these lines have appeared in (Kiffner and 
Hartmann 2011 Leib and Hartmann 2010 ). A detailed 



comparison of non-equilibrium many-photon physics in 
Bose-Hubbard and Jaynes-Cummings-Hubbard models 



was recently reported in (Grujic et al. 2012) using so- 



phisticated numerical techniques. 

Quantum transport of strongly interacting photons 
along a continuous one-dimensional fiber was studied 
in ( )Hafezi et al. 2012 2011a I. In contrast to most pre- 
vious work, photons are assumed to freely move along 
the fiber in the absence of any lattice potential: the the- 
oretical description is carried out in a first-quantization 
picture based on a set of n-photon wavefunctions ipn- 
For each value n of the total photon number, the corre- 
sponding wavefunction ipn{xi, . . . , is a function of the 
n photon coordinates Xi, . . . , a;„ and is connected to the 
neighboring ones V'n±i by terms describing the photon 
injection and photon loss processes. In the present fiber 
geometry, these processes reduce to suitable boundary 
conditions at the two ends of the fiber. The significant en- 
ergy shift of the two particles states due to interactions is 
responsible for significant bunching and anti-bunching ef- 
fects depending on the sign of interactions and the pump 
laser detuning. For attractive interactions, propagating 
two-photon bound states can be observed on the trans- 
mission resonances. 



E. Topological states and artificial gauge fields 

The motion of charged particles in magnetic fields is 
a key paradigm of quantum mechanics and underlies a 
number of intriguing phenomena in very different con- 
texts, from magnetohydrodynamics in astro- and geo- 
physics to the fractional quantum Hall effect in solid- 
state physics p3as Sarma and Pinczuk| |1997| |Yoshioka[ 
1992 ): in all these systems, the charged particles are sub- 



ject to a real magnetic field generated, e.g., by an elec- 
tric current or by permanent magnets. An interesting 
development of the geometric ideas underlying the theo- 
retical explanation of the quantum Hall effect has led to 
the discovery of a number of other materials with exotic 
properties due to the topological structure of the electron 
bands, the so-called topological insulator and topological 
superconductors ( |Qi and Zhang 2011). 

Independently from these advances of solid state 
physics, a related concept of gauge structure was intro- 



duced in the concept of molecular physics (Mead 1992 



Mead and Truhlar 1979 Resta 2000): the geometric 



phase that may appear in Born-Oppenheimer calcula- 
tions as a result of an adiabatic motion across config- 
uration space can be described in terms of an effective 
(abelian or non-abelian) gauge field. From a more ab- 
stract point of view, this is an example of the Berry 
phase (rotation) that is accumulated by a quantum sys- 
tem during an adiabatic sweep of some external parame- 
ter ( [Berry 1984) and can be described in terms of a U{1) 
(or higher) gauge field (Wilczek and Zee 1984). 

This idea was transposed to the field of atomic physics 



in ( Dum and Olshanii 1996 ) : when an atom is coher- 



ently dressed by a spatially and/or temporally varying 
optical and/or magnetic field, the adiabatic motion of 
its dark states experiences an effective gauge field due 
to the Berry phase. The high fiexibility of atomic sys- 
tems gave a novel twist to this research, as it became 
possible to engineer artificial (or synthetic) gauge fields 



for neutral atoms with almost arbitrary shapes ( Dalibard 



et al. 2011 ). An alternative way of study ing this physic s 
involves rapidly rotating atomic clouds (Cooper 20081: 



modulo the centrifugal force, the Coriolis force that ap- 
pears in the rotating frame of reference is in fact mathe- 
matically equivalent to a Lorentz force in a uniform mag- 
netic field. 

A striking experimental implementation of these ideas 
was reported in (|Lin et al. 2009), where the nucleation 



of a few quantized vortices in a Bose-Einstein condensate 
under the effect of an artificial gauge field was demon- 
strated. More recent experimental investigations have 
addressed the synthetic electric field generated by a time- 
dependent synthetic magnetic field (Lin et al. 2011 1, the 



atomic ba nds in an optical lattice under a strong mag- 
netic field ( Aidelsburger et al. 2011 ), a s uperfluid version 
of the Hall effect (LeBlanc et al. 2012), and a universal 
method to create artificial gauge fields by driving in time 
the lattice potential (Struck et al. 2012). When the arti- 



ficial gauge field is combined with strong atom-atom in- 
teractions, the atomic gases is expected to form strongly 
correlated gases that closely remind quantum Hall liq- 



uids, see e.g. (Chang et al. 2005 



Hafezi ar| |2007[ [Palmer et al. 



Oktel 2007) 



Cooper et al. 2001 



2008[ lUmucahlar and[ 



In parallel to these exciting advances, analogous stud- 
ies have been undertaken in the photonic context: the 
basic idea was to look for photonic systems where the 
orbital motion of the photon is subject to an artificial 
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FIG. 34 Topological photonic states and artificial gauge fields for light. From left to right. Two left-most panels: Photonic 
chiral edge states and effects of a large scatterer. Finite-element simulation of field propagation (left) and transmission spectra 
upon inclusion of the obstacle (right). Panels adapted from Wang et al. 2009 Third panel: Scheme of the circuit QED 



architecture to generate an artificial gauge field for microwave photons. Panel from ( Koch et al.\\2010\. F ourth panel: Sketch 



Pane 



of the photonic device to generate an artificial gauge field for visible photons, 
panel: Hofstadter butterfiy of states in the transmission spectrum of a 10 x 10 two-dimensional array of cavities as a function 



from (Hafezi et al. 201 lb I. Right-most 



of the pump frequency Aujp and of the magnetic flux per plaquette a. Panel from ( |Umucalilar and Carusotto 2011 1 



gauge field. The first proposal in this direction (Hal- stead happen in a standard waveguide. Soon later (Poo 



dane and Raghu 



sidered the two-dimensional photonic bands in an hexag- 
onal array of dielectric rods showing a sizable Faraday 
efi^ect breaking time-reversal invariance. Analogously to 
the electronic case, every photonic band is characterized 
by a topological invariant known as the Chern number: 
when two materials with difi'erent total Chern number 
are juxtaposed, a chiral state appears, which propagates 
along the interface with a unique directionality and is 
immune from scattering against disorder. In two di- 
mensional electron gases, these states play an impor- 
tant role in the integer quantum Hall eff'ect. In the 
photonic case, they form a sort of one-way waveguide 
along which light can propagate in one direction only, 
with no possibility of back-scattering at bends or imper- 
fections. In particular, the absence of backwards propa- 
gating states and of bulk states at the same energy pre- 
vents light from being scattered by the obstacle, so that 
its net effect is limited to a phase shift. Other propos- 
als to generate photonic topological insulators without 
breaking time-reversal made use of arrays of toroidal mi- 
crocavities (Hafezi et al. 2011b I and and metamaterials 



2008 Raghu and Haldane 2008) con- et al.. 2011), an important step forward was the demon- 



stration of unidirectional propagation at the interface of 
a magnetic photonic crystal with air, without the require- 
ment an ancillary cladding layer. Closely related studies 
of one-way waveguiding effects and non-reciprocal behav- 
iors appeared in ( Yu et al. 2008 ) using plasmon modes at 



the surface of a metal under a strong magnetic field and 
in (Lira et al. 2012) using a pair of single-mode waveg- 
uides under an AC-electrical drive. Other features due 
to non-reciprocal behavior in magneto-optical microwave 
devices in the presence of dissipation effects arc discussed 



in the series of works ( jPietz et al.\ |2007| |201lj ). 

Unfortunately, macroscopic photonic crystal devices as 



the ones in (Wang et al. 2009 ) seem unable to reach the 



regime of strong photon-photon interactions required for 
a fractional quantum Hall liquid of photons. The first 



work that has addressed this question is (Cho et al. 



2008), where two-dimensional arrays of coupled optical 



with strong magneto-electric coupling (Khanikaev et al. 
20T2|. 



cavities confining single atoms were theoretically consid- 
ered. A gas of impenetrable bosons is encoded in the 
spin state of the atoms in their s = 1/2 ground state; the 
artificial gauge field is obtained by suitably dressing the 
atoms with laser fields so to obtain a non-trivial hopping 
phase in the generalized Bose-Hubbard model. In 



The idea of topologically protected one-way waveg- Particular, it was shown how to control 0,, via the rela 



uides was experimentally implemented in (Wang et al. 



2009) using a photonic crystal device in the microwave 
range with a square lattice geometry as proposed 



in (Wang et al. 20081: time- reversal invariance was bro- 



ken by the strong gyrotropic permeability of ferrite rods 
in the vicinity of the ferromagnetic resonance. As it is 
illustrated in Fig |34[ one-way propagation of microwaves 
was assessed by measuring the transmission of light be- 
tween two antennas possibly separated by a macroscopic 
obstacle (left-most panel). Transmission in the two direc- 
tions (second panel from the left) differs by several orders 
of magnitudes and is almost unaffected by the presence of 
an obstacle: the microwaves are able to circumnavigate 
the obstacle without being scattered back, as it would in- 



tive phase of the dressing beams. As an exciting applica- 
tion of this novel example of artificial gauge field, a viable 
protocol to prepare Laughlin states was also proposed in 
the presence of strongly repulsive interactions U/ J — > oo. 
A suitable adiabatic ramp of the system parameters (in- 
cluding an external potential) allows to bring the system 
from a ground state where each of the photons is local- 
ized on a separated site, to a different ground state which 
has an excellent overlap with a Laughlin state. Of course, 
this method of preparing the Laughlin state assumes a 
quasi-equilibrium regime in which the slow adiabatic fol- 
lowing of the ground state is not disturbed by any loss 
process. 

A different architecture to break time reversal symme- 
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try and then obtain an artificial gauge field for photons 
was proposed in (Koch et al. 20101. The idea is to use a 



lattice of cavities where tunnel coupling between neigh- 
bor sites occurs via a circulator element: an implementa- 
tion of such a device in a circuit QED architecture with 
an array of coplanar waveguide resonators coupled by 
three junction Josephson rings connected in series and 
threaded by a static magnetic fields is sketched in the 
third panel of Figj34| Signatures of the gauge field are 



predicted to be visible in the single-photon hopping dy- 
namics in a few site geometry as well as in the photon 
bands in a spatially extended Kagome lattice geometry. 
Further developments along these lines are summarized 
in (Petrescu et al. 2012 1. Another way of generating 
an artificial gauge field for photons by dynamically tun- 
ing the individual frequency of the each cavity was dis- 



cus sed in ( Hayward et al . 
als ([Eckardt e t al. 



2012 ) along the lines of propos- 



iolovsky 2011 1 and an experi 



ment (Struck et al. 20121 in the atomic physics context. 

The interplay of the artificial gauge field with the 
driven-dissipative nature of photons was first theoreti- 
cally investigated for a non-interacting photon system 
in (Hafezi et al. 2011b). The considered system was an 



array of toroidal microcavities coupled by tunnel events 
via the connecting waveguides as sketched in the fourth 



panel of Fig 34 in this configuration, a non-trivial tun- 
neling phase (f)ij can be obtained by means of a suitable 
tayloring of the length of the connecting waveguides. Al- 
ternative strategies to generate artificial gauge fields in 
photonic devices based on a planar microcavity architec 
ture in the visible domain were proposed in ( Umucalilar 



and Carusotto 2011 1. In these systems, a nontrivial tun- 
neling phase between neighboring lattice sites can be in- 
duced by a suitable coupling of the orbital and polar- 
ization degrees of freedom in an optically active medium 
or by making the photon to perform a closed path in 
the space of polarizations during a tunneling event. This 
latter scheme generalizes to evanescent waves the geo- 
metric Pancharatnam phase that appear when photons 
propagate across a sequence of birefringent slabs ( Pan- 



charatnam 1956 1 or the Berry phase in twisted optical 
fibers (Tomita and Chiao 1986) and non-planar Mach- 
Zender interferometers (Chiao et al. [1988 ). Its impor- 
tant advantage is the possibility of scaling the transverse 
patterning of the planar microcavity to a micrometer 
scale so to integrate the artificial gauge field with photon 
blockade. 

In these two last proposals, the driven-dissipative na- 
ture of the photons is exploited to propose optical spec- 
troscopic techniques to experimentally assess the effect of 
the artificial gauge field. Remarkably, a regime of large 
magnetic field appears to be under reach, where the loop 
integral J2a 4>ij around a plaquette is of the order of 27r. 
In standard electronic systems such a regime would cor- 
respond to a magnetic flux per unit lattice cell of the 
order of the flux quantum h/e, which requires huge val- 
ues of the magnetic field. In such a strong magnetic field 
regime, the one-particle eigenstates of the conservative 



Hubbard model form the so-called Hofstadter butterfly 
of states (iHofstadterl |1976|): as it is illustrated in the 



right-most panel of Fig|34| many features of this peculiar 
self-similar structure turn out to be observable in quan- 
tities as simple as reflection and transmission spectra in 
relatively small systems. Local access to the optical field 
provides full information on the spatial structure of the 
photon wavefunction. 

In analogy to the earlier work on gyromagnetic pho- 
tonic crystals (Wang et al. 2009), the use of photonic 



lattices under strong magnetic fields to observe unidirec- 
tional propagation of light around the edge of afinite- 



sized structure was proposed in (Hafezi et al. 2011b I. 
As a consequence of the nontrivial Chern number of the 
Hofstadter bands, edge states exist at energies within 
the gap and are topologically protected against imper- 
fections and disorder: in the presence of a defect, light 
wraps around it and then keeps propagating in the same 
direction. These edge states are visible in the transmis- 
sion spectra shown in the right- most panel of Fig[34]as a 
comb of faint lines crossing the gaps. 

The interplay of an artificial gauge field with sizable 
photon-photon interactions was first theoretically inves- 
tigated in (Nunnenkamp et al. 2011 ) with a special eye to 



implementations in circuit QED devices. In spite of the 
simplicity of the three-site geometry considered, an inter- 
esting Schrodinger cat structure was pointed out in the 
few photon eigenstates as well as novel two-photon reso- 
nances in the transmission spectra beyond linear regime, 
along lines very similar to ( jCarusotto et al. 2009j). Ex- 
tension of this work in the direction of the observation 
of Laughlin states when strong artificial gauge fields are 
combined to photon blockade in a square lattice geom- 
etry has been recently undertaken in (Umucalilar and 



Carusotto 2012) 



An alternative way of generating topologically non- 
trivial states of a photon gas was proposed in (Bardyn 



and Imamoglu 2012 1 in analogy to the physics of p-wave 
superconductors: when an open chain of cavities in the 
impenetrable photon regime is pumped by parametric 
pumps injecting a squeezed vacuum in adjacent cavities 



Carusotto et al. 


20091 


Read and Green 


2000 



are expected to appear in the fermionized photon gas and 
can be detected via the mutual second-order coherence 
of the emitted light from the first and last cavities. 



F. Emerging systems for strongly correlated photons 

In this last subsection, we complete the discussion on 
the prospects of strongly correlated photon gases by giv- 
ing a few more details on two physical systems which 
are emerging as most prominent candidates to achieve 
strongly correlated photon phases in arrays of nonlin- 
ear cavities: semiconductor microstructures embedding 
quantum wells or quantum dots, and superconducting 
quantum circuits. Besides these solid-state systems, it is 
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important to remind that a very recent experiment ( Pey- electronic wavefunction in all three-dimensions ( Bastard 



ronel et al. 2012) has demonstrated the potential of 
atomic clouds in the Rydberg-EIT regime to obtain huge 
optical nonlinearities and induce strong correlations in 
photon streams. 



1. Semiconductor micro- and nano-resonators 

In the first Sections of this article, we have reviewed 
the rich physics of dilute quantum fluids of light in spa- 
tially extended two-dimensional planar microcavity ge- 
ometries. Several techniques have been developed to add 
some in-plane potential to confine polaritons in all three 
dimensions or even to create arrays of single mode boxes. 
Some among the most advanced techniques are summa- 
rized in Sec |II.D| Photonic resonator obtained by later- 
ally etching a planar semiconductor microcavity were pi- 
oneered in ( Bloch et al. 1998 ) . More recent samples fab- 



ricated along these lines are illustrated in the left panels 
of Fig. [35] a direct image of the resulting pillar micro- 
cavities is shown in the inset, while the spatial structure 
of the lowest modes is highlighted in the lower panel. 

An alternative concept of photonic resonator in the in- 
frared/visible range is obtained by introducing a defect 
in an otherwise periodic photonic crystal membrane: as 
it happens to electrons in solids, a isolated mode appear 
within the photonic gap, strongly localized in space in 
the vicinity of the defect. An example is shown in the 
central panels of Fig. [35j Most remarkable features of 
these photonic crystal cavities are the extremely high val- 
ues of the Q-factor in the 10^ range and the extremely 
small mode volumes of the order of (A/2)'^ that can be 
nowadays achieved ( Joannopoulos et al. 20081. 



Optical coupling of neighboring photonic resonators 
via the overlap of the evanescent tails has been demon- 



strated for both systems 


Bayer et al. 


1998 Dousse et al. 


2010 


Intonti et al. 


2011 


de Vasconcellos et al. 2011 


). A 



sketch and a review of the main spectroscopical proper- 
ties of a pair of coupled pillar microcavities are illustrated 
in the right panels of Fig|35] Some novel features related 
to Bose-Einstein condensation and Josephson dynamics 
in double cavity geometries have recently been experi- 



mentally addressed in (Galbiati et al. 20111 



In standard photonic devices based on semiconduc- 
tor technology, the optical nonlinearity of the material 
medium used to confine the photon is moderate and the 
device operation is far from the strong interaction regime. 
The most common strategy to approach this regime is to 
embed in the cavity some active material showing a sharp 
material resonance in the spectral neighborhood of the 
photonic mode. In this way, the excitonic content of the 
resulting polariton modes introduces strong interactions 
that can eventually lead to a photon blockade behavior. 

A first successful route in this direction has been to 
insert an electronic quantum dot in a photonic crystal 
resonator: a quantum dot consists of a nanoscopic vol- 
ume of a different material that is able to confine the 



1991 ): as a result, the electronic levels are quantized as in 
a standard atom and discrete optical transitions can be 
isolated. When one of these transitions is strongly cou- 
pled to the photonic mode, a Jaynes-Cummings system 
is recovered: provided the loss rate is weak enough, the 
optical nonlinearity that stems from the anharmonicity 
of the levels shown in the left panel of Figj29] can result 
in an efficient photon blockade effect. 

First experimental evidences of photon blockade by a 
quantum dot embedded in a microcavity driven by a 



coherent pump was reported in (Faraon et al. 20081. 
Later works have shown significantly clearer signatures 
of strong photon interactions (Reinhard et al. 2011), as 



illustrated in the right panels of Fig 29 (i) photon anti- 
bunching upon resonant excitation of the lowest-energy 
polariton state; (ii) photon bunching when the laser field 
is in two-photon resonance with the polariton eigenstates 
of the second Jaynes-Cummings manifold. A subsequent 
work ( Volz et al. 2011 ) has applied this giant optical non- 



linearity to all-optical switching on a ultrafast time scale 

Other inter- 



at the single photon level ( Volz et al. 



esting experimental investigations of 



2011) 



ow-photon- number 



optical switching in microcavities embedding quantum 



dots have recently appeared in (Bose et al. 2012 En 



glund et al. 2012 Loo et al. 2012 K A sophisticated 



four-wave mixing spectroscopy technique to study the 
nonlinear dynamics of these systems was implemented 



in (Kasprzak et al. 2010) 



In standard architectures with a single cavity mode 
coupled to a single quantum dot, an efficient pho- 
ton blockade effect requires to be deeply in the strong 
emitter-cavity coupling regime, i.e., the vacuum Rabi fre- 
quency g largely exceeding the loss rate. This means that 
the cavity Q-factor has to be kept as high as possible 
and, simultaneously, that the quantum dot is correctly 
positioned within the cavity mode profile. Experimen- 
tally, state-of-the-art samples display a coupling to loss 
rate of the order of 3, to be compared to 10 for single 



atoms in macroscopic optical cavities (Birnbaum et al. 



2005 



2011 



and 50 for superconducting circuits (Lang et al. 



The strong coupling condition has been partially 



softened in (Bamba et al. . 2011 1 where the use of a two- 



cavity geometry has been predicted to lead to efficient 
photon blockade even for loss rates comparable to the 
emitter-cavity coupling strength. 

From the point of view of realizing strongly correlated 
states of photons in a many-cavity geometry, quantum 
dot-based architectures suffer from the significant draw- 
back that neither the spatial position of individual dots 
nor the exact position of the electronic transitions can be 
controlled during the growth stage, as the different dots 
self-organize at random positions with random sizes. In- 
deed, in all quoted experiments, an ensemble of quantum 
dots was first prepared in the sample, then the photonic 
crystal cavity was fabricated in a careful way so to have 
the photonic mode located right on top of the quantum 
dot of interest with a frequency precisely tuned to the 
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FIG. 35 Left panels: semiconductor resonators obtained by lithographic patterning of a semiconductor planar microcavity. 
The photoluminescence spectra (upper panel) are characterized by a series of resonances associated to the fully confined modes 
of the squared resonator indicated in the scanning electron micrograph shown in the inset. The bottom panel shows the 
measured and calculated far-field patterns of the photon modes. Figure from (Ferrier et al. [2010 \. Central panels: photonic 



crystal nanocavity embedding a single quantum dot. Figure from Hennessy et al 
formed by etching a planar microcavity with a suitable mask 



2007 Right panels: (a) photonic molecule 
(b) the measured (symbols) and calculated (lines) energies of 
the first two optical modes of the photonic molecule for various diameters are shown. Figure from |de Vasconcellos et aL\\201l\ 



quantum dot exciton transition. A promising route to 
overcome this difficulty using spatially localized emitters 
may involve pyramidal quantum dots that can be grown 



on demand at specific spatial positions (Carron et al. 



2011 ) and a subs equent fine tuning of t he exciton energy 
by Stark effect (IBennett et al.l 120101). Other researcli 



lines that are active with similar objectives involve im- 
purities in diamond crystals (Aharonovich et al. 2011) 
or single atoms ( [Dayan et al. 2008) located on the sur- 
face of a monolithic microcavity and strongly coupled to 
a high-Q cavity mode. 

An alternative strategy to realize a strong interaction 
condition is to use a quantum well as the active ma- 
terial (Verger et al. [2006 ): excitons are then strongly 
confined only along the growth axis, while the pattern- 
ing required for photon confinement affects their in-plane 
motion on a much longer micrometer scale. As the quan- 
tum well extends throughout the whole cavity and its 
thickness along the growth axis can be controlled at the 
level of a single atomic layer, the experimental problem 
of correctly positioning the emitter is completely elimi- 
nated and the electronic contribution to inhomogeneous 
broadening can be dramatically suppressed. 

While the nonlinear dynamics of polaritons in con- 
fined systems embedding quantum wells has been very 
recently addressed in ( ^Galbiati et al. 2011 1, it is still un- 
clear whether the nonlinearity stemming from collisions 
between delocalized polaritons as described by (31 1 is 



strong enough to penetrate the U/j ^ 1. A back-of-the- 
envelope approach based on the U oc 1 / scaling (here L 
is the characteristic length of lateral confinement) would 
suggest that very tightly confined polariton boxes down 
to the sub-micron scale should allow to reach the photon 
blockade regime (Verger et al. 2006). As the quality fac- 



tor of the resonator may also get quenched when the spa- 
tial size of the resonator is reduced to sub-micron scales, 
full numerical calculations of the photonic mode profile 
in realistic devices are required. Alternative strategies to 
overcome this difficulty may be to take advantage of the 
reinforced interaction potential between the hybrid di- 
rect/indirect excitons first studied in ( Cristofolini et aTj 
2012 ) or of the quantum interference mechanism reviewed 



in Sec lIX.B.41 



2. Superconducting quantum circuits and circuit QED systems 

Superconducting quantum circuits based on Josephson 
junctions arc another most promising playground where 
to investigate the physics of strongly correlated photons. 
Recent reviews of superconductor quantum circuits can 



be found in ( Schoelkopf and Girvin 2008 You and Nori 



2011). An example of architecture used in recent exper- 



iments is shown in the top panel of Fig |36| 

A simple theoretical description of these systems is ob- 
tained by representing the macroscopic electromagnetic 
variables such as currents and voltages (or charges and 
fluxes) as non-commuting quantum operators and then 
describing the evolution of the system in terms of an 
Hamiltonian. In the simplest case of a LC circuit, the 
Hamiltonian has the quadratic form of a harmonic oscil- 
lator. 

The addition of Josephson junctions to the circuit is 
essential for the emergence of quantum features and for 
quantum applications: thanks to its nonlinear current- 
voltage relationship, a Josephson junction introduces a 
non-dissipative nonlinearity into the circuit and produces 
a non-harmonic energy spectrum with non-equispaced 
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FIG. 36 Top panels: Sketch of a superconducting copla- 
nar waveguide resonator (green) interacting with a Joseph- 
son transmon qubit (left inset) that can be biased with 
large-bandwidth flux and charge gate lines (top inputs), (b) 
Schematic of resonator with two-level system generating a 
photon in mode a, emitting it into mode c onto a beam split- 
ter with modes d, e and f. Figure from (Bozyigit el al. 2011 1. 



Bottom panels: Experimental demonstration of resonant pho- 
ton blockade in a circuit-QED system, (a) Second-order corre- 
lation function measurements (dots) compared to theoretical 
calculations (lines) for different values of the drive amplitude. 
Same plot for a thermal incident field (b) and for a coherent 
drive (c). Figure from (Lang el al. 20111 



energy levels. As a result, the superconducting quan- 
tum circuit can be considered as an artificial atom with 
a discrete energy spectrum that depends on the specific 
connections and elements of the circuit. Under appro- 
priate conditions, the two lowest states of the quan- 
tum circuit are much closer in energy than all higher 
states, so that this effectively behaves as a two-level ar- 
tificial atom. By changing the type of circuit connec- 
tions and the relative values of the capacitances, induc- 
tances and Josephson energies, a wide variety of artifi- 
cial atoms can be obtained, among which the Cooper 
pair box(Nakamura tt al. 1999), the quantronium( Vion 



fluxonium (Manucharyan et al. 2009) 



et al.\ I2002D, the transmon" (|Koch et al.\ |2007[ ) or the 



By inserting a Josephson two-level artificial atom in a 
superconducting transmission line resonator, it is possi- 



ble to create a Jaynes-Cummings system where the ar- 
tificial atom is coupled to the bosonic field of the mi- 
crowave resonator. Given the analogy with usual cavity 



QED with real atoms (Berman et al. 1993 Raimond 



et al. 2001 ), this exciting field of research field has been 



called circuit QED (Blais et al. 2004 Wallraff et al. 



2004). With respect to standard cavity QED implemen 
tations using atoms in optical or microwave cavities, the 
superconducting circuit QED offer better perspectives of 
scalability: lithographic techniques analogous to the ones 
used in semiconductor technology can be in fact used to 
create large arrays of resonators and Josephson artificial 
atoms on a single solid-state chip. In addition, circuit 
QED offer the possibility to have very large Rabicou- 
plings comparable to the transition frequency (Devoret 



et al. 2007) and to efficiently control in-situ and in a 
dynamical way the properties of the artificial atom with 
external fields. This latter property was indeed at the 
heart of the recent observation of an analog dynamical 
Casimir effect (Wilson et al. 2011). 



Compared to semiconductors, the main practical limi- 
tation of superconducting quantum circuits is that they 
must operate in the microwave frequency range and at 
mK temperatures of a dilution fridge in order to have 
UbT <C <SC a, where A is the superconducting gap 
energy and uj is the resonator frequency (typically in the 
GHz range). On the other hand, the figures of merit of- 
fered by these circuit QED systems are quite spectacular 
as compared to previous quantum optical systems: for in- 
stance, a superconducting Jaynes-Cummings system can 
display a value of the 17/7 ratio between the vacuum Rabi 
frequency and the loss rate larger than 10^: a comparable 
value 3/7 ~ 300 was obtained so far only with Rydber g 
atoms in superconducting cavities ( [Raimond et al.\ 2001 1, 
but these systems are hardly scalable to many-cavity ge- 
ometries. This remarkable large value of (7/7 was cru- 
cial in the recent experiments that have used Josephson 
atoms with tunable coupling to syntetie arbitrary quan- 
tum states of the coupled qubit-cavity field system and 
then perform a quantum tomography measurement of the 
state pofheinz et a/.[|2009| ). 

While single photon detectors in the microwave range 



are still in the exploratory phase (Chen et al. 2010 



Romero et al. 20091, two-channel heterodyne schemes 
using linear amplitude detectors have allowed for accu- 
rate measurements of the quadratures of the microwave 



field down to the quantum limit (Bozyigit et al. 2011 



Gabelli et al. 2004 Lang et al. 2011): an example of 



application of this technique is illustrated in the bottom 
panel of Fig |36[ where we show how a strong photon anti- 
bunching and a marked sub-poissonian statistics can be 
obtained via a photon blockade mechanism under a co- 



herent drive (Lang et al. 2011) 



The possibility of coupling in a controllable way sev- 
eral resonators has been demonstrated in (Mariantoni 



et al. 2008 2011 Wang et al. 20111, paving the way 



to the study of arrays of resonators in the near future. 
These impressive experimental advances, together with 
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the proposed implementations of artificial gauge fields 
for microwave photons (Koch et al. 2010 Nunnenkamp 



et al. 2011 1 puts superconducting circuits among the 



most promising candidates for the study of many-body 
physics in photon gases. 



X. CONCLUSIONS AND PERSPECTIVES 

In this review we have summarized recent develop- 
ments in the theoretical and experimental study of gases 
of interacting photons in solid-state systems. An effec- 
tive photon mass appears as a result of spatial confine- 
ment along the growth axis, while the nonlinear optical 
susceptibility of the material medium induces sizable bi- 
nary interactions between photons. Strong coupling of 
the photon with some long-lived electronic excitation in 
the medium is a succesful strategy to reinforce binary in- 
teraction between the dressed bosonic particles that arise 
from the mixing of light with the matter excitation, the 
so-called polaritons. Recent experiments demonstrating 
Bose-Einstein condensation and superfluidity effects in 
luminous gases were indeed performed using polariton 
gases in semiconductor planar microcavities embedding 
quantum wells. 

In contrast to standard many-body systems such as 
liquid Helium, electron liquids in solid-state systems and 
ultracold atomic gases, polaritons have a finite lifetime as 
a consequence of radiative and non-radiative dissipative 
processes and some external pump is required to replen- 
ish the gas and compensate for the particle loss. As a 
result, the non-equilibrium steady state of the system 
is no longer determined by a standard thermodynami- 
cal equilibrium condition, rather by a dynamical balance 
of driving and dissipation. On one hand, the driven- 
dissipative nature of the photon fluid is responsible for 
novel phenomena in the collective dynamics; on the other 
hand it provides (ironically) a powerful tool for real-time 
diagnostics of the many-body system just by looking at 
the emitted light. 

After a brief survey of the historical development of 
the concept of photon fluids, we started our presenta- 
tion with a general theoretical summary of the quantum 
field description of photons and polaritons in solid-state 
devices: even if our focus was concentrated on planar ge- 
ometries, the concepts are readily transferred to other 
geometries, including one-dimensional polariton wires, 
harmonic traps, and even periodic lattices. We have 
then illustrated the generalized non-equilibrium Gross- 
Pitaevskii equation that can be used to describe the 
steady-state and the collective dynamics of a dilute pho- 
ton gas at the mean-field level. A solution of this equa- 
tion is provided for the simplest geometries in the differ- 
ent pumping regimes. This general theoretical results are 
then used to explain recent experimental observations of 
the condensate shape and the momentum distribution. 

Linearization of the generalized Gross-Pitaevskii equa- 
tion around the steady state provides the spectrum of 



elementary excitations of the fiuid. Based on the general 
features of this spectrum, we have reviewed the recent 
demonstration of superfluid hydrodynamic behaviors in 
polariton fluids. Under a coherent pumping, a variety of 
regimes was observed depending on the flow speed com- 
pared to the Landau critical velocity, from superfluid flow 
to the Bogoliubov-Cerenkov emission of phonons and sin- 
gle particle excitations in the wake of a weak defect; for 
large and strong defects, the strong local perturbation 
of the fluid has been shown to lead to the nucleation 
of topological excitations such as vortices and dark soli- 
tons. Intriguing and so far unexplained experimental ob- 
servations of superfluid behaviors under a OPO pumping 
regime have been critically discussed. 

In the last part of the review, we have presented an 
overview of the emerging field of strongly correlated pho- 
ton gases. For strong photon-photon interactions, the 
mean-field description based on a Gross-Pitaevskii equa- 
tion for the macroscopic wavefunction breaks down and 
the many-photon wavefunction starts displaying peculiar 
quantum correlations. Exception made for very simple 
single-mode cavities for which effective photon blockade 



has been already observed in a number of systems ( Birn- 
baum "aH |2005| [Faraon et al. 2008 Lang et al. 



2011 



Reinhard et al. 2011), the first experimental claim of 



strongly interacting photons in spatially extended geome- 
try has appeared very recently using an atomic gas in the 
Rydberg-EIT configuration (Peyronel et al. 2012 1. In 



the meanwhile, an intense activity is devoted to the theo- 
retical study of systems where quantum states of photon 
matter can be generated, for instance Mott insulators, 
Tonks- Girardeau gases, Laughlin states of quantum Hall 
physics. Also in this case, the driven-dissipative nature of 
the photon gas is responsible for a wealth of unexplored 
features due to the interplay of non-equilibrium statisti- 
cal mechanics, quantum optics, many-body physics, pos- 
sibly with long-term applications for all-optical quantum 
information processing. 

In spite of the great achievements that the research on 
photon gases has experienced in the last few years, an 
even larger number of questions are still awaiting exper- 
imental and/or theoretical answer. From the experimen- 
tal point of view, the big challenge is in our opinion the 
development of devices in the visible and/or in the mi- 
crowave range of frequencies where strong photon-photon 
interactions can be associated with a non-trivial spatial 
dynamics, either in free space geometries or in lattice 
ones. A specific discussion of two systems that we be- 
lieve most promising in this direction is given in the last 
section of the article, namely semiconductor micro- and 
nano-resonators for visible or infrared light and circuit- 
QED systems for microwaves. 

Important questions are still open also in the dilute 
gas regime. In addition to a complete explanation of ob- 
servations, more experimental work to fully characterize 
the elementary excitation spectrum of a dilute photon 
gas in the different pumping regimes is required and is 
expected to shine light on novel aspects of superfiuidity. 
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Of course, all this discussion is not restricted to systems 
in the strong coupling regime where the elementary ex- 
citations have a polaritonic nature: for instance, verti- 
cal cavity surface emitting laser devices are promising 
alternative candidates for the study of the collective dy- 
namics of the photon gas. In turn, a better understand- 
ing of these many-body features of the non-equilibrium 
Bose-Einstein transition will shine new light on the laser 
threshold of spatially extended devices such as VCSELs: 
in particular, very little appear to be known yet on the 
critical fluctuations in the vicinity of the transition point. 

Another research direction that is expected to have 
a tremendous impact on both fundamental science and 
on applications is the one on the generation of artifi- 
cial gauge fields for photons: on one hand, quantum me- 
chanics of particles in the presence of gauge fields has 
been predicted to show a number of novel fascinating 
effects. On the other hand, we expect that their full 
exploitation in practical photonic devices will allow to 
take advantage of topological features to elaborate opti- 
cal information with unprecedented robustness towards 
experimental imperfections. Inclusion of optical nonlin- 
earities in the arena is an almost unexplored field that 
may considerably enlarge the spectrum of possibilities 
with a number of novel unexpected effects and possibly a 
bright future for the physics of luminous quantum fluids. 
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